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ABSTRACT 

Heat  transfer  tests  were  carried  out  using  a  small  tube  bundle  of 
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instrumenting  five  tubes  within  the  bundle,  both  the  convective  and 
nucleate  boiling  regions  were  studied  in  detail,  with  emphasis  on  the 
'bundle  effect'  (ie.  the  effect  of  the  lower  tubes  in  operation  on  the 
upper  tubes  within  the  bundle) .  In  addition,  the  influence  of  increased 
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I.  INTRODUCTION 

A.    BACKGROUND 

One  of  today's  major  environmental  concerns  is  the  depletion  of  the 
earth's  protective  ozone  layer.  In  1987,  an  international  conference  was 
held  in  Montreal,  Canada  to  address  the  problems  caused  by 
Chlorofluorocarbons  (CFCs)  to  the  earth's  ozone  layer.  CFCs  are  manmade 
chemicals  of  chlorine,  fluorine,  and  carbon  and  are  unique  in  that  they 
have  a  combination  of  desirable  properties:  low  in  toxicity,  non- 
flammable, non-corrosive,  non-explosive,  extremely  stable  and  compatible 
with  many  other  materials.  This  extreme  stability  is  what  causes  problems 
to  the  ozone  layer  due  to  the  fact  that  CFCs  only  break  down  in  the  upper 
atmosphere  when  subjected  to  intense  ultraviolet  radiation.  This  break 
down  produces  chlorine  which  has  been  linked  to  the  depletion  of  the 
earth's  ozone  layer. 

In  September  1987,  24  nations  representing  the  United  Nations 
Environment  Program  (UNEP)  met  and  signed  the  Montreal  Protocol.  They 
discussed  the  substances  that  deplete  the  ozone  layer  [Ref.l]  and  called 
for  a  near-term  freeze  on  the  production  and  consumption  of  these 
substances.  The  agreement  required  production  of  these  chemicals  to  be 
cut  back  to  1986  levels  followed  by  a  two-phased  reduction  culminating  in 
cutbacks  of  50%  by  mid-1998;  this  came  into  effect  on  July  1,  1989.  In 
1990,  a  progress  meeting  was  held  in  London  where  UNEP  delegates  agreed  to 
completely  phase  out  all  CFCs  by  the  year  2000  [Ref.  2].   In  the  spring  of 


1992,  President  Bush  pushed  up  the  complete  phase  out  od  CFC's  by  the  year 
1995. 

The  U.S.  Navy  uses  a  number  of  different  CFCs  (designated  by  "R"  for 
refigerants)  for  various  refrigeration  and  air  conditioning  (AC)  needs. 
Presently,  the  U.S.  Navy  has  approximately  1850  shipboard  AC  plants  using 
both  R-12  (in  reciprocating  compressor)  and  R-114  (in  centrifugal 
compressor)  plants.  To  comply  with  the  Montreal  Protocol  and  U.S. 
legislation,  the  Mechanical  Systems  Branch  (Code  2772)  at  the  Naval 
Surface  Warfare  Center  (NSWC)  is  pursuing  research  in  the  elimination  of 
shipboard  use  of  CFCs.  As  mentioned  by  Chilman  [Ref.  3]  this  research  is 
to  be  completed  in  three  phases: 

1.  To  identify  in  the  short  term  suitable  alternative  ozone— safe 
chemicals  to  replace  R-114  and  R-12.  To  accomplish  this  task,  the  heat 
transfer  characteristics  must  be  similar  to  the  existing  refrigerants  in 
place  and  hence  the  need  for  a  database  exists  for  current  refrigerants 
(R-114  and  R-12)  so  that  they  can  be  compared  to  the  new  proposed 
refrigerants  (HFC-124  and  HFC-134A  respectively) . 

2.  In  the  longer  term,  to  research,  develop,  and  test  substitute 
chemicals  and  alternative  technologies  to  replace  existing  CFC  uses. 

3.  To  implement  new  cooling  system  technologies  into  the  fleet  which 
do  not  depend  on  CFCs  or  their  replacement. 

This  thesis  is  a  continuation  of  the  previous  work  at  NPS  and 
supplements  NSWC's  research  on  alternatives  to  CFCs  by  establishing 
baseline  nucleate  pool  boiling  data  of  pure  R-114  and  R— 114/oil  mixtures 
from  a  small  bundle  of  enhanced  tubes  representing  a  section  of  a  flooded 
evaporator.   Emphasis  is  placed  on  the  natural  convection  and  boiling 


regions,  hystersis  phenomena,  and  analysis  of  various  oil  concentrations 
on  the  overall  heat  transfer  performance. 

B.    OBJECTIVES 

The  objectives  of  this  thesis  are  as  follows: 

1.  Understand  in  greater  detail  both  the  convection  and  nucleate  pool 
boiling  phenomena  and  hysteresis  effects  within  a  small  Turbo— B  tube 
bundle . 

2.  Obtain  data  using  a  Turbo— B  tube  bundle  for  increasing  and 
decreasing  heat  flux  for  R-114/oil  mixture  with  oil  concentrations  of  0, 
1,  2,  3,  6,  and  10  percent. 

3.  Compare  data  with  earlier  studies  at  the  Naval  Postgraduate  School 
using  R-114/oil  mixtures  for  other  enhanced  tube  surfaces. 


II.  LITERATURE  SURVEY 

A.    SMOOTH  TUBE  BUNDLES 

In  recent  years  significant  progress  has  been  made  in  understanding 
nucleate  boiling  heat  transfer  phenomena  on  the  shell  side  of  flooded 
evaporators.  Extensive  work  on  smooth  tube  bundles  has  been  reported  by 
Cornwell  (Leong  and  Cornwell  [Ref.  4],  Cornwell  et  al .  [Ref.  5],  Cornwell 
and  Scoones  [Ref.  6],  Cornwell  [Ref.  7]).  Cornwell  and  Schuller  [Ref.  8] 
conducted  a  photographic  study  of  boiling  R-113  in  a  smooth  tube  bundle  at 
one  atmosphere.  One  of  their  conclusions  was  that  bubbles  leaving  the 
lower  tubes  in  the  bundle  impacted  and  caused  a  sliding  motion  around  the 
upper  tubes.  Cornwell  and  Schuller  observed  the  two— phase  flow  patterns 
and  deduced  that  sliding  bubbles  from  lower  tubes  on  upper  tubes  could 
account  for  significant  heat  transfer  in  the  top  part  of  the  bundle. 
Cornwell  [Ref.  7]  later  found  that  in  the  nucleate  boiling  region,  sliding 
bubbles  and  liquid  forced  convection  could  account  for  all  the  heat 
transfer  in  the  top  of  the  bundle. 

The  influence  of  tube  position  within  a  bundle  of  smooth  tubes  using 
R-ll,  R-12,  R-22  and  R-113  has  been  studied  extensively  by  Wallner  [Ref. 
8],  Fujita  et  al .  [Ref.  9],  Chan  and  Shourki  [Ref.  10],  Rebrov  et  al . 
[Ref.  11],  and  Marto  and  Anderson  [Ref.  12].  Using  both  in-line  and 
staggered  tube  arrangements  with  various  tube  pitch-to-diameter  ratios 
between  1.2  and  2.0,  their  work  verified  that  the  influence  of  the  lower 
tubes  in  a  bundle  can  significantly  increase  the  heat  transfer  performance 


of  upper  tubes  at  low  heat  fluxes  due  to  two-phase  convective  effects.  At 
high  heat  fluxes  (typically  >  50  kW/m2)  in  the  fully  developed  nucleate 
boiling  region,  the  data  for  all  the  tubes  merged  onto  a  single  curve. 
This  is  representative  of  a  single  smooth  tube  and  shows  that  there  is  no 
'bundle  effect'  (ie.  no  improvement  over  a  single  tube  under  similar 
conditions)  in  the  high  heat  flux  (nucleate  boiling)  region. 

Chan  and  Shoukri  [Ref.  10]  studied  the  boiling  characteristics  of  a 
smooth  in-line  tube  bundle  in  R-113.  They  concluded  that  at  lower  heat 
fluxes,  the  heat  transfer  process  is  strongly  influenced  by  two-phase 
convection  effects,  resulting  in  higher  heat  transfer  coefficients  on  the 
upper  tubes.  At  high  heat  fluxes,  however,  they  found  that  the  dominant 
mode  of  heat  transfer  was  nucleate  boiling  from  the  upper  tubes  and  that 
convective  effects  from  below  were  insignificant.  At  these  high  fluxes, 
the  bundle  performance  was  similar  to  the  trends  of  a  single  tube  in  a 
single  tube  apparatus.  Fujita  et  al .  [Ref.  9]  also  found  that  the  heat 
transfer  at  low  heat  fluxes  using  a  smooth  tube  bundle  in  R— 113  was 
enhanced  by  convection  induced  by  rising  bubbles  (ie.  a  steady  increase  in 
performance  of  the  upper  tubes  as  additional  lower  tubes  were  activated) . 
They  attributed  this  enhancement  to  the  "positive  bundle  effect".  A  high 
heat  fluxes,  this  enhancement  disappeared. 

Anderson  [Ref.  13]  found  similar  effects  as  above  for  a  smooth  tube 
bundle  in  pure  R-114.  Furthermore,  he  reported  that  the  presence  of  up  to 
3%  oil  (by  mass)  actually  improved  the  heat  transfer  performance.  This  is 
similar  to  data  reported  for  a  single  smooth  tube  by  Wanniarachchi  et  al . 
[Ref.  14]  .  Furthermore,  at  an  oil  concentration  of  10%,  only  a  slight 
degradation  in  the  heat  transfer  (compared  to  pure  R— 114)  was  found.   He 


obtained  a  maximum  heat  transfer  performance  for  the  bundle  at  an  oil 

concentration  of  2%. 

S 

B.    ENHANCED  TUBE  BUNDLES 

Much  less  work  has  been  done  on  enhanced  tubes  (Enhanced  means  any 
surface  that  is  not  smooth).  However,  the  effects  at  high  and  low  heat 
fluxes  mentioned  above  are  similar  to  those  obtained  for  finned  tube 
bundles  by  Yilmaz  and  Palen  [Ref .  15] ,  Muller  [Ref .  16] ,  and  Hahne  and 
Muller  [Ref.  17].  Stephen  and  Mitrovic  [Ref.  18]  looked  at  R-12  and  R-114 
boiling  from  a  GEWA-T  tube  bundle.  Apart  from  the  magnitude  of  the  heat 
transfer  coefficient  varying  with  fluid,  the  trends  are  very  similar  to 
those  mentioned  above  for  smooth  and  finned  tube  bundles. 

For  porous  coated  surfaces,  Czikk  et  al .  [Ref.  19]  found  no  'bundle 
effect'  over  a  wide  range  of  heat  flux  (1—100  kW/m2)  and  the  bundle  data 
agreed  closely  with  single  tube  data.  Arai  et  al.  [Ref.  20]  found  that 
the  'bundle  effect'  for  a  Thermoexcel  tube  bundle  was  smaller  than  that 
found  for  a  smooth  or  finned  tube  bundle.  However  as  before,  any  'bundle 
effect'  was  eliminated  at  high  heat  fluxes  where  the  data  for  all  the 
tubes  agreed  closely  with  single  Thermoexcel  tube  results.  These  effects 
are  similar  to  those  found  by  Czikk  et  al.  [Ref.  19]  for  the  porous 
coated. 

Chilman  [Ref.  3]  reported  experiments  with  a  Turbo-B  tube  bundle  in 
pure  R-113,  conducting  both  increasing  and  decreasing  heat  flux  tests.  He 
concluded  that  in  the  natural  convection  region,  heated  lower  tubes  do  not 
have  any  influence  on  the  heat  transfer  from  upper  tubes.   Also,  Chilman 


reported  the  presence  of  heated  lower  tubes  within  a  bundle  reduced  the 
incipient  boiling  point. 

Stephan  and  Mitrovic  [Ref.  18]  reported  the  influence  of  oil  on  the 
boiling  heat-transfer  coefficient  of  R-12  using  a  T-shaped  finned  tube 
(Gewa-T)  bundle.  They  reported  that  the  ratio  of  oil  to  no  oil  heat- 
transfer  coefficients  decreased  with  the  mass  fraction  of  oil  for  all 
except  the  highest  heat  flux  (22  kW/m2)  where  an  increase  in  heat  transfer 
was  noted  for  oil  concentrations  between  1  and  6%  .  They  concluded  that 
the  influence  of  oil  on  heat  transfer  was  mainly  due  to  the  thermal 
properties  of  the  specific  oil  used  in  the  experiments  and  its  interaction 
with  the  refrigerant. 

Schlager  et  al.  [Ref.  21]  summarized  the  influence  of  oil  on 
refrigerant  in  pool  boiling.  They  stated  that  under  certain  conditions 
(typically  low  pressure  and  high  heat  flux) ,  the  heat  transfer  coefficient 
increased  at  low  oil  concentration.  Stephan  [Ref.  22]  first  pointed  this 
out  and  attributed  the  phenomenon  to  foaming.  Burkhardt  and  Hahne  [Ref. 
23]  for  a  finned  tube  bundle  found  that  the  maximum  heat  transfer 
coefficient,  which  was  10%  to  15%  above  the  oil-free  value,  occurred  at  a 
concentration  of  about  4%. 

Heimbach  [Ref.  24]  conducted  experiments  with  R— 12/oil  mixtures  on  a 
finned  tube  bundle.  He  reported  that  the  presence  of  up  to  2%  oil,  did 
not  affect  the  heat  transfer  performance  significantly.  However  at  higher 
concentrations  (3%  to  7%),  an  increase  in  the  heat  transfer  was  observed. 
He  also  attributed  this  to  foaming  and  postulated  that  changes  in  the 
properties  of  the  mixture  might  facilitate  the  formation  of  bubbles. 


Anderson  [Ref .  13]  and  Akcasayar  [Ref .  25]  conducted  experiments  with 
finned  (19  fpi)  and  High  Flux  (porous  coated)  tube  bundles  in  pure  R-114 
and  R-114/oil  mixtures  in  the  same  apparatus.  Anderson  reported  a  maximum 
heat  transfer  performance  at  an  oil  concentration  of  3%  for  the  finned 
tube  bundle.  Akcasayar  also  reported  that  the  finned  tube  bundle 
performance  increased  1.65  times  with  3%  oil  concentration  (compared  with 
pure  R— 114)  at  the  maximum  heat  flux  level.  For  6%  and  10%  oil 
concentrations,  the  performance  of  the  bundle,  when  compared  with  lower 
oil  concentrations,  decreased.  This  was  especially  significant  at  a  10% 
oil  concentration.  Compared  with  the  finned  tube  bundle,  the  High  Flux 
tube  bundle  had  a  1.5  times  better  heat  transfer  performance  at  a  heat 
flux  of  30  kW/m2  in  pure  R-114.  However,  these  performance  ratios 
decreased  with  increased  oil  concentrations  such  that  the  finned  bundle 
outperformed  the  High  Flux  bundle  at  6%  oil  concentration.  This  was 
especially  true  at  the  highest  heat  fluxes  where  the  High  Flux  bundle 
performance  was  not  much  better  than  a  smooth  tube  bundle  of  similar  size. 


III.  EXPERIMENTAL  APPARATUS 

A.    TEST  APPARATUS  OVERVIEW 

The  experimental  apparatus  including  the  auxiliary  equipment  and  the 
evaporator/condenser  is  shown  in  Figure  1.  The  following  is  only  a 
general  description  of  the  whole  experimental  apparatus.  A  more  detailed 
look  at  the  condenser  and  evaporator  is  provided  in  section  C.  Further 
information  about  the  apparatus  is  provided  by  Murphy  [Ref.  26]  and 
Anderson  [Ref.  13]. 

The  experimental  apparatus  is  essentially  made  up  of  three  closed 
loops.  The  first  loop  consists  of  an  8  ton  refrigeration  unit  located 
outside  the  laboratory  which  is  used  to  cool  an  ethylene  glycol/water 
mixture.  The  second  loop  is  the  ethylene  glycol/water  mixture  flowing 
through  the  condenser.  This  mixture  is  contained  within  a  large  sump 
within  the  laboratory.  The  flow  rate  through  the  condenser  is  delivered 
by  two  pumps  which  can  be  operated  independently  or  together;  this  coolant 
mixture  condenses  the  refrigerant  vapor  in  the  condenser  and  maintains 
system  pressure  and  temperature.  Pump  number  one  provides  coolant  flow 
through  the  four  test  condenser  tubes  as  well  as  to  one  of  the  auxiliary 
condenser  coils  (bottom  coil) .  Pump  number  two  provides  coolant  through 
a  manifold  which  distributes  the  coolant  to  the  remaining  four  auxiliary 
condenser  coils  within  the  condenser.  The  third  loop  is  the  evaporator 
and  condenser  itself  designed  for  reflux  operation.   The  vapor  generated 


in  the  evaporator  flows  upward  and  condenses  in  the  condenser;  the 
condensate  then  returns  to  the  evaporator  via  gravity. 

B.    AUXILIARY  EQUIPMENT 

1.   28  kW  Refrigeration  Unit 

This  unit  was  used  to  cool  a  1.8  m3  reservoir  sump  of  an 
ethylene  glycol/water  mixture  (coolant)  to  the  desired  temperature  needed 
to  condense  the  refrigerant  vapor.  For  R-114,  the  temperature  control  was 
set  to  maintain  the  sump  at  -15  °C.  The  refrigeration  unit  had  a  cooling 
capacity  of  28  kW  (8  tons). 

2  .   Ethylene  Glycol/Water  Mixture 

The  coolant  used  was  a  54%  (by  weight)  ethylene  glycol/water 
mixture.  This  refrigerated  mixture  was  used  to  control  the  system 
pressure  and  temperature  by  circulating  coolant  through  the  auxiliary 
condenser  coils  and/or  condenser  test  tubes  at  different  flow  rates. 

3 .  Pumps 

Two  pumps  were  available  to  circulate  the  coolant  from  the  sump 
through  the  condenser.  Pump  number  one  fed  four  test  condenser  tubes  and 
one  of  the  auxiliary  condenser  coils.  Pump  number  two  fed  the  other  four 
auxiliary  condenser  coils.  Pump  number  one  was  the  primary  pump  used  at 
low  heat  fluxes;  pump  number  two  was  started  (as  necessary)  at  high  heat 
fluxes  to  maintain  the  desired  saturation  pressure  in  the  evaporator. 

4.  Flowmeters 

Five  calibrated  float-type  flowmeters,  connected  to  pump  number 
one,  were  used  to  measure  the  flow  rates  passing  through  the  four  test 
condenser  tubes  and  one  auxiliary  condenser  coil.    One  additional 
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flowmeter  (connected  to  pump  number  two)  was  used  to  measure  the  total 
flow  to  the  remaining  four  auxiliary  condenser  coils.  Each  of  the  five 
auxiliary  condenser  coils  had  a  globe  valve  to  regulate  (or  shut  off) 
flow  as  desired.  For  the  four  test  condenser  tubes,  the  coolant  flow  was 
regulated  by  a  flowmeter  valve. 

C .    EVAPORATOR/CONDENSER 

An  overall  view  of  the  evaporator  and  condenser  is  shown  in  Figure  2 . 
The  evaporator  was  designed  to  simulate  a  small  portion  of  a  refrigerant 
flooded  evaporator.  Front  and  side  views  of  the  evaporator  are  shown  in 
Figures  3  and  4.  It  was  fabricated  from  stainless  steel  plate  and  formed 
into  a  short  cylinder,  610  mm  in  diameter  and  241  mm  long.  Electrically- 
heated  tubes  were  cantilever-mounted  from  the  back  wall  of  the  evaporator 
to  permit  viewing  along  the  axis  of  the  tubes  through  the  lower  of  two 
viewing  windows.  A  plexiglas  plate  was  attached  to  the  front  of  the  tube 
bundle  to  ensure  tube  alignment  during  experiments.  Each  viewing  window 
had  a  layer  of  glass  and  plexiglas,  the  glass  being  used  on  the 
refrigerant  side  in  order  to  prevent  surface  cracking  of  the  plexiglas. 
The  plexiglas  gave  the  glass  added  strength  and  served  as  a  safety  barrier 
in  case  the  glass  broke  during  pressurization. 

The  electric  power  can  be  applied  separately  to  each  set  of  heaters 
using  a  STACO  240  V,  23.5  kVA  rheostat  controller  shown  in  Figure  5. 
Also,  the  desired  number  of  instrumented  tubes,  active  tubes,  simulation 
or  auxiliary  heaters  can  be  individually  activated  by  using  circuit 
breakers.  The  five  simulation  heaters,  each  with  a  maximum  rating  of  4 
kW,  were  mounted  below  the  test  bundle  in  order  to  simulate  15  additional 
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tube  rows  in  a  larger  bundle  and  to  provide  an  inlet  vapor  quality  into 
the  bottom  of  the  test  bundle  as  suggested  by  Webb  [Ref .  27] .  The  liquid 
pool  was  maintained  at  2.2  °C  (corresponding  to  a  saturation  pressure  of 
1  atmosphere)  by  passing  coolant  through  the  condenser. 

Figure  6  is  a  schematic  sectional  view  of  the  evaporator  that  shows 
the  four  kinds  of  heated  tubes  installed  in  the  evaporator.  These  were 
instrumented,  active,  auxiliary,  and  simulation.  For  this  study  only  the 
instrumented,  active,  and  simulation  heaters  were  used;  the  auxiliary 
heaters  are  needed  for  experiments  either  at  higher  pressures  or  for  other 
refrigerants  which  have  a  higher  normal  boiling  point  (such  as  R-113) . 
Table  1  gives  the  power  rating  for  these  heaters  and  the  number  used  in 
the  evaporator. 

The  tube  bundle  itself  consists  of  instrumented,  active,  and  dummy 
tubes.  The  location  of  each  tube  is  represented  by  the  respective  letter 
I,  A,  and  D  as  shown  in  Figure  6.  The  test  bundle  consists  of  two  types  of 
heated  tubes:  12  active  tubes  (marked  "A")  which  contained  1  kW  cartridge 
heaters,  and  5  instrumented  tubes  (marked  "I")  which,  in  addition  to  the 
1  kW  cartridge  heaters,  contained  six  wall  thermocouples  each. 

In  measuring  boiling  heat  transfer  coefficients,  great  care  must  be 
exercised  with  the  cartridge  heater  and  temperature  measuring 
instrumentation  to  ensure  good  accuracy.  The  instrumented  test  tubes  were 
fabricated  using  the  same  method  as  that  used  by  Hahne  and  Muller  [Ref. 
17]  and  Wanniarachchi  et  al .  [Ref.  14].  The  exact  procedure  can  be  found 
in  Eraydin  [Ref.  28].  Figure  7  is  a  cross-sectional  sketch  of  an 
instrumented  tube,  showing  the  construction  details  and  the  location  of 
the  wall  thermocouples.   The  thermocouples  were  embedded  in  the  wall  at 
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different  circumferential  and  longitudinal  positions  along  the  heated 
section  of  the  tube  shown  in  Figure  7. 

The  five  instrumented  tubes  were  located  along  the  centerline  of  the 
tube  bundle,  forming  a  vertical  in-line  column.  All  the  instrumented  and 
active  tubes  were  Turbo-B  tubes  made  by  Wolverine  Tube  Co.  (see  Section 
E) .  These  tubes  were  nominally  15.9  mm  in  outside  diameter  and  were 
arranged  in  an  equilateral  triangular  pitch  (ie.  centerline— to— centerline 
spacing)  of  19.1  mm,  giving  a  pitch-to-diameter  ratio  of  1.35.  The 
thickness  of  the  Turbo-B  enhancement  was  0.85  mm  giving  a  diameter  to  the 
base  of  the  enhancement  of  14.2  mm. 

The  bundle  also  contained  a  number  of  unheated  dummy  smooth  tubes 
(marked  "D" )  that  were  used  to  guide  the  two  phase  mixture  through  the 
bundle.  The  dummy  tubes  were  made  from  commercially  available  15.9  mm  OD 
smooth  copper  tubing.  Two  vertical  baffle  plates  made  of  aluminum  were 
used  on  either  side  of  the  bundle  to  restrict  circulation  into  and  out  of 
the  bundle  at  the  sides.  A  dummy  rack  (Figure  8)  consisting  of  12  solid 
rods  made  of  aluminum  (15.9  mm  OD  and  spaced  19.1  mm  from  centerline-to- 
centerline)  was  placed  below  the  tube  bundle.  This  rack  had  a  triangular 
pitch  arrangement  with  vapor  retaining  plates  on  the  sides  and  was 
designed  for  two  purposes:  to  collect  all  rising  two  phase  flow  generated 
by  the  simulation  heaters  and  direct  it  into  the  test  bundle  and  to 
simulate  vapor  passing  through  a  large  bundle  before  reaching  the 
instrumented  tubes.  A  small  open  space  (approximately  5  mm  in  height)  was 
left  between  the  bundle  and  dummy  tube  rack.  This  space  allowed  some 
refrigerant  from  below  to  enter  the  bundle  and  replace  the  vapor  being 
generated  in  the  bundle.  However,  there  was  also  a  space  below  the  dummy 
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rack  that  allowed  the  majority  of  the  circulation  to  occur.  Thus, 
liquid/vapor  circulation  was  vertically  upward  over  the  five  instrumented 
test  tubes  with  no  net  horizontal  component.  Most  of  the  liquid-vapor 
mixture  after  passing  through  the  bundle  was  separated  when  it  reached  the 
pool  surface.  However,  due  to  the  strong  circulation  patterns  set  up 
within  the  liquid  pool,  some  vapor  bubbles  remained  trapped  in  the  liquid 
and  circulated  around  the  pool. 

The  condenser  included  four  test  tubes  (each  1.219  m  in  length  and 
15.9  mm  OD)  in  a  vertical  in-line  column  and  five  auxiliary  copper  coils. 
For  the  boiling  experiments,  these  tubes  were  used  to  regulate  the 
pressure  and  temperature  in  the  evaporator.  The  condenser  was  designed  to 
permit  independent  condensation  studies  of  small  in-line  tube  bundles, 
using  the  evaporator  as  a  source  of  vapor.  Details  of  the  condenser  can 
be  found  in  Mazzone  [Ref .  29] . 

D.    DATA  ACQUISITION  SYSTEM/INSTRUMENTATION 

As  described  by  Akcasayar  [Ref.  25]  ,  a  Hewlett  Packard  HP-3497A  Data 
Acquisition  System,  HP-9125  computer  and  HP-701  printer  were  used  for  data 
acquisition,  data  reduction  and  data  printing  respectively.  Although  an 
HP-9826  computer  and  HP-7470A  plotter  can  be  used  for  final  graph 
printing,  a  Macintosh  computer  (using  Criketgraph)  was  utilized.  As 
described  by  Anderson  [Ref.  13]  ,  type-T  copper-constantan  thermocouples 
were  used  to  make  temperature  measurements  on  the  HP— 3497A  using  a  relay 
multiplexer  assembly  equipped  with  thermocouple  compensation.  A  20- 
channel  relay  multiplexer  card  was  used  to  measure  the  voltage  output  from 
both  the  voltage  and  current  sensors.  The  voltage  measurements  were  taken 
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from  separate  sensors  that  measured  the  voltage  going  to  the  tube  bundle, 
the  simulation  heaters  and  the  auxiliary  heaters.  The  total  current 
going  through  the  auxiliary  and  simulation  heaters  was  measured  using  an 
American  Aerospace  Control  (ACC)  current  sensor.  The  current  to  each 
instrumented  tube  heater  was  measured  using  five  identical  current 
sensors.  The  voltage  supplied  to  the  other  active  tubes  was  also 
measured,  but  the  current  of  each  active  tube  was  not.  Instead,  the  total 
current  for  a  pair  of  active  tubes  was  measured,  and  this  was  sufficient 
since  these  tubes  each  had  the  same  power  output  (1000  W)  as  the 
instrumented  test  tube  heaters  and  there  was  no  apparent  reason  to  monitor 
each  active  tube  individually.  Computer  channel  assignments  for  data 
acquisition  and  array  assignments  are  given  in  Table  2. 

E.    GEOMETRY  OF  TURBO -B  TUBE 

The  Turbo-B  tube,  manufactured  by  Wolverine  Tube  Inc.,  contained  an 
enhanced  surface  geometry.  The  exterior  boiling  enhancement  is  made  by 
raising  low  integral  fins,  cutting  diagonally  across  these  fins,  and  then 
rolling  the  fins  to  compress  them  to  form  mushroom-like  pedestals  [Ref. 
30].  This  process  forms  numerous  re-entrant  passageways.  Figure  9  shows 
the  surface  of  the  tube  at  25  times  its  actual  size.  The  tube  is 
currently  available  in  copper,  cupro-nickel,  and  low  carbon  steel. 

The  relative  dimensions  of  the  tube  used  in  this  study  are  as  follows: 
Tube  material  -  Copper 
Nominal  Outside  diameter  =  15.9  mm 
Enhanced  surface  length  =  203.2  mm 
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Thickness  of  Enhancement  -  0.85  mm 
Diameter  to  Base  of  Enhancement  =  14.2  mm 
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Figure   1.    Schematic  View  of  the   Experimental  Apparatus 
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Figure  2 .  Evaporator/Condenser  Schematic 
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Figure  3.  Front  View  of  Evaporator 


19 


b 

E 

E 


Z 

X 

o 

u 

V) 

o 

t/> 

_1 

UJ 

CO 

DC 

Q. 

K 

2 

oc 

O 

o 

o 

a 

a 

ouj 

3 

zb 

to 

?3 

UJ 

•   q. 

-1 

o 

— 

o 

K\ 

_) 

3 

CD 

o 

O 

zce 

l 
UJ 

ItOin 

(D 

E 

ujQ-o 

3 

fe 

►"9:2 

K 

<£ 

23- 

—  <O0T 

K 


E 
E 

«-> 
O 


E 
E 
ro 
O 


to 

t- 
< 

a 


0> 


5^S 

E 

E 

.—  OD  _l 

<"oa 

N. 

f»  UJ  — 

c\i 

ox- 

_o~- 

eS| 

E       E 

oi     m 

—      <\j 

o 

UJ? 

<•- 

_iO 
u.  2 


E  E 
E  E 

CDtD 
NO 

—  tr 


Figure  4.  Side  View  of  Evaporator 
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Figure  5.  Photograph  of  208  V,  75  A,  Variable  Transformers 
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Figure  6.  Sectional  View  of  Evaporator  Showing  Tube  Bundle 
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Figure  7.  Thermocouple  Locations  on  an  Instrumented  Boiling  Tube  and 
Tube  Section  View 
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Figure   8 .    Photograph  of  Dummy  Rack 
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Figure  9.  Close-up  View  of  Turbo-B  Tube  Surface  (25  X) 
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Table  1.  EVAPORATOR  HEATERS 


Heater  Type 


Instrumented  Tube  Heaters 


Active  Tube  Heaters 


Auxiliary    Heaters 


Simulation   Heaters 


Number 


12 


Power  Rating  per  Heater 


1000W 


1000W 


4000W 


4000W 


Table   2.      COMPUTER/DATA 
ASSIGNMENT 


ACQUISITION 


Amperage  Sensor  Description 

Channel 

Array 

Tube  1 

30 

Amp(0) 

Tube  2 

31 

Amp(1) 

Tube  3 

32 

Amp(2) 

Tube  4 

33 

Amp(3) 

Tube  5 

34 

Amp(4) 

Active  Heater  Group  1 

35 

Amp(5) 

Active  Heater  Group  2 

36 

Amp(6) 

Active  Heater  Group  3 

37 

Amp(7) 

Active  Heater  Group  4 

38 

Amp(8) 

Active  Heater  Group  5 

39 

Amp(9) 

Auxiliary    Heaters 

25 

Amp(10) 

Simulation   Heaters 

26 

Amp(11) 

Voltage  Sensor  Description 

Channel 

Array 

Instrumented/Active 

27 

Volt(0) 

Simulation   Heaters 

28 

Volt(1) 

Auxiliary    Heaters 

29 

Volt(2) 

26 


Table  2.  COMPUTER/DATA  ACQUISITION  ASSIGNMENT  (CONT. ) (cont . ) 


Thermocouple    Description 

Channel 

Array  in  code 

Vapor  1-Top  of  Condenser 

00 

T(0) 

Vapor  2-Top  of  Condenser 

01 

T(1) 

Vapor  3-Top  of  Evaporator 

02 

T(2) 

Liquid  1-Top  of  bundle 

03 

T(3) 

Liquid  2-Top  of  bundle 

04 

T(4) 

Liquid  3-Bottom  of  bundle 

05 

T(5) 

Tube  1,No.  1 

40 

T(6) 

Tube  1,No.  2 

41 

T(7) 

Tube  1  ,No.  3 

42 

T(8) 

Tube  1,No.  4 

43 

T(9) 

Tube  1,No.  5 

44 

T(10) 

Tube  1  ,No.  6 

45 

T(11) 

Tube  2,No.  1 

46 

T(12) 

Tube  2,No.  2 

47 

T(13) 

Tube  2,No.  3 

48 

T(14) 

Tube  2,No.  4 

49 

T(15) 

Tube  2,No.  5 

50 

T(16) 

Tube  2,No.  6 

51 

T(17) 

Tube  3,No.  1 

52 

T(18) 

Tube  3,No.  2 

53 

T(19) 

Tube  3,No.  3 

54 

T(20) 

Tube  3,No.  4 

55 

T(21) 

Tube  3, No.  5 

56 

T(22) 

Tube  3,No.  6 

57 

T(23) 

Tube  4, No.  1 

58 

T(24) 

Tube  4, No.  2 

59 

T(25) 

Tube  4, No.  3 

60 

T(26) 

Tube  4, No.  4 

61 

T(27) 

Tube  4,No.  5 

62 

T(28) 

Tube  4,No.  6 

63 

T(29) 

Tube  5,No.  1 

64 

T(30) 

Tube  5,No.  2 

65 

T(31) 

Tube  5, No.  3 

66 

T(32) 

Tube  5, No.  4 

67 

T(33) 

Tube  5,No.  5 

68 

T(34) 

Tube  5, No.  6 

69 

T(35) 
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IV.  EXPERIMENTAL  PROCEDURES 

A.    REMOVAL  OF  THE  TUBE  BUNDLE  AND  BUNDLE  DISASSEMBLY 

Before  starting  the  removal  of  the  tube  bundle  from  the  evaporator, 
the  front  glass  viewing  windows  were  carefully  removed.  Next,  all 
thermocouple  wires  and  tube  heater  electrical  connections  were 
disconnected.  After  this  was  completed,  the  nuts  securing  the  backing 
plate  and  support  block  were  removed,  and  the  tube  bundle  was  taken  out 
from  the  back  of  the  evaporator. 

When  the  bundle  needed  to  be  disassembled,  it  was  ensured  that  there 
was  a  clean  working  surface.  The  first  task  was  to  remove  the  plexiglas 
plate  attached  to  one  end  of  the  aluminum  baffle  plates  (ie.  the  front  of 
the  bundle  assembly)  by  four  screws.  The  ten  screws  on  the  side  of  each 
aluminum  baffle  plate  were  then  removed  (these  were  attached  to  the  dummy 
tubes  down  each  side  of  the  bundle)  .  The  aluminum  plates  were  then  pulled 
off  the  bundle.  The  four  corner  dummy  tubes  (two  top  and  two  bottom) 
remained  attached  to  the  tube  bundle  support  block  as  they  were 
countersunk  into  the  block.  The  six  outer  smooth  tubes  (three  per  side) 
could  be  easily  pulled  from  the  aluminum  plates  as  they  were  attached  only 
by  the  screws  already  removed.  The  other  ten  smooth  tubes  were  then 
unscrewed  from  the  tube  bundle  support  block  as  seen  in  Figure  10.  The 
smooth  tubes  were  engraved  to  ensure  proper  identification  during 
reassembly.  With  these  tubes  and  aluminum  baffle  plates  removed,  only  the 
instrumented  and  active  enhanced  heater  tubes  remained.   These  tubes  were 
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removed  from  the  support  block  by  loosening  the  outer  0— ring  compression 
plate,  disconnecting  the  active  heater  tube  wired  pairs,  and  pulling  the 
tubes  from  the  block.  Reassembly  of  the  tube  bundle  is  done  by  reversal 
of  this  procedure. 

B.    SYSTEM  CLEAN-UP 

If  the  system  had  been  previously  operated  with  refrigeration  oil  (or 
contaminated  from  some  other  source)  it  had  to  be  thoroughly  cleaned.  To 
accomplish  this,  the  entire  apparatus  had  to  be  taken  apart  and  cleaned  in 
the  following  manner. 

After  removal  of  the  refrigerant  (R-113  by  directly  draining  into  5 
gallon  drums  via  drain  valve  R-5  (see  Figure  2  in  Chapter  3)  at  the  bottom 
of  the  evaporator  or  R-114  by  boiling  off  into  the  storage  tank  by  opening 
R-l  and  R— 7)  and  with  the  system  at  atmospheric  pressure,  all  electrical 
connections  to  the  bundle  were  disconnected  and  the  front  viewing  glass 
windows  were  removed.  An  electric  fan  was  used  for  safety  to  ensure 
proper  ventilation.  The  tube  bundle  was  then  removed  and  disassembled  as 
described  in  section  A;  the  dummy  tube  rack  was  also  removed. 

Having  removed  the  tubes  from  the  tube  bundle,  they  were  individually 
washed  with  warm  water,  rinsed  and  then  wiped  down  with  acetone.  The 
smooth  tubes  were  cleaned  with  Copper  Brite  (a  commercial  copper  cleaning 
product)  to  remove  any  oxidation.  They  were  also  wiped  down  with  warm 
water  and  then  with  acetone.  The  same  procedure  was  followed  for  the 
Turbo— B  tubes  except  they  were  not  cleaned  with  Copper  Brite  for  fear  of 
clogging  the  channels.  During  the  cleaning  process,  a  soft  bristled 
toothbrush  was  used  to  ensure  the  enhanced  surface  was  cleaned  properly, 
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exercising  care  not  to  interfere  with  the  tube  surface.   The  evaporator 
shell  was  cleaned  in  a  similar  manner,  using  warm  water  and  acetone. 

C.  INSTALLATION  OF  THE  TUBE  BUNDLE 

Once  the  tube  bundle  had  been  cleaned  and  reassembled  (see  section  A)  , 
and  before  tightening  the  backing  plate  nuts,  the  whole  assembly  was 
carefully  guided  back  into  the  evaporator  section,  ensuring  the  plexiglas 
viewing  cover  of  the  tube  bundle  was  not  damaged.  After  the  bundle  was  in 
position,  it  was  ensured  that  the  dummy  tube  rack  was  properly  positioned 
below  the  bundle  and  that  the  vapor  thermocouple  positions  were  still  1.75 
cm  above  the  bundle.  Then,  all  the  nuts  were  tightened  equally  on  opposite 
sides  to  give  equal  compression  on  the  gasket.  To  replace  the  front 
window,  very  small,  equal  torques  (using  a  torque  wrench)  were  applied 
circumferentially  to  each  nut  on  the  outer  ring  support  in  turn.  After 
the  window  was  in  place ,  each  tube  (which  extended  through  the  outer  0— 
ring  compression  plate)  was  lightly  tapped  forward  so  as  to  touch  the 
front-viewing  window.  The  backing  plate  was  then  tightened  and  the 
individual  tube  0-rings  compressed,  providing  a  good  seal  for  the  system. 
The  compression  plate  had  grooves  for  the  tube  0— rings  to  sit  in  to  help 
with  proper  alignment  and  ensure  a  good  seal. 

D.  SYSTEM  LEAKAGE  TEST 

After  the  system  was  isolated  from  the  atmosphere  and  system  integrity 
was  restored,  a  Seargent  Welch  10  SCFM  vacuum  pump  was  connected  to  the 
apparatus  (via  valves  R-l  and  R-8)  and  the  pressure  taken  down  to  25  inHg 
vacuum.   Valves  R— 1  and  R— 8  were  then  secured  and  the  system  was  left 
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untouched  for  at  least  24  hours  to  see  if  there  was  any  air  leakage  in. 
If  there  was  significant  leakage  (>1  inHg  over  24  hours) ,  then  the  vacuum 
was  broken  by  cracking  open  valve  R— 2  slowly  (this  ensured  that  no 
moisture  entered  the  system) .  The  system  was  pressurized  (with  air)  to  15 
psig  through  valve  R-2 .  Large  leaks  could  then  be  detected  by  simply 
listening  to  the  air  issuing  from  the  system;  small  leaks  were  detected  by 
spraying  a  soapy  water  solution  to  all  surfaces  where  leaks  were  most 
likely  to  occur  (front  viewing  glass  gaskets,  backing  plate  gasket,  all 
fittings/valves  coming  off  the  condenser/evaporator,  0— ring  seals  of  the 
bundle  tubes  etc) .  Extreme  care  must  be  taken  to  ensure  no  moisture 
enters  the  inside  of  the  heated  tubes  where  the  heater  wires  protrude. 
After  all  leaks  were  detected  and  corrected,  the  system  was  again 
subjected  to  a  vacuum  for  another  24  hour  period.  If  the  vacuum  held, 
then  the  system  was  ready  to  receive  refrigerant.  If  not,  the  above  leak 
correction  test  was  repeated. 

E .    REFRIGERANT 
1.   Fill 

a.   From  System  Storage  Tank 

A  refrigerant  storage  tank  was  used  to  store  R-114  during 
modification/repairs  to  the  system.  The  storage  tank  prevented  discharge 
of  the  R-114  into  the  atmosphere  and  made  the  experimentation  less  costly. 
To  fill  the  evaporator  with  R-114  from  the  storage  tank,  the  ethylene 
glycol/water  coolant  temperature  was  first  reduced  to  —15  °C.  The  system 
pressure  was  then  maintained  below  the  storage  tank  pressure  (vapor 
pressure  of  R-114  at  20  °C  is  approximately  15  psig)  by  circulation  of  the 
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coolant  through  the  condenser  test  tubes  and  auxiliary  coils.  Valves  R-6 
and  R-4  were  then  opened  to  draw  the  R-114  from  the  storage  tank  to  the 
evaporator.  The  amount  of  refrigerant  that  was  transferred  was  controlled 
by  throttling  valve  R-6  to  obtain  the  desired  level  in  the  evaporator.  If 
required,  additional  R-114  could  be  transferred  from  a  68  kg  storage 
cylinder  to  the  system  using  valve  R-2  (see  section  l.b). 

b.       From  Refrigerant  Storage  Cylinder 

To  fill  the  apparatus  from  the  68  kg  storage  cylinder,  the 
system  pressure  was  reduced  in  the  same  way  as  above.  A  hose  assembly 
containing  a  Drierite  gas  purifier  was  connected  between  the  storage 
cylinder  and  valve  R— 2 .  A  gas  purifier  was  used  not  only  to  remove  all 
impurities,  but  also  to  remove  any  water  from  the  refrigerant.  Once  in 
place,  both  the  storage  cylinder  valve  and  R-2  were  opened  until  the 
desired  refrigerant  level  was  reached  in  the  evaporator. 
2 .   Removal  to  the  Storage  Tank 

For  tube  replacement,  system  maintenance  or  system  clean  up,  the 
R-114  was  transferred  to  the  storage  tank.  The  ethylene  glycol/water 
coolant  temperature  flowing  through  the  storage  tank  was  cooled  to  -15  °C; 
valves  R-7  and  R-8  were  opened  and  the  vacuum  pump  was  turned  on  to  put 
the  storage  tank  under  vacuum.  Once  the  storage  tank  was  under  a  20  inHg 
vacuum,  valve  R-8  was  shut  and  valve  R-l  was  opened.  Next  R— 114  was  slowly 
boiled  off  to  the  storage  tank  by  using  the  tube  bundle,  simulation  and 
auxiliary  heaters  at  a  heat  flux  of  600  kW/m2  (slow  boiling  is  important 
to  ensure  minimum  transfer  of  oil).  As  the  refrigerant  level  decreased, 
individual  heaters  were  turned  off  to  ensure  none  were  uncovered.   Once 
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below  the  level  of  the  heaters,  the  final  few  cm  of  R-114  was  boiled  off 
using  heat  from  the  atmosphere.  Once  all  of  the  R-114  was  transferred, 
valves  R— 1  and  R— 7  were  shut. 

F .    OPERATION 

1.  System  Startup,  Securing  and  Emergency  Procedures 
See  Appendix  D 

2.  Normal  Operation 

The  evaporator  was  filled  with  R-114  to  a  level  of  approximately 
10  cm  above  the  top  tubes  in  the  bundle.  Prior  to  operating  the  system, 
the  8  ton  refrigeration  unit  was  run  for  approximately  an  hour  to  reduce 
the  ethylene  glycol/water  coolant  in  the  sump  to  a  temperature  of  —15  °C . 
The  pressure  in  the  evaporator/condenser  was  usually  12  to  15  psig  if  the 
system  had  been  secured  overnight.  As  the  sump  was  brought  to 
temperature,  the  data  acquisition  system  and  computer  were  turned  on. 
This  allowed  the  temperature  in  the  system  to  be  monitored  during  cool- 
down  to  saturation  conditions.  With  this  and  pump  number  one  running, 
one  auxiliary  condenser  coil  and  the  four  condenser  test  tubes  were  used 
to  bring  the  pool  down  to  a  subcooled  condition  (for  R-114,  approximately 
1  °C  on  all  three  pool  thermocouples) .  Subcooling  of  the  refrigerant  was 
done  to  ensure  the  pool  had  an  evenly  distributed  temperature  prior  to 
starting  a  run.  After  reaching  this  subcooled  condition,  all  coolant 
supply  to  the  condenser  was  secured.  The  pool  was  then  allowed  to  'heat 
up'  by  conduction  from  the  surroundings.  Once  a  saturation  temperature  of 
2.2  °C  was  reached,  the  instrumented  tube(s)  (and  simulation  heaters  for 
test  7)  was/were  switched  on  and  set  to  the  desired  heat  flux  value. 
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This  lengthy  procedure  was  done  to  prevent  the  tubes  from  prematurely 
nucleating.  The  heat  flux  of  the  instrumented  tubes  was  then  slowly 
increased  at  desired  intervals  by  adjusting  the  rheostat.  For  increasing 
heat  flux,  the  data  was  taken  with  very  small  heat  flux  increments  (every 
1000  kW/m2)  ,  waiting  at  least  5  minutes  to  attain  steady  state  conditions 
at  each  heat  flux.  At  all  regions  of  the  boiling  curve  (and  especially 
near  the  onset  of  nucleate  boiling) ,  two  readings  were  taken  at  each  heat 
flux  to  ensure  accuracy.  The  bundle  was  continuously  visually  monitored 
through  the  observation  windows.  Figure  11  shows  the  tube  bundle 
arrangements  used  during  the  experimentation.  Test  one  was  with  only  one 
instrumented  tube  turned  on  at  any  position  within  the  bundle.  Test  two 
was  with  instrumented  tubes  one  and  two  active  in  the  bundle.  Test  three 
was  with  instrumented  tubes  one,  two  and  three  active  in  the  bundle.  Test 
four  was  with  instrumented  tubes  one,  two,  three,  and  four  active  in  the 
bundle.  Test  five  was  with  all  five  instrumented  tubes  active  in  the 
bundle.  Test  six  was  all  five  instrumented  tubes  plus  all  five  pairs  of 
active  enhanced  heater  tubes  active  in  the  bundle.  Test  seven  was  the 
same  as  test  six  with  the  addition  of  all  five  simulation  heaters  active. 
For  each  data  set,  the  five  simulation  heaters  had  the  same  heat  flux  as 
the  tubes  within  the  bundle. 

G.    OIL  ADDITION 

During  the  bundle  experiments,  successive  amounts  of  York-C  oil  were 
added  into  the  evaporator.  Since  the  weight  of  the  refrigerant  in  the 
evaporator  was  60.3  kg,  the  amount  of  oil  corresponding  to  1%  by  weight 
was  measured  as  670  ml,  2%  1340  ml  etc.   The  oil  was  syphoned  into  the 
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evaporator  via  a  funnel/hose  connection  through  valve  R-3  by  reducing  the 
pressure  in  the  evaporator  to  less  than  15  inHg  vacuum.  Ensuring  that  no 
air  entered  the  system,  valve  R-3  was  promptly  shut  when  the  desired 
amount  of  oil  had  been  added. 

H.    DATA  REDUCTION  PROCEDURES 

The  data  reduction  program  "DRP4RH"  was  used  during  the  experiments 
for  processing  the  data  collected  (  see  Appendix  E  for  listing)  .  The 
program  was  written  in  HP  Basic  3.01  and  run  on  an  HP-9000  series 
computer.  The  characteristics  and  capabilities  of  this  software  are 
similar  to  those  provided  by  Anderson  [Ref .  13] .  The  following 
modifications  were  made: 

1.  Correction  for  pool  height  by  Chilman  [Ref.  3] 

2.  Installation  of  new  thermocouple  at  the  bottom  of  the  liquid  pool 
(bundle  inlet  temperature)  by  Chilman  [Ref.  3] 

3.  The  ability  to  obtain  data  from  one  instrumented  tube  at  any  position 
in  the  bundle . 

4.  The  natural  convection  correlation  of  Churchill  and  Chu  [Ref.  31]  for 
a  single  horizontal  cylinder  in  an  'infinite'  liquid  pool  was  added  for 
comparison  with  experimental  data. 
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Figure  10.  Photograph  of  Tube  Bundle  Support  Block 
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Figure  11.  Tube  Bundle  Arrangements  used  During  Experimentation 
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V.  RESULTS  AND  DISCUSSION 


A .    INTRODUCTION 

The  results  are  presented  in  four  sections  with  subs— sections  as 
appropriate.  The  first  section  discusses  the  preliminary  experiments 
which  led  to  modifying  the  experimental  start— up  procedure  to  include 
subcooling.  The  second  section  discusses  the  natural  convection  effects, 
nucleate  pool  boiling  phenomena  and  hysteresis  effects  within  the  tube 
bundle  in  pure  R-114.  The  third  section  discusses  similar  phenomena  for 
R-114/oil  mixtures  and  their  effects  on  the  above.  The  fourth  section 
shows  comparisons  of  data  taken  during  this  thesis  with  previously 
obtained  data  at  the  Naval  Postgraduate  School. 

A  list  of  data  files  taken  during  this  investigation  may  be  found  in 
Appendix  A.  All  data  files  used  in  this  thesis  use  the  following  filename 
sequence.  Each  file  is  composed  of  five  sets  of  alpha— numeric  characters 
used  to  describe  the  experiment. 

First  set   (2  char.)   Tube  Type 

Heat  Flux 


TB  (Turbo-B) 


Second  set  (1  char.) 
Third  set   (2  char.) 


Oil  Percent 


Fourth  set  (2  char.)    Test  Type 


I  (Increasing) 
D  (Decreasing) 

00  (0%) 

01  (1%) 

02  (2%) 

03  (3%) 
06  (6%) 
10  (10%) 

01  (test  1) 

02  (test  2) 

03  (test  3) 

04  (test  4) 
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05  (test  5) 

06  (test  6) 

07  (test  7) 

Fifth  Set  (1  char.)  Additional  tests  A-Z  (If  conducted) 
To  give  an  example,  the  filename  TBI0107  means  "Turbo-B  tube  bundle, 
increasing  heat  flux,  1%  R-114/oil  mixture  and  test  number  7".  If  more 
detail  is  desired  about  a  specific  data  set,  see  Appendix  A.  All  plot 
filenames  are  similar  to  data  filenames  except  they  start  with  the  letter 
"P" .   The  test  numbers  are  shown  in  Figure  11  in  Chapter  IV. 

All  graphs  are  plotted  showing  heat  flux  (W/m2)  along  the  ordinate  (y 
axis)  as  a  function  of  wall  superheat  (K)  along  the  abscissa  (x  axis). 
The  wall  superheat  is  defined  as  the  difference  between  the  'corrected' 
average  tube  wall  temperature  (ie.  having  accounted  for  depth  of 
thermocouple  burial)  and  the  local  liquid  saturation  temperature 
(corrected  for  hydrostatic  head  within  the  bundle)  .  The  heat  flux  was 
corrected  to  account  for  the  heat  lost  through  the  unheated  tube  ends. 
The  heat  flux  was  varied  from  600  to  100,000  W/m2  for  increasing  heat  flux. 
To  ensure  greater  detail  at  the  point  of  incipience,  the  heat  flux  was 
increased  in  small  steps;  these  settings  varied  from  test  to  test.  The 
heat  flux  values  for  decreasing  experiments  were  taken  at  prescribed 
settings  for  easy  comparison  with  past  experiments  and  future  reference. 
These  heat  flux  settings  were  lxlO5,  7.5xl0\  5xl0\  3xl0\  2xl0\  1 .  5xl0\ 
lxl0\  7xl03,4xl03,  2xl03,  lxlO3  W/m2.  Approximately  30-40  data  points  were 
taken  for  each  increasing  heat  flux  run  and  20-25  points  were  taken  for 
each  decreasing  heat  flux  run. 
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B.    PRELIMINARY  EXPERIMENTS 

After  the  pure  R— 113  was  removed  following  Chilman's  [Ref.  3] 
experiments  and  the  apparatus  and  system  cleaned,  pure  R-114  was  added 
from  the  storage  cylinders.  Five  tests  were  conducted  using  this  R-114 
and  these  are  shown  in  Figure  12.  The  first  test  (TBI0001A)  was  test  one 
with  the  top  tube  activated.  The  procedure  described  by  Eraydin  [Ref.  28] 
was  followed,  but  the  plot  of  test  TBI0001A  and  Eraydin' s  data  for  test 
one  (also  shown  in  Figure  12)  produced  significantly  different  results  in 
the  natural  convection  (NC)  region.  The  data  of  Eraydin  show  a  greater 
heat  transfer  coefficient  (lower  wall  superheat)  than  test  TBI0001A  which 
show  results  closer  to  the  Churchill/Chu  [Ref.  31]  correlation  (C/C)  for 
natural  convection.  The  only  difference  in  the  apparatus  between  test 
TBI0001A  and  Eyradin's  experiments  was  the  addition  of  a  third 
thermocouple  at  the  bottom  of  the  pool.  For  test  TBI0001A,  all  three  pool 
temperature  readings  were  within  +/-  0.1  °C  prior  to  recording  data.  For 
the  data  of  Eraydin,  only  the  temperature  at  the  top  of  the  pool  could  be 
checked. 

Test  TBI0001B  was  thought  to  be  a  repeat  of  test  TBI0001A,  but  for  the 
bottom  tube  in  the  bundle  (tube  5  only).  However,  upon  observation  and 
investigation,  it  was  found  that  two  tubes  (tube  one  and  tube  five)  were 
activated  due  to  the  way  the  program  DRP4  was  set  up.  Hence  the  data 
presented  is  for  tube  one  (the  top  tube)  with  tube  five  (the  bottom  tube) 
activated  as  well.  The  program  was  then  modified  to  obtain  data  for  a 
single  tube  (test  one)  at  any  position  within  the  bundle. 

Test  TBI0001C  was  conducted  using  the  bottom  tube  (tube  five)  as  a 
single  tube.   Again  following  Eraydin' s  [Ref.  28]  procedure,  partial 
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nucleation  was  observed  immediately,  explaining  why  the  data  lie  well  to 
the  left  of  the  Churchill/Chu  correlation  (ie.  a  higher  heat  transfer 
coefficient)  .  It  can  be  seen  that  test  TBI0001C  is  in  good  agreement  with 
the  single  tube  of  Eraydin  (top  tube).  The  reason  for  this  is  probably 
due  to  the  start— up  conditions,  which  were  not  carefully  monitored  for  any 
of  these  tests. 

Test  TBI0001D  was  conducted  using  only  tube  three  and  also  followed 
the  procedure  of  Eraydin.  This  displayed  similar  behavior  as  test 
TBI0001A  (tube  one)  except  that  nucleation  was  delayed,  ie.  occurred  at  a 
higher  wall  superheat. 

All  experiments  thus  far  were  conducted  following  Eraydin' s  procedure 
for  pure  R— 114.  It  was  next  decided  to  vary  his  procedure  slightly.  For 
test  TBI0001E,  the  pool  was  first  subcooled  slightly  to  1°C.  This  ensured 
that  any  nucleation  sites  were  deactivated  and  that  the  whole  pool  was  at 
an  even  temperature.  The  pool  was  then  brought  up  to  the  required 
saturation  temperature.  However,  partial  nucleation  was  still  observed 
during  the  run.  The  plot  of  TBI0001E  is  similar  to  test  TBI0001C  and 
Eraydin' s  test  one. 

With  all  of  these  confusing  results,  it  was  decided  to  empty,  clean 
and  recharge  the  evaporator  with  fresh  R-114.  Upon  boiling  off  the  'old' 
R— 114,  a  small  quantity  of  oil  contamination  was  found  in  the  bottom  of 
the  evaporator.  This  could  have  been  either  vacuum  pump  oil  (which  may 
have  leaked  into  the  apparatus)  or  refrigerant  oil  that  entered  the  system 
with  the  R-114  charge  (which  should  be  minimal) .  A  third  possibility 
could  be  oil  from  previous  refrigerant/oil  mixture  experiments.  However, 
this  is  not  likely  since  the  system  had  been  completely  stripped  and 
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cleaned  twice  since  the  last  mixture  tests  had  been  conducted.  Samples  of 
each  type  of  oil  (vacuum  pump  oil,  miscible  refrigerant  oil  (York-C)  and 
a  sample  of  the  contamination)  were  sent  away  for  analysis  and  the  results 
were  still  forthcoming  at  the  time  of  writing  this  thesis.  From  the  color 
of  the  contamination,  it  would  appear  to  have  originated  from  the  vacuum 
pump.  This  would  seem  to  be  more  logical  since  Chilman  had  experienced 
problems  with  the  vacuum  pump. 

The  system,  including  the  bundle,  were  cleaned  thoroughly  and  fresh 
R-114  was  added.  However,  a  gas  purifier  was  utilized  to  ensure  that  only 
pure,  clean  R-114  (with  no  moisture)  was  added  (see  Chapter  IV  for  further 
details).  Test  TBI0001F  then  was  conducted  following  Eraydin's  [Ref.  28] 
procedure  with  no  subcooling  and  'ignoring'  the  bottom  pool  thermocouple 
value .  The  data  showed  good  agreement  with  Eraydin  in  the  natural 
convection  region,  but  a  lower  heat-transfer  coefficient  in  the  boiling 
region.  A  possible  explanation  for  this  lower  heat  transfer  is  the  fact 
that  each  individual  tube  was  fully  cleaned  prior  to  adding  the  new  R-114 
and  the  surface  characteristics  may  have  been  modified  in  some  way. 

Test  TBI0001G  was  conducted  using  the  same  procedure  except  the  pool 
was  initially  subcooled  down  to  1°C  for  30  minutes.  This  ensured  the  pool 
had  an  even  temperature  distribution  throughout.  The  data  were  now  much 
closer  to  the  Churchill/Chu  [Ref.  31]  correlation.  However,  some 
premature  nucleation  was  still  observed.  Test  TBI0001J  repeated  the  above 
test  with  the  pool  temperature  subcooled  to  1°C  for  at  least  one  hour  to 
further  deactivate  any  remaining  nucleation  sites  within  the  bundle  prior 
to  starting  experiments.  The  data  then  agreed  with  the  Churchill/Chu 
correlation  as  seen  in  Figure  12.   It  became  apparent  that  premature 
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nucleation  could  affect  the  natural  convection  data  significantly. 
Therefore,  for  all  subsequent  tests,  this  same  procedure  was  adopted  with 
strict  observation  of  the  bundle  to  ensure  no  premature  nucleation 
occurred. 

C.    PURE  R-114  TURBO -B  TUBE  BUNDLE  EXPERIMENTS 
1.   Test  One  for  Different  Tube  Positions 

The  first  set  of  experiments  conducted  were  performance  tests  in 
pure  R-114.  Figure  13  shows  increasing  heat  flux  for  a  single  tube  within 
the  bundle  at  different  positions  (positions  1,  3,  and  5)  while  Figure  14 
shows  corresponding  data  (including  typical  uncertainties)  for  decreasing 
heat  flux.  All  three  tube  positions  agree  closely  with  the  Churchill/Chu 
[Ref.  31]  correlation  in  the  natural  convection  region.  The  difference 
between  position  1  and  position  3  and  5  may  be  that  position  1  is  affected 
by  the  fact  that  the  flow  is  'free'  to  expand  after  leaving  the  bundle. 
This  difference  may  also  be  due  to  wall  temperature  uncertainty  due  to 
differences  in  the  fabrication  process  (see  uncertainty  analysis  Appendix 
C) ;  however  in  the  natural  convection  region,  this  uncertainty  is  low  due 
to  relatively  high  values  of  wall  superheat.  Figure  13  also  shows  that 
tube  position  within  the  bundle  may  influence  the  point  of  incipience. 
Bergles  and  Rohsenow  [Ref.  32]  have  studied  the  incipient  point  in  more 
detail.  They  concluded  that  nucleation  was  controlled  by  a  nucleation 
parameter,  N,  given  by 


(a)(Ts)(vfg) 

N  -  


(PvXhfg) 
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which  creates  an  incipient  boiling  superheat  given  by 

(TH  -  Ts)  =  2(N)/r. 
In  the  above  expression,  r  is  the  local  bubble  radius.  Calculation  of  the 
nucleation  parameter  for  pure  R— 114  showed  that  as  saturation  pressure 
increases,  the  nucleation  parameter  decreases.  Assuming  that  the  radius 
of  curvature  of  the  forming  bubbles  is  constant  (which  is  reasonable  for 
a  Turbo-B  surface  which  has  large,  regularly  spaced  cavities),  then  (TH  — 
Ts)  also  decreases  and  nucleation  may  be  expected  to  occur  earlier.  This 
certainly  seems  to  be  verified  in  Figure  13  where  the  incipient  point 
occurs  earlier  (lower  wall  superheat)  for  a  lower  tube  ie .  where  there  is 
an  increase  in  the  local  saturation  pressure.  You  et  al .  [Ref.  33]  also 
showed  a  decrease  (approx.  30%)  in  the  average  incipient  superheat  as 
pressure  was  increased  from  1  to  1.5  bar  for  pool  boiling  of  FC-72  on  a 
single  tube,  offering  some  other  experimental  verification  for  this 
conclusion. 

Chilman  [Ref.  3]  conducted  test  one  for  R-113  using  the  top  tube 
only  and  varied  the  local  saturation  pressure  by  varying  the  pool  height 
in  the  evaporator.  He  found  that  the  point  of  incipience  was  delayed  when 
increasing  the  liquid  pool  height  (ie.  the  hydrostatic  pressure  head). 
Boundary  layer  effects  due  to  different  liquid  circulation  patterns  may 
have  caused  this  delay  in  nucleation  and  more  research  is  certainly  needed 
to  fully  understand  the  influence  of  pressure  on  the  point  of  incipience. 
One  experiment  that  could  be  conducted  would  be  to  vary  the  pool  height, 
but  keep  the  local  pressure  at  each  tube  constant  by  simultaneously 
varying  the  vapor  pressure  above  the  pool. 
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Once  nucleation  occurs,  Figure  13  shows  that  the  single  tube 
experiments  merge  onto  a  single  boiling  curve.  In  this  region,  there 
appears  to  be  no  effect  of  hydrostatic  pressure  head.  Figure  14  shows  the 
corresponding  decreasing  heat  flux  data  for  a  single  tube  at  the  same 
three  positions  within  the  bundle.  It  shows  no  significant  influence  of 
tube  position  in  the  bundle  for  decreasing  heat  flux.  Note  that  at  low 
heat  flux,  the  experimental  uncertainty  in  heat  flux  and  wall  superheat  is 
the  largest  (see  Appendix  C) . 

2.    Test  Two  to  Test  Seven 

Figures  15  to  20  show  test  two  to  test  seven  for  increasing  heat 
flux  with  pure  R-114.  Also  shown  for  comparison  in  each  figure  is  the 
Churchill/Chu  [Ref .  31]  correlation  although  this  is  only  truly  valid  for 
a  single  tube  in  an  infinite  pool.  Figure  15  shows  good  agreement  with 
the  Churchill/Chu  correlation  in  the  natural  convection  (NC)  region  and 
shows  no  effect  of  the  lower  tube  on  the  upper  tube.  The  incipient  point 
occurs  approximately  at  the  same  wall  superheat  for  both  tubes;  once 
boiling  the  lower  tube  has  the  higher  heat-transfer  coefficient.  This  is 
contrary  to  the  results  obtained  by  Chilman  [Ref.  3]  and  Anderson  [Ref. 
13]  with  pure  R— 113,  where  the  higher  tube  had  the  better  heat  transfer. 
The  reason  for  this  difference  is  not  known,  but  may  be  due  to  the 
explosive  nature  of  incipience  for  R-114  compared  to  the  more  gradual 
partial  incipience  for  R— 113.  For  R-114  experiments,  the  pool  was 
subcooled  by  1  °C  for  only  about  hour  while  for  the  R-113  experiments,  the 
pool  was  left  in  a  subcooled  state  since  the  previous  experiment  (for  R- 
113  experiments,  the  pool  is  'heated  up'  to  saturation  conditions).   This 
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difference  in  subcooling  may  significantly  affect  the  nature  of  the 
observed  incipience.  Further  research  should  be  conducted  at  the 
incipient  point  to  address  some  of  the  questions.  Figure  16  for  three 
tubes  activated  shows  similar  behavior  as  Figure  15  (ie.  during  nucleate 
boiling,  the  lowest  tube  has  the  highest  heat— transfer  coefficient). 
Also,  Figure  16  shows  that  the  lowest  tube  seems  to  nucleate  last. 
Figures  17  and  18  for  four  and  five  tubes  activated  show  similar  trends  in 
both  the  NC  and  boiling  regions.  It  appears  that  the  tubes  nucleate  in 
order  down  the  bundle  (ie.  the  top  tube  nucleates  at  the  lowest  wall 
superheat  and  the  bottom  tube  nucleates  at  the  highest  wall  superheat) . 
For  tests  six  and  seven  (Figures  19  and  20  respectively)  ,  the  maximum 
controllable  heat  flux  was  less  than  tests  one  through  five  due  to  the  use 
of  smooth  tubes  in  the  condenser  limiting  the  condensate  rate  (and  hence 
pressure)  in  of  the  vapor  space  (if  enhanced  tubes  had  been  used,  this 
could  have  been  increased).  Figure  19  (test  six)  is  consistent  with  the 
above  trends.  Furthermore,  the  effect  of  activating  the  whole  bundle 
seems  to  cause  the  lowest  tube  to  nucleate  at  a  lower  heat  flux  and  wall 
superheat.  It  also  appears  that  there  may  be  some  influence  of  tube 
position  in  the  boiling  region;  however,  this  is  probably  due  to 
inaccuracy  in  the  wall  temperature  measurements  (see  Appendix  C) .  When 
the  simulation  heaters  are  also  activated  (test  seven,  Figure  20)  the 
trends  in  the  NC  region  are  similar  (ie.  no  effect  of  lower  tubes  on  upper 
tubes).   Full  nucleation  of  the  bundle  however,  seems  to  occur  earlier. 

Figure  21  shows  the  data  from  tube  one  for  all  seven  tests  with 
increasing  heat  flux.  This  figure  is  of  more  fundamental  interest  as  it 
shows  the  same  tube  under  different  bundle  conditions.   It  therefore  gives 
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a  better  direct  comparison  of  the  effect  of  heated  lower  tubes  as  any 
uncertainty  in  the  tube  wall  measurements  (the  largest  error  in  the 
experimental  data)  is  the  same  for  each  test  (ie.  any  effects  seen  in  the 
data  are  bundle  effects).  In  the  NC  region,  tube  one  alone  (test  one)  is 
somewhat  different.  This  may  be  due  to  'expansion'  of  the  flow  as  it 
leaves  the  top  of  the  bundle  (ie.  where  the  velocity  of  the  flow  has 
slowed  down)  or  due  to  convective  effects  by  the  addition  of  another  tube. 
For  test  two,  Figure  21  shows  an  effect  of  the  lower  tube  on  tube  one 
performance  in  the  NC  region.  For  test  three  and  all  subsequent  tests  no 
further  improvement  is  seen.  It  appears,  therefore  that  in  the  natural 
convection  region,  an  upper  tube  is  affected  by  a  lower  tube  directly 
below;  however,  when  additional  lower  tubes  are  heated,  there  is  no 
further  increase  in  performance  of  the  top  tube.  There  is  also  no  effect 
on  the  incipient  point  (apart  from  test  one  mentioned  above) .  In  the  high 
heat  flux  boiling  region,  there  is  also  little  enhancement  due  to  the 
lower  tubes.  This  is  to  be  expected  in  an  enhanced  tube  bundle,  where  the 
total  heat  transfer  at  high  heat  fluxes  is  primarily  due  to  nucleation 
from  the  tube  surface  itself,  rather  than  from  convection  around  the 
surface  from  the  tubes  below. 

Figures  22  to  27  show  tests  two  to  seven  for  decreasing  heat  flux 
in  pure  R-114.  When  comparing  tube  one  to  tube  two,  Figure  22  shows  a 
increase  in  heat  transfer  performance  of  tube  one  by  tube  two  in  the 
boiling  region  while  Figure  15  (increasing  heat  flux)  showed  the  opposite 
effect.  The  most  probable  reason  for  this  crossover  is  that  these  two 
experiments  (TBI0002)  and  (TBD0002)  were  conducted  on  different  days  and 
startup  procedures  were  slightly  different.   For  TBI0002  (increasing  heat 
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flux),  the  test  was  conducted  as  outlined  in  Chapter  IV  section  F  (ie. 
subcooled  to  1  °C ,  gradually  heated  up  with  data  taken  over  a  period  of 
approximately  4  hrs) .  For  TBD0002  (decreasing  heat  flux),  the  pool  was 
subcooled  to  1  °C,  and  then  tube  one  and  two  were  turned  on  to  the  highest 
heat  flux  (105  W/m2)  and  allowed  to  heat  up  for  30  minutes.  The  total  time 
boiling  for  TBI0002  at  the  highest  heat  flux  was  therefore  less  than  for 
TBD0002 .  For  all  other  experiments,  increasing  followed  by  decreasing 
runs  were  conducted  on  the  same  day  approximately  15  minutes  apart.  More 
research  should  be  conducted  to  investigate  nucleation  site 
activation/deactivation. 

At  high  heat  fluxes  for  test  three,  Figure  23  shows  no  heat 
transfer  performance  improvement  between  the  tubes.  However,  at  low  heat 
fluxes,  there  does  appear  to  be  an  improvement  on  tube  one  and  two  from 
tube  three.  Figures  24  to  27  show  similar  trends  at  high  and  low  heat 
fluxes  (ie.  the  top  tubes  are  further  enhanced  by  lower  tubes  at  low  heat 
fluxes)  .  It  should  also  be  noted  that  the  lowest  tube  in  any  specific 
decreasing  heat  flux  test  had  the  worst  performance. 

Figure  28  compares  tube  one  only  for  tests  one  to  seven  for 
decreasing  heat  flux.  As  stated  above,  there  appears  to  be  a  definite 
tube  enhancement  at  low  heat  fluxes  with  little  or  no  enhancement  at  high 
heat  fluxes.  This  is  probably  due  to  convective  effects  which  tend  to 
increase  the  heat  transfer  performance  of  the  upper  tubes  due  to  the 
presence  of  lower  tubes  (ie.  bubbles  coming  from  the  lower  tubes  impinge 
and  slide  over  the  upper  tubes  and  increase  the  heat  transfer) .  At  high 
heat  fluxes,  on  the  other  hand,  all  the  tubes  are  nucleating  so  vigorously 
that  these  'sliding'  bubbles  have  little  or  no  noticeable  affect  on  the 
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overall  performance.  This  supports  the  hypothesis  of  Cornwell  [Ref.  7] 
that  total  heat  transfer  in  a  bundle  is  due  to  a  summation  of  convective 
and  nucleation  heat  transfer  phenomena.  Similar  trends  were  found  by 
Anderson  [Ref.  13]  and  Akcasayar  [Ref.  25]  for  smooth  and  finned  tube 
bundles  respectively  (using  the  same  apparatus)  and  also  by  Arai  et  al . 
[Ref.  20]  for  a  Thermoexcel-E  tube  bundle.  However,  Akcasayar  [Ref.]  did 
not  find  such  an  enhancement  effect  for  a  High  Flux  tube  bundle  indicating 
that  at  low  heat  fluxes,  a  porous  coated  surface  already  has  a  significant 
number  of  active  nucleation  sites  such  that  impinging  bubbles  from  below 
have  little  or  no  added  effect.  Turbo— B  is  more  similar  to  a  Thermoexcel— 
E  surface  and  at  low  heat  fluxes,  these  two  types  of  surface  obviously 
exhibit  different  nucleation  characteristics  to  those  of  a  porous  coated 
tube . 


D.    R-114/0IL  MIXTURES  TURBO- B  TUBE  BUNDLE  EXPERIMENTS 
1.   Tests  with  1%  and  2%  oil 

Only  four  experiments  were  conducted  with  a  1%  and  2%  R-114/oil 
mixture.  These  were  tests  one  and  seven  for  both  increasing  and 
decreasing  heat  flux;  no  experiments  were  conducted  for  tests  two  through 
six.  Figure  29  shows  test  one  at  1%  oil  concentration  for  increasing  and 
decreasing  heat  flux,  clearly  showing  a  hysteresis  pattern.  Compared  with 
pure  R-114  (Figures  13  and  14) ,  Figure  29  shows  no  apparent  effect  of  oil 
on  the  heat  transfer  in  either  the  NC  or  boiling  regions. 

Figure  30  shows  test  seven  for  increasing  heat  flux  (1%  oil). 
Compared  with  Figure  20  (pure  R-114),  there  are  similar  trends  (ie.  no 
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effect  in  the  NC  or  boiling  regions) .  The  tubes  again  appear  to  be 
nucleating  'in  order'  down  the  bundle,  as  found  with  test  seven  in  pure  R- 
114.  For  decreasing  heat  flux,  Figure  31  shows  similar  trends  to  test 
seven  in  pure  R-114  (Figure  27) .  Thus  a  1%  oil  concentration  appears  to 
have  little  or  no  effect  on  bundle  performance  for  both  increasing  and 
decreasing  heat  flux. 

Figure  32  shows  test  one  for  a  2%  oil  concentration  for 
increasing  and  decreasing  heat  flux.  Figures  33  and  34  show  test  seven 
for  increasing  and  decreasing  heat  flux  respectfully  for  2%  oil 
concentration.  All  three  graphs  (Figures  32  to  34)  are  similar  to  those 
for  pure  and  1%  oil  concentrations,  showing  that  2%  oil  also  has  little 
effect  on  overall  bundle  heat  transfer  performance. 
2.   Tests  with  3%  oil 

Nine  experiments  were  conducted  with  a  3%  R-114/oil  mixture.  In 
addition  to  test  one  and  seven  (conducted  for  both  increasing  and 
decreasing  heat  flux  as  before)  ,  tests  two  through  six  were  conducted  for 
decreasing  heat  flux  only.  Figure  35  shows  test  one  with  3%  oil  for  both 
increasing  and  decreasing  heat  flux.  Again,  the  figure  clearly  shows 
hysteresis  effects  between  the  increasing  and  decreasing  experiments.  As 
with  previous  oil  percentages  (Figures  13,  14,  29,  and  32),  it  shows  there 
is  no  apparent  effect  of  oil  in  the  NC  region,  similar  to  previous  test 
one  data  for  other  oil  percentages.  However,  at  the  highest  heat  flux 
(100  kW/m2)  there  is  an  increase  in  the  heat  transfer  of  about  10%.  This 
is  similar  to  the  increases  found  by  Burkhardt  and  Hahne  [Ref.23]  in  a 
finned  tube  bundle  and  Arai  et  al .  [Ref.  20].   Figure  36  shows  test  seven 
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for  increasing  heat  flux.  The  Churchill/Chu  [Ref.  31]  correlation  for 
pure  refrigerant  is  plotted  for  comparison  only.  Agreement  is  good, 
demonstrating  that  in  the  NC  region,  in  addition  to  there  being  no  effect 
of  tube  position,  there  is  also  no  apparent  effect  of  oil  concentration  on 
the  heat-transfer  coefficient.  As  before,  the  tubes  appear  to  be 
nucleating  'in  order'  (ie.  top  tube  nucleates  first  with  the  bottom  tube 
nucleating  last) . 

Figures  37  to  42  show  data  from  tests  two  to  seven  for  decreasing 
heat  flux  only.  All  show  no  effect  of  lower  tubes  on  upper  tubes  (within 
the  bundle)  in  the  boiling  region  at  high  heat  fluxes  (the  small  amount  of 
scatter  is  probably  due  to  inaccuracies  in  the  wall  temperature 
measurements) .  Each  successive  figure  shows  that  the  lowest  tube  has  the 
lowest  heat-transfer  coefficient;  this  tube  is  then  enhanced  by  the 
activation  of  tubes  below  it.  Again  it  should  be  noted  that  the 
experimental  uncertainty  is  larger  at  low  heat  fluxes.  At  all  oil 
concentrations,  tube  five  is  seen  to  have  the  lowest  heat  transfer 
performance.  According  to  Chilman  [Ref.  3],  tube  five  had  the  highest 
uncertainty  in  the  wall  temperature  measurements  and  this  might  be  the 
cause  of  this  discrepancy. 

If  one  compares  Figure  42  for  3%  oil  with  Figure  28  for  pure  R- 
114,  it  can  again  be  seen  that  there  is  a  small  increase  in  the  bundle 
heat-transfer  coefficient  for  the  R-114/oil  mixture  at  the  highest  heat 
fluxes.  For  all  tests  with  oil  added,  significant  foaming  was  observed  at 
the  pool  surface  and  this  may  be  the  cause  of  this  increase  in  heat 
transfer.  Schlager  et  al .  [ref.  21]  in  their  review  article  point  out 
that  for  certain  conditions  (typically  low  pressure  and  high  heat  flux), 
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the  heat-transfer  coefficient  increases  at  low  oil  concentrations;  they 
attributed  this  to  foaming.  Figure  43  compares  test  one  to  seven  for  tube 
one  for  decreasing  heat  flux.  As  before  with  pure  R-114  (Figure  28), 
there  appears  to  be  a  definite  increase  in  performance  of  the  upper  tubes 
by  lower  tubes  at  the  low  heat  fluxes  due  to  convection  effects,  with 
little  or  no  such  increase  at  high  heat  fluxes. 
3.    Tests  with  6%  oil 

The  same  nine  tests  as  with  3%  oil  were  conducted  with  a  6%  R- 
114/oil  mixture.  Figure  44  shows  test  one  for  both  increasing  and 
decreasing  heat  flux.  It  clearly  shows  a  hysteresis  'loop'  between 
increasing  and  decreasing  experiments.  In  comparison  with  other  oil 
concentrations,  the  point  of  incipience  occurs  at  a  slightly  lower  heat 
flux.  There  also  appears  to  be  a  small  degradation  in  performance  (10-15% 
compared  with  3%  oil  concentration)  at  the  highest  heat  flux  (100  kW/m2) 
due  to  the  oil,  but  there  is  no  apparent  effect  in  the  NC  region.  Figure 
45  shows  test  seven  for  increasing  heat  flux.  As  before,  there  is  no 
apparent  effect  of  oil  on  the  heat  transfer  in  the  NC  region  and  the  tubes 
appear  to  be  nucleating  'in  order'.  The  point  of  incipience  also  seems 
unaffected  by  the  presence  of  the  oil. 

Figures  46  to  51  show  tests  two  to  seven  for  decreasing  heat 
flux.  At  the  highest  heat  fluxes,  there  is  a  similar  small  degradation 
in  the  heat  transfer  as  found  with  test  one  (10-15%)  when  compared  with  a 
3%  oil  concentration  (Figures  37  to  42) .  When  compared  with  pure  R-114 
(Figures  22  to  27),  there  is  neither  enhancement  nor  degradation, 
indicating  that  any  enhancement  provided  by  3%  oil  is  offset  by  6%  oil. 
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At  low  heat  fluxes,  the  data  are  not  only  very  similar  to  that  for  pure 
refrigerant,  but  also  to  the  other  R-114/oil  mixtures  (ie.  at  low  heat 
fluxes,  there  is  no  effect  on  heat  transfer  at  any  oil  concentration). 
Figure  52  compares  tests  one  to  seven  for  tube  one  for  decreasing  heat 
flux.  As  before,  there  appears  to  be  the  same  convective  enhancement  at 
low  heat  fluxes  with  no  enhancement  (due  to  the  successive  activation  of 
lower  tubes  within  the  bundle)  at  high  heat  fluxes. 
4.    Tests  with  10%  oil 

The  same  nine  tests  were  repeated  for  an  R-114/oil  mixture  with 
10%  oil.  Figure  53  shows  test  one  for  both  increasing  and  decreasing 
heat  flux.  Incipience  occurred  at  a  slightly  heat  higher  flux  than  both  3% 
and  6%  oil  concentrations  indicating  that  there  appears  to  be  no 
systematic  increase  or  decrease  in  this  point  with  increase  in  oil 
concentration.  More  importantly,  there  is  a  significant  decrease  in  the 
heat  transfer  at  the  highest  heat  fluxes  (20%)  when  compared  with  pure  R- 
114.  This  is  probably  due  to  the  re— entrant  channels  becoming  'clogged' 
with  oil  as  the  R-114/oil  mixture  is  'transported'  to  the  surface  at  a 
high  rate.  Figure  54  shows  test  seven  for  increasing  heat  flux.  As 
before,  this  shows  that  the  NC  region  is  unaffected  by  either  oil 
concentration  or  lower  tubes  in  the  bundle.  At  the  highest  heat  fluxes 
available  (40  kW/m2)  there  appears  to  be  little  decrease  in  the  bundle 
performance  (when  compared  to  pure  R-114)  due  to  the  oil.  This  indicates 
that  at  typical  evaporator  operating  heat  fluxes,  the  presence  of  oil  does 
not  significant  effect  the  heat  transfer  enhancement  process.   At  higher 
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heat  fluxes,  however,  the  effect  of  oil  appears  to  be  very  significant  as 
seen  in  Figure  53. 

Figures  55  to  60  show  test  two  to  test  six  for  decreasing  heat 
flux  with  10%  oil.  At  low  heat  fluxes,  there  seems  to  be  no  effect  of  the 
oil  on  the  local  heat  transfer  performance.  However,  at  high  heat  fluxes, 
there  is  a  significant  decrease  in  performance.  Interestingly,  if  one 
compares  Figures  53,  and  Figures  55  to  58  at  high  heat  fluxes,  the  lowest 
activated  tube  in  the  bundle  is  significantly  degraded.  The  effect  of 
activating  a  lower  tube  significantly  enhances  the  heat  transfer  from  the 
tube  directly  above  and  (to  a  lesser  degree)  the  tubes  even  higher  in  the 
bundle.  This  may  be  due  to  the  vigorous  boiling  action  of  lower  tubes 
partly  'scouring'  the  oil  rich  layer  which  'blankets'  the  upper  tubes. 
This  effect  was  also  noticeable  with  the  High  Flux  bundle  (Akcasayar  [Ref . 
25]). 

Figure  61  compares  tube  one  for  tests  one  to  seven  for  decreasing 
heat  flux.  If  one  compares  Figures  26  (0%),  43  (3%),  52  (6%)  and  61  (10%) 
for  tube  one  for  all  seven  tests,  it  is  clear  that  at  low  heat  fluxes,  the 
heat  transfer  coefficient  is  similar,  regardless  of  oil  concentration. 
Furthermore,  convective  effects  are  consistent  and  provide  similar 
enhancements  in  heat  transfer  performance  for  all  concentrations.  At  high 
heat  fluxes  at  0,  3,  and  6%  oil  concentrations,  there  is  little 
enhancement  due  to  activation  of  lower  tubes.  However,  at  10%,  there  does 
appear  to  be  a  small  heat  transfer  enhancement  due  to  activation  of  lower 
tubes.  This  was  attributed  above  to  increased  'scouring'  of  the  oil  from 
the  vicinity  of  the  Turbo-B  surface  by  the  increase  in  bubble  activity  as 
more  tubes  are  activated  within  the  bundle.   However,  for  a  practical 
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operating  heat  flux  range  between  15  and  30  kW/m2,  there  is  no  significant 
degradation  in  heat  transfer  performance  for  an  oil  concentration  of  up  to 
10%. 

E.    COMPARISON  OF  R-114/0IL  MIXTURE  EXPERIMENTS 

Figures  62  and  63  compare  tube  one  from  test  one  for  increasing  and 
decreasing  heat  flux  for  all  oil  concentrations.  Figure  62  shows  no 
effect  of  oil  in  the  NC  region,  but  some  degradation  in  the  boiling 
region.  The  correlation  of  Churchill/Chu  is  included  for  comparison.  The 
incipient  point  appears  relatively  random,  indicating  no  early  or  delayed 
nucleation  caused  by  the  presence  of  oil.  Figure  63  shows  similar 
degradation  with  a  significant  effect  of  the  oil  (20%  decrease  in  the  heat 
transfer  from  0%  to  10%  oil)  at  the  highest  heat  fluxes. 

Figures  64  and  65  compare  tube  one  from  test  seven  for  increasing 
and  decreasing  heat  flux  for  all  oil  concentrations.  As  compared  to  tube 
one  test  one  (Figures  62  and  63) ,  the  presence  of  oil  has  some  degradation 
effect  on  the  heat  transfer  performance  (15%)  in  the  NC  region;  this  may 
be  due  to  a  change  in  the  mixture  properties  which  would  tend  to  increase 
the  wall  superheat  slightly  as  shown.  As  with  Figure  62,  Figure  64  shows 
that  the  incipient  point  appears  relatively  random.  Figure  65  shows 
similar  trends  to  Figure  63  (ie.  no  effect  of  oil  at  high  heat  fluxes), 
but  shows  convective  effects  at  low  heat  fluxes.  This  was  expected  and 
previously  reported  (Figures  28,  43,  52,  and  61). 

Figures  66  and  67  show  the  average  bundle  heat-transfer  coefficient 
(ie.  an  average  of  all  five  instrumented  tubes)  as  a  function  of  heat  flux 
for  test  seven  for  increasing  and  decreasing  heat  flux  respectively  at  all 
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oil  concentration.  The  data  are  from  the  same  data  set  as  that  shown  in 
Figures  64  and  65.  Comparing  Figures  64  and  66,  degradation  is  seen  in 
the  NC  region  (15%)  due  to  the  change  in  mixture  properties  as  mentioned 
above.  However,  at  a  practical  operating  heat  flux  range  between  15  and 
30  kW/m2,  the  presence  of  up  to  10%  oil  causes  no  degradation  in  bundle 
performance  as  seen.  These  trends  over  this  heat  flux  range  were  similar 
to  that  found  for  a  High  Flux  tube  bundle  (Akcasayar  [Ref.  25]). 
Comparing  Figures  65  and  67,  similar  trends  (ie.  no  significant  effect  of 
oil)  are  found.  However,  due  to  the  limit  in  maximum  controllable  heat 
flux  for  test  seven,  data  at  'higher'  heat  fluxes  (  up  to  105  W/m2)  could 
not  be  obtained  (as  mentioned  earlier).  As  shown  in  Figures  62  and  63, 
there  may  be  a  significant  degradation  in  the  heat  transfer  performance  at 
these  higher  heat  fluxes. 

F.    COMPARISON  WITH  PREVIOUS  NPS  DATA 

Figure  68  shows  a  comparison  between  the  present  data  for  a  Turbo— B 
bundle,  the  data  of  Anderson  [Ref.  13]  for  a  smooth  tube  bundle  and  the 
data  of  Akcasayar  [Ref.  25]  for  both  a  19  fpi  and  High  Flux  tube  bundle  in 
R-114.  For  clarity,  only  test  one  (tube  one)  for  a  decreasing  heat  flux 
in  pure  R-114  has  been  shown.  Figure  68  shows  that  the  Turbo-B  tube  has 
a  significantly  lower  heat-transfer  coefficient  than  the  High  Flux  tube  at 
all  heat  flux.  This  is  surprising  since  Sugiyama  [Ref.  34]  showed  that  in 
the  single  tube  apparatus,  the  Turbo-B  tube  was  the  best  performer.  The 
reason  for  this  difference  in  behavior  is  not  known.  Also,  all  of  the  data 
appear  to  be  parallel  to  each  other;  one  may  expect  the  enhanced  tubes  to 
have  a  different  slope  to  a  smooth  tube  due  to  the  greater  amount  of 
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nucleation.  At  the  highest  heat  fluxes,  Turbo-B  and  the  19  fpi  tube 
appear  to  have  a  similar  heat  transfer  performance.  However,  the  finned 
tube  heat  flux  is  based  on  the  root  diameter;  if  the  actual  finned  area 
had  been  used,  then  the  heat  flux  would  be  significantly  lower.  Figure  68 
shows  that  the  heat  transfer  enhancement  given  by  the  Turbo— B  tube  when 
compared  to  a  smooth  tube  is  about  three  at  high  heat  fluxes  and  increases 
to  about  five  at  low  heat  fluxes. 

Figures  69  and  70  compare  the  average  bundle  heat-transfer  coefficient 
(ie.  test  seven)  for  a  given  oil  percentage  to  that  with  no  oil  for  all 
four  tube  bundles  at  heat  fluxes  of  15  and  30  kW/m2  respectively.  15  and 
30  kW/m2  were  chosen  as  being  representative  of  the  lower  and  upper  limits 
of  heat  fluxes  used  in  practical  Naval  evaporators.  At  15  kW/m2,  Figure 
69  shows  large  enhancements  for  the  smooth  and  finned  tube  bundles  for  all 
oil  concentrations,  especially  at  the  lower  oil  concentrations.  However, 
the  Turbo-B  and  High  Flux  tube  bundles  show  a  degradation  in  the  heat- 
transfer  coefficient  at  all  oil  concentrations  (approximately  5-10%  at  low 
oil  concentrations  dropping  to  nearly  25%  at  10%  oil  for  the  High  Flux 
bundle).  At  the  higher  heat  flux  (30  kW/m2) ,  Figure  70  shows  similar 
trends  as  Figure  69.  However,  the  High  Flux  bundle  now  exhibits  a  40% 
decrease  in  the  average  bundle  heat— transfer  coefficient  at  10%  oil. 

Figures  71  and  72  compare  the  average  bundle  heat-transfer  coefficient 
(ie.  test  seven)  for  each  enhanced  tube  to  that  for  the  smooth  tube  bundle 
(tested  by  Anderson  [Ref.  13])  for  all  oil  concentrations  at  heat  fluxes 
of  15  and  30  kW/m2  respectively.  At  15  kW/m2  with  pure  R-114,  Figure  71 
shows  an  enhancement  factor  of  3.7  for  the  Turbo-B  tube  bundle.  This 
enhancement  decreases  slowly  with  increasing  oil  percentage  to  a  factor  of 
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about  2.5  at  10%  oil.  This  agree  very  closely  with  the  19  fpi  bundle. 
The  High  Flux  bundle  exhibits  much  larger  enhancements,  from  over  6  at  0% 
oil  to  just  over  3  at  10%  oil.  At  30  kW/m2  with  pure  R-114,  Figure  72 
shows  an  enhancement  of  3.8  for  the  Turbo-B  tube  bundle,  decreasing  to 
about  3  at  10%  oil.  This  again  agrees  closely  with  the  19  fpi  bundle.  It 
should  be  noted  that  the  High  Flux  bundle  enhancement  has  decreased  to  a 
similar  value  and  gets  worse  than  the  other  bundles  as  the  heat  flux  is 
further  increased. 
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Figure  12.  Performance  of  Test  One  For  Preliminary  Experiments 


59 


1  \_ri_»l_     1 

o     TUBE  3 
+     TUBE  5 
C/C[REF.31] 

g 

V  " 

< 
V    « 

4 

P 

c 

0 

• 
0 

• 

1j 

1     ^ 

C 

k       ' 

. 

\       ► 

1 

> 

\5 

i 

\ 

• 
via 

I 

f 

i 

1 

.  — _  _ 

V< 

„J, 

: 

o 
o 
o 
o 
o 


o 
o 
o 
o 


o 
o 
o 


(21/Y/M)  Xmd  1V3H 


< 


< 


o 
o 


Figure  13.   Performance  of  Test  One  at  Various  Tube  Positions  for 
Increasing  Heat  Flux 
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Figure  14.  Performance  of  Test 
Decreasing  Heat  Flux 
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Figure   15.    Performance  of  Tubes  1  and  2  for  Increasing  Heat  Flux  in  Pure 
R-114 
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Figure  16.  Performance  of  Tubes  1,2,  and  3  for  Increasing  Heat  Flux  in 
Pure  R-114 
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Figure  17.  Performance  of  Tubes  1,  2,  3,  and  4  for  Increasing  Heat  Flux 
in  Pure  R-114 
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Figure  18.  Performance  of  All  Five  Tubes  for  Increasing  Heat  Flux  in 
Pure  R-114 
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Figure  19.   Performance  of  All  Five  Tubes  with  Active  Pairs  for 
Increasing  Heat  Flux  in  Pure  R-114 


66 


4a44-UU4-U 


OT 


o"o*o#o* 


IE 

.t/> 


x> 


CO 
Li.' 
t-  C\J  CO  ^-  10  UJ 
Uj  UJ  LiJ  qj  UJ  EL 
CD  CD  CD  cq  CO  0 
D     3     P     3     3    V=::- 

p  i-  h  p:  H  O 


o 
o 
o 
o 
o 


o 
o 
o 
o 


o 
o 
o 


(ZIN/M)  Xmd  1V3H 


l- 
< 

CO 


< 


o 
o 


Figure  20.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  Pure  R-114 
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Figure   21.    Comparison   of  Tests   One    to   Seven   for  Tube    1   for   Increasing 
Heat   Flux   in  Pure  R-114 


68 


A    >M     -W    JJAV    UAW    WVMWVUV  HIWWAMVWAI 


La 


°o°o§ogo°o 
o°o0o08°o0 


o 
o 
o 
o 
o 


o 
o 
o 
o 


o 
o 
o 


< 

CO 

■ 


< 


o 
o 


(3IAI/M)    Xmd    1V3H 

Figure   22.    Performance  of  Tubes  1  and  2  for  Decreasing  Heat  Flux  in  Pure 
R-114 
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Figure  23.  Performance  of  Tubes  1,  2,  and  3  for  Decreasing  Heat  Flux  in 
Pure  R-114 
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Figure  24.  Performance  of  Tubes  1 ,  2 ,  3,  and  4  for  Decreasing  Heat  Flux 
in  Pure  R-114 
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Figure  25.  Performance  of  All  Five  Tubes  for  Decreasing  Heat  Flux  in 
Pure  R-114 
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Figure  26.   Performance  of  All  Five  Tubes  with  Active  Pairs  for 
Decreasing  Heat  Flux  in  Pure  R-114 
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Figure  27.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  Pure  R-114 
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Figure   28.    Comparison   of  Tests   One    to   Seven   for   Tube    1    for   Decreasing 
Heat   Flux   in  Pure  R-114 
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Figure  30.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  R-114  with  1%  Oil 


77 


< 


o 
o 
o 
o 
o 


o 
o 
o 
o 


o 
o 
o 


o 
o 


(ZIAI/M)  Xmd  ±V3H 


Figure  31.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  R-114  with  1%  Oil 
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Figure  32   Performance  of  Tube  1  for  Increasing/Decreasing  Heat  Flux  in 
R-114  with  2%  Oil 
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Figure   33.      Performance     of     the     Bundle     with     Simulation     Heaters     for 
Increasing  Heat  Flux   in  R-114  with  2%   Oil 
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Figure  34.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  R-114  with  2%  Oil 
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Figure   35.    Performance  of  Tube  1   for  Increasing/Decreasing  Heat  Flux  in 
R-114  with   3%   Oil 
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Figure  36.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  R-114  with  3%  Oil  neacers  tor 
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Figure  37.  Performance  of  Tubes  1  and  2  for  Decreasing  Heat  Flux  in  R- 
114  with  3%  Oil 
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Figure    38.    Performance  of  Tubes   1,    2,    and  3   for  Decreasing  Heat  Flux   in 
R-114  with   3%   Oil 
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Figure  39.  Performance  of  Tubes  1,  2,  3,  and  4  for  Decreasing  Heat  Flux 
in  R-114  with  3%  Oil 


86 


< 
(f) 


< 


(ZKI/M)  Xmd  1V3H 


Figure  40.  Performance  of  All  Five  Tubes  for  Decreasing  Heat  Flux  in  R- 
114  with  3%  Oil 
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Figure  41.   Performance  of  All  Five  Tubes  with  Active  Pairs  for 
Decreasing  Heat  Flux  in  R-114  with  3%  Oil 
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Figure  42.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  R-114  with  3%  Oil 
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Figure  43.  Comparison  of  Tests  One  to  Seven  for  Tube  1  for  Decreasing 
Heat  Flux  in  R-114  with  3%  Oil 
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Figure  44.  Performance  of  Tube  1  for  Increasing/Decreasing  Heat  Flux  in 
R-114  with  6%  Oil 
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Figure  45.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  R-114  with  6%  Oil 


92 


-■ 


IT 


•ojogogogg 


, — 





« 


I 


6 

• 


„ — 





1' 


.  —  — 





, ^^^  * 




,  ««  »«s—  f 


■•»  ***  w^  wv<  wm  -. 





VW<  AW  WW  WkWI  ><wv 


i-  CM 

LU  LU 

CD  CD 

Z)  D 


•     o 


u. 


o 
o 
o 
o 
o 


o 

O 
O 
O 


O 
O 
O 


o 
o 


< 

C/) 


< 


o 
o 


(3IAI/M)    Xmd    1V3H 


Figure   46.    Performance   of  Tubes   1  and  2   for  Decreasing  Heat  Flux  in     R- 
114  with   6%   Oil 
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Figure  49.  Performance  of  All  Five  Tubes  for  Decreasing  Heat  Flux  in  R- 
114  with  6%  Oil 
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Figure  50.   Performance  of  All  Five  Tubes  with  Active  Pairs  for 
Decreasing  Heat  Flux  in  R-114  with  6%  Oil 
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Figure  51.   Performance  of  the  Bundle  with  Simulation  Heaters 
Decreasing  Heat  Flux  in  R-114  with  6%  Oil 
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Figure  52.  Comparison  of  Tests  One  to  Seven  for  Tube  1  for  Decreasing 
Heat  Flux  in  R-114  with  6%  Oil 
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Figure   53.    Performance  of  Tube  1   for  Increasing/Decreasing  Heat  Flux  in 
R-114  with   10%  Oil 
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Figure  54.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  R-114  with  10%  Oil 


101 


o 
o 


o 
o 
o 
o 
o 


o 
o 
o 
o 


o 
o 
o 


< 

CO 


< 


o 
o 


(3IAI/M)  Xnid  1V3H 

Figure  55.  Performance  of  Tubes  1  and  2  for  Decreasing  Heat  Flux  in  R- 
114  with  10%  Oil 
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Figure  56.  Performance  of  Tubes  1,  2,  and  3  for  Decreasing  Heat  Flux  in 
R-114  with  10%  Oil 
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Figure  57.  Performance  of  Tubes  1,  2,  3,  and  4  for  Decreasing  Heat  Flux 
in  R-114  with  10%  Oil 
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Figure    58.    Performance  of  All  Five  Tubes   for  Decreasing  Heat  Flux  in  R- 
114  with   10%  Oil 
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Figure  59.   Performance   of  All  Five  Tubes  with  Active   Pairs   for 
Decreasing  Heat  Flux  in  R-114  with  10%  Oil 
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Figure  60.   Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  R-114  with  10%  Oil  "eaters  tor 
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Figure  61.  Comparison  of  Tests  One  to  Seven  for  Tube  1  for  Decreasing 
Heat  Flux  in  R-114  with  10%  Oil 
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Figure  62.  Comparison  of  Test  One  for  Increasing  Heat  Flux  in  R-114  /Oil 
Mixtures 
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Figure  63.  Comparison  of  Test  One  for  Decreasing  Heat  Flux  in  R-114  /Oil 
Mixtures 
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Figure  64.  Comparison  of  Tests  One  to  Seven  Tube  One  for  Increasing  Heat 
Flux  in  R-114/0il  Mixtures 
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Figure  65.  Comparison  of  Tests  One  to  Seven  Tube  One  for  Decreasing  Heat 
Flux  in  R-114/0il  Mixtures 
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Figure  66.  Mean  Bundle  Heat-Transfer  Coefficient  for  Increasing  Heat 
Flux  in  R-114/0il  Mixtures 
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Figure  68.  Test  One  Comparison  of  Turbo-B,  Smooth,  Finned,  and  High  Flux 
Tube  Bundles  for  Decreasing  Heat  Flux  in  Pure  R-11A 
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Figure  69.  Ratio  of  Mean  Bundle  Heat-Transfer  Coefficient  With  Oil  to 
Heat-Transfer  Coefficient  Without  Oil  for  Different  Oil  Percentages  at 
a  Heat  Flux  of  15  kW/m2 
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Figure  70.  Ratio  of  Mean  Bundle  Heat-Transfer  Coefficient  With  Oil  to 
Heat-Transfer  Coefficient  Without  Oil  for  Different  Oil  Percentages  at 
a  Heat  Flux  of  30  kW/m2 
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Figure  71.  Ratio  of  Mean  Bundle  Heat-Transfer  Coefficient  of  Enhanced 
Tube  to  Heat-Transfer  Coefficient  of  Smooth  Tube  for  Different  Oil 
Percentages  at  a  Heat  Flux  of  15  kW/m2 
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Figure  72.  Ratio  of  Mean  Bundle  Heat-Transfer  Coefficient  of  Enhanced 
Tube  to  Heat-Transfer  Coefficient  of  Smooth  Tube  for  Different  Oil 
Percentages  at  a  Heat  Flux  of  30  kW/m2 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.    CONCLUSIONS 

Nucleate  boiling  data  of  R-114  at  atmospheric  pressure  were  obtained 
using  a  small  bundle  of  Turbo-B  copper  tubes.  The  data  were  obtained  for 
both  increasing  and  decreasing  heat  flux  and  at  different  oil 
concentrations.  Based  upon  the  results  pertaining  to  this  particular 
bundle  and  apparatus,  the  following  conclusions  may  be  made: 

1.  Natural  Convection  Region 

a.  For  a  single  upper  tube,  a  second  lower  tube  directly  below 
when  turned  on  does  increase  the  heat  transfer  performance  of  the  upper 
tube,  however  when  additional  lower  tubes  are  heated  no  net  increase  in 
performance  occurs . 

b.  The  presence  of  heated  lower  tubes  in  the  bundle  reduces  the 
incipient  boiling  point  of  the  upper  tubes  and  the  tubes  tend  to  nucleate 
'in  order'  (ie.  top  tube  first,  bottom  tube  last). 

c.  The  effect  of  adding  oil  to  the  refrigerant  (up  to  10%) 
reduces  the  heat-transfer  coefficient  slightly  (approximately  10-15%)  due 
to  changes  in  the  fluid  properties. 

2.  Boiling  Region 

a.  For  pure  R-114,  the  presence  of  heated  lower  tubes  on  the  top 
tube  causes  no  enhancement  at  high  heat  fluxes  (>  20  kW/m2)  ,  but  at  low 
heat  fluxes  (<  20kW/m2)  ,  there  is  a  significant  enhancement  due  to 
convective  effects. 
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b.  At  very  low  heat  fluxes  (<2  kW/m2)  ,  the  presence  of  oil  has 
little  effect  on  the  heat  transfer  performance  of  the  top  tube  in  the 
bundle.  At  higher  heat  fluxes  (>  2kW/m2)  ,  the  performance  is  enhanced  by 
10-15%  at  low  concentrations,  but  is  degraded  up  to  20%  at  10%  oil 
concentration  at  the  highest  heat  fluxes. 

c.  At  typical  operating  heat  fluxes  (15-30  kW/m2)  ,  the  bundle 
performance  is  reduced  between  5-15%  with  oil. 

B.    RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  Conduct  experiments  with  varying  pool  height,  but  keep  the  local 
pressure  at  each  tube  constant  by  simultaneously  varying  the  vapor 
pressure  above  the  pool. 

2.  Additional  experiments  with  R-113  and  R-114  should  be  conducted  to 
investigate  explosive  (R-114)  and  partial  (R-113)  incipience  at  the  onset 
of  nucleation  varying  the  time  at  the  incipience. 

3.  Some  instrumentation  should  be  added  such  that  the  flowrates 
through  the  bundle  can  be  determined.  From  these  measurements,  vapor 
quality  can  be  determined. 

4.  Metal  guide  plates  should  be  manufactured  and  placed  on  each  side 
between  the  simulation  tube  bundle  and  the  tube  bundle  itself.  This 
further  channels  the  flow  of  refrigerant  thru  the  bundle  at  high  heat 
fluxes . 

5.  Attention  needs  to  be  given  to  the  question  of  refrigerant 
disposal.  There  are  reclamation  projects  undertaken  by  most 
manufacturers;  however,  a  method  still  needs  to  be  found  to  remove  the 
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refrigerant  from  the  apparatus  into  a  container  suitable  for  such 
reclamation. 

6.  A  high  speed  camera  should  be  used  to  study  the  nucleation  process 
and  circulation  patterns  in  more  detail  in  the  bundle.  Neutrally  buoyant 
particles  might  be  placed  in  the  pool  to  facilitate  study  of  circulation 
patterns  within  the  bundle. 
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APPENDIX  A:  LIST  OF  DATA  FILE 


Table  3.  DATA  FILE  NAMES  FOR  TURBO-B 
TUBE  BUNDLE  EXPERIMENTS 


FILE  NAME 

NUMBER  OF 
DATA  POINTS 

NUMBER  OF 
INSTRUMENTED  TUBE 

PERCENT 
OF  OIL 

NUMBER  OF 

ACTIVE 

PAIRS 

NUMBER  OF 

SIMULATION 

HEATERS 

TBI0001H 

38 

1  (15) 

0 

0 

0 

TBD0001H 

38 

1  (15) 

0 

0 

0 

TBI0001I 

25 

1  (13) 

0 

0 

0 

TBD0001I 

22 

1  (13) 

0 

0 

0 

TBI0001J 

34 

1 

0 

0 

0 

TBDOOOU 

25 

1 

0 

0 

0 

TBI0002A 

38 

2 

0 

0 

0 

TBD0002 

25 

2 

0 

0 

0 

TBI0003A 

33 

3 

0 

0 

0 

TBD0003 

27 

3 

0 

0 

0 

TBI  0004 

36 

4 

0 

0 

0 

TBD0004 

24 

4 

0 

0 

0 

TBI0005 

34 

5 

0 

0 

0 

TBD0005 

24 

5 

0 

0 

0 

TBI  0006 

32 

5 

0 

5 

0 

TBD0006 

26 

5 

0 

5 

0 
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FILE  NAME 

NUMBER  OF 
DATA  POINTS 

NUMBER  OF 
INSTRUMENTED  TUBE 

PERCENT 
OF  OIL 

NUMBER  OF 

ACTIVE 

PAIRS 

NUMBER  OF 

SIMULATION 

HEATERS 

TBI0007 

31 

5 

0 

5 

5 

TBD0007 

24 

5 

0 

5 

5 

TB10101 

33 

1 

1 

0 

0 

TBD0101 

23 

1 

1 

0 

0 

TBI0107 

29 

5 

1 

5 

5 

TBD0107 

26 

5 

1 

5 

5 

TBI0201 

38 

1 

2 

0 

0 

TBD0201 

22 

1 

2 

0 

0 

TBI  0207 

25 

5 

2 

5 

5 

TBD0207 

20 

5 

2 

5 

5 

TBI0301 

34 

1 

3 

0 

0 

TBD0301 

25 

1 

3 

0 

0 

TBD0302 

24 

2 

3 

0 

0 

TBD0303 

22 

3 

3 

0 

0 

TBD0304 

24 

4 

3 

0 

0 

TBD0305 

24 

5 

3 

0 

0 

TBD0306 

21 

5 

3 

5 

0 

TBI  0307 

26 

5 

3 

5 

5 

TBD0307 

21 

5 

3 

5 

5 

TBI  0601 

32 

1 

6 

0 

0 
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FILE  NAME 

NUMBER  OF 
DATA  POINTS 

NUMBER  OF 
INSTRUMENTED  TUBE 

PERCENT 
OF  OIL 

NUMBER  OF 

ACTIVE 

PAIRS 

NUMBER  OF 

SIMULATION 

HEATERS 

TBD0601 

22 

1 

6 

0 

0 

TBD0602 

22 

2 

6 

0 

0 

TBD0603 

24 

3 

6 

0 

0 

TBD0604 

22 

4 

6 

0 

0 

TBD0605 

24 

5 

6 

0 

0 

TBD0606 

21 

5 

6 

5 

0 

TBI0607 

27 

5 

6 

5 

5 

TBD0607 

21 

5 

6 

5 

5 

TBI1001 

33 

1 

10 

0 

0 

TBD1001 

22 

1 

10 

0 

0 

TBD1002 

27 

2 

10 

0 

0 

TBD1003 

24 

3 

10 

0 

0 

TBD1004 

24 

4 

10 

0 

0 

TBD1005 

23 

5 

10 

0 

0 

TBD1006 

16 

5 

10 

5 

0 

TB1 1007 

24 

5 

10 

5 

5 

TBD1007 

17 

5 

10 

5 

5 

TBI0001A 

24 

1 

0 

0 

0 

TBI0001C 

27 

1(15) 

0 

0 

0 

TBI0001D 

33 

KI3) 

0 

0 

0 
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FILE  NAME 

NUMBER  OF 
DATA  POINTS 

NUMBER  OF 
INSTRUMENTED  TUBE 

PERCENT 
OF  OIL 

NUMBER  OF 

ACTIVE 

PAIRS 

NUMBER  OF 

SIMULATION 

HEATERS 

TBI0001E 

33 

1 

0 

0 

0 

TBI0001F 

38 

1 

0 

0 

0 

TBI0001G 

38 

1 

0 

0 

0 
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APPENDIX  B:  SAMPLE  CALCULATIONS 


Data  set  number  1  Tube  1  of  experiment  TBD1005  (Turbo-B  tube, 
decreasing  heat  flux,  10%  oil  concentration,  test  5)  was  used  for  the 
sample  calculations  in  order  to  validate  the  program  used  for  data 
acquisition  DRP4RH .   The  working  fluid  was  R-114. 

1.  Test  tube  dimensions 

Dtc  =   11.60  mm 
D0  -   14.15   mm 
Dj   -   12.70  mm 
L     =  203.2  mm 
Ly  =   25.4  mm 

2.  Measured  Parameters 

Tl  -  10.62  °C 
T2  =  10.94  °C 
T3  -  9.96  °C 
T4  ■=  10.90  °C 
T5  -=  10.04  °C 
T6  •=  9.07  °C 
Tldl  -  2.27  °C 
Tld2  -  2.21  °C 
Aas   -   3.513   V 
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Vas  -  3.189  V 


3.   Calculations 

The  heaters  power  is  first  calculated  for 
q  =  Vas(V)xAas(V)x60(V/V)xl(A/V) 
Note:  The  multiplication  factors  of  volts  and  amp  sensors  are  60  and 
1,  respectfully. 
Therefore : 

q  =   (3.189) (3.513) (60V/V)(1A/V) 
q  =    672.19  Watts 
The  tube  inside  wall  temperature  is  obtained  from  the  average  of  all 
six  thermocouple  readings. 

-  1/6(10.62  +  10.94  +  9.96  +  10.90  +  10.04  +  9.07) 

-  10.25  °C 

The  tube  outside  temperature  is  calculated  by  knowing  the  inside  wall 
temperature  using  Fourier's  Conduction  Law.  Uniform  radial  conduction  is 
assumed. 

g[ln<-^2_)] 

r_  -  r„,  - 


vo  vi 


2ii{kcu)  (L) 


where  the  second  term  on  the  right  hand  side  is  the  Fourier  conduction 
term.   If  we  define  this  term  as 
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g[ln(-^2-)] 

<f> ZI£ 

*  2ti  (*cu)  (L) 


and 


6*  -   Tuo  -  Tsatc 


where    kcu    is     the     thermal    conductivity    of    copper    and    is    calculated    as 
follows 

kcu  -  434.0   -    [0.112(7^)] 
kcu  -  434.0   -    [0.112(283.25)] 

kcu  -  402.28   W/mK 
now 


672.19  [In  (lf4|)] 
v0  "      wi      2n  (402.28)  (.2032) 


r„0  -    (10.25   -    .2501)  °  C 


T„0  -  9.98*  C 


The  liquid  saturation  temperature  at  the  top  of  the  tube  bundle  is 
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Tsat  _     tldl   +   Tld2 


Tsat  -    2'27    +  2'21 


Tsat  -  2.24 °C 

In  order  to  calculate  the  local  saturation  temperature  for  each  tube, 
correction  factors  are  needed  to  account  for  hydrostatic  pressure 
differences  between  the  tube  locations  and  the  liquid  free  surface.  This 
difference  is  calculated  by: 

AP  -  pig)  (he) 

For  the  top  tube  in  the  bundle  which  is  0.124  m  below  the  thermocouple 
measuring  pool  temperature. 

LP  -   1523.  12  (9. 81)  (0.124) 

AP  -  1852.78  Pa 

For  1852.78  Pa  pressure  difference,  corrected  saturation  temperature  is 
obtained  by  adding  0.04  °C  (from  standard  tables  for  R-114)  to  Tsat. 
Corrected  Tsat  is: 

Tsatc  -    (2.24  +  0.04)  °  C 
Tsatc  -   2.28°  C 
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Therefore,  the  wall  superheat  can  be  obtained  by  the  following: 


9h  -  Two  -  Tsatc 


9h  -  (9  .98  -  2.28)  °  C 


bb  -   7 .70°  C 


Now  that  the  wall  superheat  is  known,  we  need  to  calculate  the  heat 
flux  and  the  heat-transfer  coefficient.  To  do  this,  we  know  that  the  tube 
is  12  inches  long  and  is  heated  in  a  eight  inch  center  portion  of  the 
tube .  The  unheated  lengths  of  the  tube  are  a  one  inch  and  a  three  inch 
section  on  opposite  ends  of  the  tube.  These  unheated  lengths  have  a  fin 
effect  during  the  heat  transfer  process  to  the  evaporating  refrigerant. 
In  order  to  account  for  this,  the  following  procedure  was  adopted  for  both 
one  and  three  inch  sections.  Calculations  are  shown  below  for  the  one 
inch  section.  Heat  transfer  from  the  unheated  end  is  calculated  as  heat 
from  the  base  of  the  fin: 

qf  -    [  (hb)   (p)  (km)   Uc)  ]  °-5  (6b)  (tanh  [  (n)  (Lc)  ]  ) 
where 

P  -  *    (D0) 

-  7T  (.01415)  m 

-  .04445  m 
now 

/^  =  V(D20  -  D2,-) 

-  tt/4(.014152  -  .012702) 
■=  3.0578  x  10"5  m2 


134 


The  corrected  length  of  unenhanced  surface  at  the  end  was  calculated  as 
follows 

Lc  -  Lu  +  (t/2) 

=  0.0254  +  [(0. 01415-0. 0127)/2] 

=  0.0258  m 

hb  is  the  natural  convection  heat  transfer  coefficient  of  the  fin  like 
ends  and  was  calculated  by  using  Churchill-Chu  [Ref.  22]  correlation  for 
natural  convection  for  a  smooth  cylinder,  as  modified  by  Pulido  [Ref. 27]. 


[g(P)  (Dp)   (6b)  (tanh(i2Lc))]   /6 

—  0.6  +.3  87 c 

Dn  r,.,  ._55_919/161*/27 

Pi 


hb  -  -£-  [0.6+.387 .,„      -]2 

D°  [l+[    -559]9/16] 


where 


(V(p)  ■■» 

Therefore  an  iterative  technique  was  necessary  to  calculate  h^,.  The 
iterative  technique  used  was  to  assume  h^,  was  190  W/m2K  and  continue  the 
iteration  until  successive  values  are  within  0.001  of  each  other.  The 
fluid  physical  properties  are  calculated  at  the  vapor  mean  film 
temperature,  given  by  the  following  equation. 


Tfilm 


Tsatc   +  Two 
2 
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„  2.28   +  9.98 

1  film   "  o 


Tfiim  -  6  .  13  °  C  -  279  .  13  °  K 


For  R-114,    the   physical  properties   are   given   in  the  program  by: 
Dynamic   viscosity.    Tf,im   in  °K 

H  =   exp[-4.4636   +    (1011 .47/Tf ilm)  ]xl0-3 

M  =  430.927xl0"6  kg/m  s 
Specific   heat.    Tfjim   in  °K 

Cp  =    [0.40188  +   1.65007xl0-3(Tfilm)    +   1 .  51494xlO_6(Tf  Um2)    - 
6 .  67853xl0-10(TfUm3)  ]xl03 

Cp  =   966.31  J/kgK 
Density.    Tfilm   in  °K 

p   -   16.0184533    (36.32   +  61 .  146414»|»1/3   + 
17  .476838\|r1/2   +    1  .  119828\|r2) 

where 

[1.8 (Tfilm) ] 


i|/  -   1 


and 


753 .95 


<|r   -    .333 


p   -   1512.09    -^2 
m3 


Thermal    conductivity   of  R-114,    Tfjlm   in  °C 
k  =  0.071  -(0.000261)  (Tfilm) 
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k  =   6.936xl(T2  W/mK 

Prandtl   Number 

Pr  =    [(Cp)/i]/k 
Pr  =   6.003 

Thermal    Expansion   Coefficient 
P  -  -    (l/p)  (Ap/Ar) 


p279.03  -   1512.395   iSf 


P279.23     "     1511.824       *f 

/n3 


P   -   -(1/1512.395)  [(.571)/(0.2)] 


P   -   1.89X10-3    (1/J0 


Kinematic  viscosity 


v  -  JL 

P 


_     430.927x10"° 


1512.09 


v   -  2.849X20"7  m2/ s 


Thermal   Diffusivity 
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a-         * 


(P)CC 


6  .936X20-2 
a  — 


(1512.09)966.31 


o  -  4  .147x10-*  m2/s 
Knowing  the  above  properties,  the  heat-transfer  coefficient  h^  can  be 

obtained  by  iteration 


hb  =  362.57  W/m2K 


Knowing  this  we  can  calculate  n 


n  _  r  (362.57)  (44.45X10'3)  ,  os 
(402.28)  (3.0578x10"*) 


n   -  36.19 

then  we   can  obtain   qf 
q,  -    [  (hb)  (p)  (kcu)  (Ac)  ]  °-5  (6b)  (tanh  [  In)  (Lc)  ]  ) 

qf  -    [(362.57)  (.04445)  (402.28)  (3  .  0578X10'5)  ]  05 
(7.7  0)  (tanh[ (36 .19  6)  (0.0258) ] ) 

g,  -  2.51  W 

The  corresponding  results  for  the  three  inch  section  are 
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hb  =  289.47  W/m2K 
qf  =  1.24  W 


Therefore,  the  heat  transfer  through  the  heated  length  of  the  tube  is 

qs  -  q  -  qf  (1  inch  section)  -  qf  (3  inch  section) 

qs  =  (672.19  -  1.24  -  2.51)  W 

qs  =  668.35  W 

and  the  heat  flux  and  the  heat  transfer  coefficient  are  as  follows 

q"  -  qsAs 

=    qs/((w)(D0)(L)) 
=    (668. 35)/((tt)(0.01415)(. 2032)) 
=   7.398x10*   W/m2 
and   finally   the   heat    transfer   coefficient 


A8{Tvo  -  Tsat) 


.    m  668.35 


9  .033xl<T3(7  .70) 


h  -  9.6  0SX103  W/m2K 
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APPENDIX  C:  UNCERTAINTY  ANALYSIS 


The  same  data  run  (TBD1005)  was  chosen  for  the  uncertainty  analysis. 
Therefore,  the  measured  and  calculated  parameters  found  in  the  sample 
calculation  were  used  in  this  section.  The  unceratainty  analysis 
performed  was  for  a  high  heat  flux,  but  the  procedure  could  be  performed 
at  any  heat  flux  to  determine  the  uncertainty  bands.  All  uncertainties 
are  presented  as  a  percentage  of  the  calculated  parameter.  The 
uncertainty  associated  with  the  experimental  parameters  is  calculated  from 
the  equation  suggested  by  Kline  and  McClintock  [Ref.  35].   For  example: 

R  =  R(xi.x2> >xn) 

then 

bR  "  [("lr8Xi)2+(Jr6x2)2+ +  <-^-***)2]0-5 

where 

6R     =  uncertainty  of  the  desired  dependant  variable 

xn  -  measured  variables 

6xn  ■=  uncertainty  in  measured  variables 
The  boiling  heat-transfer  coefficient  is  given  by 


AB(Tvo-Tsat) 
where 


qr[ln(i>] 

T      -  T — 

"°    Vi      2*(/cr„)  (L) 
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In  the  above  equation,  the  second  term  on  the  right  hand  side  is 
usually  called  the  Fourier  heat-transfer  conduction  term.  If  we  define 
this  as 


QT[ln(^2-)] 

4> 


2n(kcv)  (L) 


and 


6b   -  Tvo-Tsati 


With  this  notation,  the  uncertainty  in  the  heat-transfer  coefficient 
is  obtained  using  the  following  equation. 


•£.  [(^)2+(^£)2+(J8I^)2+(M)2+(l^)2]o.! 
h  q  As  db  db  6b 

where 

q=  Vxl 

q=  V(V)xI(V)x60(V)xl(A/V) 

and  the  uncertainty  is 


M  -  t(^)2+(^)2]o.5 

q  V  I 

The  accuracy  in  the  voltage  and  current  sensors  are  as  follows 
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6  Vas   -  ±0.05  V  bAas  -   ±0.025  A 


From  the  sample  calculation  section 

Vas  =  3.189  V  Aas  -  3.513  V 

Therefore, 


_8g   _    .     bVas  ,2+  ^  bAas  ,2,  0.5 
q  Vas  Aas 


M   -    \(    °-05    )2  +  (  0-025)21o.5 
q  3.189  3.513 

A2   -  i.7  2  percent 
Q 

Calculation  of   the  uncertainty  of   the   surface   area   is   as   follows 
As   =   tt(D0)(L) 


A.  D'  L 


Knowing  the  dimensions  of  the  tube  from  the  manufacturer  and  estimated 
inaccuraccies  from  work  shop  tools  and  human  error,  the  unceratinty  was 
calculated. 
Dimensions 

D0  =  14.15  mm  L  =  203.2  mm 

Inaccuracies  in  measurements 
<5D0=  0.1  mm  6L  =  0.2  mm 


Uncertainty  analysis  performed 
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A*  Dr,  L 


M±   _    [(     0,1     }»+(     0^2     ,'jo.. 

Ae  14.15  203.2 


5<^ 
-   -  0.713  5  percent 


The   uncertainty  calculation  for   the   Fourier  conduction  term  given  below 

<t>  or  kcu  L 

kcu  was    calcluated  using 

kcu  -  434.0-  [0.112  (Tvl)] 

kcu  -  434.0- [0.112(283.25)] 
km  -  4  02.28   W/mK 

and   its   uncertainty 

bkcu  -    [  (0.112  (5T,,))2]0-5 

5TH1-  and  <5Tsat  are  obtained  using  uncertainties  in  the  thermocouple 
readings.  Average  wall  inside  temperature  TH1-  was  obtained  taking  the 
average  of  six  thermocouple  readings  inside  the  tube  wall.  The 
uncertainty  associated  with  this  variable  is 

bTwi   -  [6  (  ^  -  V1  )    ]°-5 
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where  6Twj  for  each  thermocouple  was  obtained  by  taking  the  difference 
between  the  measured  wall  temperature  and  the  average  wall  temperature. 
Using  this  method,  it  has  been  attempted  to  try  and  take  into  acount  some 
the  uncertainty  introduced  by  the  fabrication  procedure  for  the  tube  (ie. 
air  gap)  .  For  this  particular  heat  flux  the  following  6TU\  were  found  for 
tube  thermocouples  1,  2,  3,  4,  5,  and  6. 


bfwi  -  t(-^l)2+(-^i)2+(-5#-)2+(-^)2+(-^)2+(-^)2]°-5 

6  6         6         6         6  6 


&Twl   -    .274°  C 


The  uncertainty  level  for  all  remaining  thermocouple  readings  (ie. 
excluding  those  in  the  tube  wall  which  was  considered  to  have  a  higher 
uncertainty)  was  estimated  to  be  +/-  0.1  °C  corresponding  to  an  emf  of 
approximately  4  /xV.  Saturation  temperature  was  obtained  by  taking  the 
average  of  two  thermocouple  readings  and  the  uncertainty  in  this 
temperature  was  calculated  from  the  following  equation. 


6Tsat   -  [2(— £)  ]05 


bTsat   -  [2(-^)2]05 


8  Tsa t  -   0 . 07  °  C 

Knowing  the  uncertainty  in  the  temperatures,  we  can  now  calculate  the 
uncertaines  in  the  following: 
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bkcl]  -    [  (0.112  (6rvi))2]0-5 
bkcu  -    [(0. 112(283. 25)>2] 


bkcu  -  31.7  24   W/mK 


Now  we  can  calculate  the  uncertainty  in  the  Fourier  conduction  term 


<t>       q  kcu  L 


-M  -    [(0.0172)2+(  l\-12*V+{^±-)>V-s 
4>  402.28      203.2 


— ?-   -  8.07  2  percent 


from  the  sample  calculations  we  know  that 
<$>  -   0.2619°  C  6$  -  0.2679(0.08072)  -  0.0216°  C 

therefore : 


^1  .  [(^)2+(^)2.(^i)2+(M)2+(A^£i)2]o.s 
h  q  as  eb     &b  eb 


M    -  [(0.0172)^(0.07135)2+(A4i)2+(-^|^)2+(^4J)2]°-5 
n  7.70       7.70       7.70 


-^  -  7  .42  percent 

/2 


Finally  the  calculation  of  wall  superheat  temperature 


e*  -  rw  -  rsat 


5eh   _    r  (  S  r*°  ) 2  +    (  ftJgat)2]o.5 


eb  eb  eb 
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&Qb    _    r,  0.274  ,2         ,   .07  >21o.5 
0^  Ll      7.7      '  ^  7.7  '    J 


-  3  .  67  percent 


6. 


Table  4  shows  the  results  of  the  uncertainty  analysis  performed.  The  high 
and  low  heat  flux  correspond  to  the  approximate  values  of  7 . 5xlOA  W/m2  and 
lxlO3  W/m2  respectively.  Note  that  the  highest  uncertainty  (over  50%)  is 
in  the  wall  superheat  at  very  low  heat  flux  (1  kW/m2)  .  This  is  due  to  the 
very  low  measured  value  of  wall  superheat  (0.54  °C)  which  can  not  be 
accurately  measured.  Thus  higher  uncertainty  occurs  at  very  low  heat 
fluxes.  However,  once  the  wall  superheat  gets  higher  (at  higher  heat 
fluxes)  the  uncertainty  in  wall  superheat  decreases  significantly  (to 
about  4%)  indicative  of  the  fact  that  the  measure  wall  temperature  is 
relatively  more  accurate. 
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Table  4.  UNCERTAINTY  ANALYSIS  RESULTS 


VARIABLE 

HIGH  HEAT   FLUX 

LOW  HEAT   FLUX 

e> 

7.7 

0.54 

Tul 

10.25 

2.80 

TSat 

2.24 

2.20 

6Vas/Vas 

1.57% 

14.6% 

5Aas/Aas 

0.712% 

7.9% 

5q/q 

1.72% 

16.6% 

5D0/D0 

0.707% 

0.707% 

5L/L 

0.098% 

0.098% 

S  As/As 

0.714% 

0.714% 

0  Kcu/K.cu 

7.89% 

7.89% 

66^/6^ 

3.67% 

52.37% 

5h/h 

7.42% 

28.93% 
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APPENDIX  D:  OPERATING  PROCEDURE 


A.    SYSTEM  STARTUP 

1.  Power  to  the  28  kW  (8  ton)  refrigeration  unit  is  provided  by  the 
breakers  located  in  the  main  distribution  panel  located  in  the  laboratory. 
These  breakers  were  never  secured.  However,  if  power  to  this  panel  was 
lost,  then  these  breakers  must  be  reset. 

2.  Turn  the  switch  on  the  refrigeration  unit  control  panel,  located 
in  front  of  the  refrigeration  unit  to  the  "auto"  position  after  passing 
through  "on"  position.  This  switch  is  always  left  on,  unless  unit  was 
taken  down  for  long  repairs . 

3.  Push  the  start  button  in  the  control  box  for  the  recirculation 
pump.  This  control  box  is  located  on  the  bulkhead  above  the  recirculation 
pump  in  the  outside  area  adjacent  to  the  refrigeration  unit. 

4.  Set  the  desired  temperature  on  the  roughly  graduated  Fahrenheit 
scale  on  the  control  panel  thermostat.  It  requires  approximately  one  hour 
to  chill  the  sump  to  -  15  °C .  The  thermometer  in  the  ethylene  glycol/water 
mixture  (sump)  must  be  monitored  to  ensure  the  desired  sump  temperature  is 
attained  and  maintained.  Slight  adjustments  in  the  refrigeration  unit 
thermostat  can  be  expected  due  to  the  coarseness  of  its  scale. 

5.  Energize  the  desired  pumps  by  switching  on  the  breakers  in  the  main 
distribution  power  panel.  Once  the  power  is  energized  in  the  main  power 
panel  and  after  ensuring  the  pump  suction  valves  are  open,  turn  on  the 
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pump  motors  by  pushing  down  on  the  arm  of  the  appropriate  breaker  box  for 
the  pumps  located  on  the  bulkhead  next  to  the  ethylene  glycol/water  sump. 
The  pumps  are  marked  "auxiliary  condenser"  (Pump  #2)  and  "instrumented 
tube  condenser"  (Pump  #1),  respectively  in  the  breaker  box.  Flow  in  the 
auxiliary  condensate  coils  can  be  controlled  with  the  individual  globe 
valves  located  at  the  coil  penetrations  on  the  apparatus.  The  auxiliary 
condenser  coils  will  produce  the  fastest  adjustments  to  the  system 
pressure . 

6.  Energize  the  heater  variac(s)  desired  by  switching  on  the  breakers 
(Bundle,  and  Simulation  (for  test  7  only))  in  the  main  distribution  power 
panel  and  individual  breakers  for  each  of  these  in  the  power  distribution 
box  (near  apparatus) . 

7.  After  ensuring  that  the  breakers  for  the  heated  tubes  desired  are 
in  the  "on"  position,  follow  the  experimental  procedures  for  normal 
operation  outlined  in  Chapter  IV. 

B.    SYSTEM  SHUTDOWN 

1.  Turn  all  variacs  to  the  zero  position  and  switch  off  all  breakers 
in  the  power  panels. 

2.  If  apparatus  will  not  be  operated  for  an  extended  period,  turn  the 
switch  on  the  refrigeration  control  panel  to  the  "off"  position  after 
passing  through  "on". 

3.  Allow  the  recirculation  pump  to  operate  for  at  least  five  minutes 
after  switching  off  the  refrigeration  pump  unit  to  dissipate  any  back 
pressure  in  the  system. 
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4.  Turn  the  breakers  for  the  pumps  to  the  off  position  at  the  switch 
boxes,  and  then  secure  the  power  at  the  main  distribution  power  panel. 

C.    EMERGENCY  SHUTDOWN 

1.  Secure  all  power  at  the  main  distribution  power  panel 

2 .  Evacuate  building 

3.  Call  Fire  Department 
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APPENDIX  E:  PROGRAM  DRP4RH 


10001  FILE  NAME   0RP4RH 

10041  DATE        November  22,  1388 

1006'  REVISED    FEB  1992  (R.  HAAS  > 

1012" 

1(216  BEEP 

1023  PRINTER  IS  1 

1024  Idp-0 

10231 

1032  PRINT  USING  "4X, ""Select  option  default  is  0  

1025  PRINT  USING  "6X,""0  Taking  data  or  re-process ing  previous  data 

1040  PRINT  USING  " 6 X ,  • " 1  Plotting  data  on  Log-Log  

1044  PRINT  USING  "6X,""2  Plotting  data  on  Linear 

1048  PRINT  USING  "BX,""3  Purge""" 

1052  PRINT  USIN6  *SX,""4  Fixup""* 

1056  PRINT  USING  "GX,"S  Move""" 

1060  PRINT  USING  "6X,""fi  Comb""" 

1064  PRINT  USING  "6X,'"7  Read  Plot""" 

10681 

1072'  IDP  IS  A  PROGRAM  VARIABLE  TO  SELECT  A  SUBROUTINE 

107S  INPUT  Idp 

1080  IF  Idp-0  THEN  CALL  Main 

1084  IF  Idp-1  THEN  CALL  Plot 

1088  IF  Idp-2  THEN  CALL  Plin 

1092  IF  Idp-3  THEN  CALL  Purg 

1096  IF  Idp-4  THEN  CALL  Fixup 

1100  IF  Idp-S  THEN  CALL  Move 

1104  IF  Idp-5  THEN  CALL  Comb 

1108  IF  Idp-7  THEN  CALL  Readplot 

1112  ENO 

11161 

1120  SUB  Main 

11241  ICAL-THERMOCOUPLE  CALIBRATION 

1128  COM  /Cc/  C(7) 

1132  DIM   Emf<  35  >,T<  3S),0la(6 ),D2a<6  >,Dia(6  ),Doa<6  >,La(6 ),Lua(S  >,Kcua<6  ),Et( 19  ). 

Ldtc(20),Uolt<2  >,Amp(  11  ),Tua(  5>,Twl5),Theta(5),Thetab<5),Q<5>,01<S),Qdp(5) 

1136  DIM   Htube(5  ),Tn(5),Tp(S  ) 

11401 

11441  THERMOCOUPLE  ARRAY  <C<  ))  INITIALIZATION 

1148  DATA  0.10086091,25727.94369,-767345.9295,78025595.81 

1152  DATA  -9247486589, 5. 97688E+11, -2. 56192E+13, 3. 94078E+14 

1156  REAO  C(«  ) 

1160> 

1164  1  PRINT  HEAOER  AND  INITIALIZE  TIME  CLOCK 

1168  PRINTER  IS  701 

1172  BEEP 

1176  INPUT  "ENTER  MONTH,  DATE  AND  TIME  ( MM : OD  HH' MM  SS ) " , DateS 

U80i  OUTPUT  OIRECTEO  TO  OATA  AQUISITION  SYSTEM  (HP  3497A  ) 

1184  OUTPUT  709. "TO". DateS 

1188  OUTPUT  709. "TD" 

1192  ENTER  709;Oate$ 

1196  PRINT 

1200  PRINT  "  Month,  date  and  time   ".DateS 

1204  PRINT 

1208  PRINT  USING  "10X,""NOTE   Program  name    0RP4RH 

1212  BEEP 

1216' 

1220'  ON  IS  THE  VARIABLE  FOR  DISC  NUMBER  FOR  RECORD  KEEPING  ONLY 

1224  INPUT  "ENTER  DISK  NUMBER", Dn 

1228  PRINT  USING  " 16X, " "Disk  number   -  "",ZZ".0n 
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1232  BEEF 

123E  l*-0 

1240  INPUT  "ENTER  INPUT  MODE  ( 8-34 97A,  1  -FILE  )  0-DEFAULT  '  ,  In 

1244  ' 

12481  INPUT  MODE  ZERO  IS  FROM  THE  DAT*  AQUISITION  SYSTEM 

1252  IF  Im-e  THEN 

1256  BEEP 

1260  INPUT  "6IVE  A  NAME  FOR  THE  RAW  DATA  FILE" ,D2f l leS 

1264  PRINT  USING  " 16X ,  ""File  name   " " , 1 4A " , D2 f i leS 

1268' 

1272  1  CREATE  BDAT  FILE  ON  THE  MASS  STORAGE  MEDIA 

1276  CREATE  BDAT  D2fileS,60 

1280'  CREATE  AN  INPUT/OUTPUT  LINK  TO  OPEN  FILES 

1284  ASSIGN  gcile2  TO  D2fileS 

1288' 

1292!  CREATE  DUMMY  FILE  UNTIL  Nrun  KNOWN 

1296  Dlf ile$-"OUMMY" 

1300  CREATE  BDAT  Dlfilel,60 

1304  ASSI6N  CFilel  TO  DlfileS 

1308  OUTPUT  SFileliDateS 

13121 

1316'  CREATE  A  PLOT  FILE 

1320  BEEP 

1324  INPUT  "6IVE  A  NAME  FOR  THE  PLOT  FILE'.PfileS 

1328  CREATE  BDAT  PfileS,30 

1332  ASSI6N  »Plot  TO  PfileS 

1336  BEEP 

1340' 

1344'  IDTC  -  NUMBER  (TOTAL)  OF  DEFECTIVE  THERMOCOUPLES 

1348  INPUT  "ENTER  NUMBER  OF  DEFECTIVE  TCS  ( 0-DEFAULT  ) ' , Idtc 

1352'  LDTC  -  LOCATION  OF  OEFECTIVE  THERMOCOUPLE 

1356' 

1360  IF  Idtc-0  THEN 

1364  PRINT  USIN6  "lBX,""No  defective  TCs  exist**" 

1368  ELSE 

1372  PRINT  USING  " 1BX ," "Defective  Thermocouple*  Indicated  by  -99.99""' 

1376  END  IF 

1380' 

1384  BEEP 

1388'  DEFECTIVE  THERMOCOUPLES  MAY  BE  IN  CHANNELS  40-E9 

1392'  THEMOCOUPLES  ARE  ENTERED  AS  DEFECTIVE  BY  COMPUTER  CHANNEL  NO. 

1396'  JDTC-COUNTER  IN  LOOP  FOR  DEFECTIVE  THERMOCOUPLES 

1400' 

1404  IF  ldtc>0  THEN 

1408  FOR  Jdtc-0  TO  Idtc-1 

1412  INPUT  "ENTEF  DEFECTIVE  TC  LOCATION  (BY  COMPUTER  CHANNEL  NUMBER) 

" ,Ldtc( Jdtc  ) 

1416  BEEP 

1420  NEXT  Jdtc 

1424  END  IF 

1428  PRINTER  IS  701 

1432  OUTPUT  GFilel;Ldtc(«  ) 

1436  i 

1440  i    Im-1  option  (THIS  OPTION  ALLOWS  DATA  ENTRY  UITH  DATA  FILE) 

1444  ELSE 

1448  BEEP 

1452  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" ,D2f i leS 

145E  PRINT  USING  "16X,""File  nane   " " , 14ft" ; D2f l le$ 

1460  ASSIGN  eFile2  TO  D2fiieS 

1454  ENTER  @File2.Nrun 
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M68 

u?: 

U7S 

1480 

M84 

1488 

1492 

149G 

1500 

1504 

1508 

1512 

1516 

1520 

1524 

1528 

1532 

1536' 

1540 

1544 

1548 

1552 

1556 

1560 

1564 

1563 

1572 

1576' 

1580 

1584 

1588 

1592 

1596' 

1600 

1604' 

1608 

1512 

1616 

1620 

1624 

1628' 

1632 

1633 

1634 

(DEF 
1B35 
1636 
1637 
1639' 
1640 

1644 

1648 

16S2 

1653 

1655 

1656 

pos 

1657 

1658 

1660' 


ENTER  SFileZ.DoldS.Ldtc 

SEEP 

INPUT  "GIVE  ft  NAME  FOR  PLOT  FILE 

CREATE  SCAT  PfileS,30 

ASSIGN  SPlot  TO  P'llei 


I tt ,Bop , Nht , Nat o.Nrt , Cor r 
.PfileS 


PRINT  USING  "16< 
SEEP 

PRINTER  IS  1 
PRINT  USING  "4X, 
PRINT  USING  "6< , 
PRINT  USING  "SX, 
PRINT  USING 
PRINT  USING 
INPUT  Itt 

END  IF 

IF  Im-1  THEN  GOTO  1763 

PRINTER  IS  1 


"SX, 
"BX, 


This  data  set  taken  on 


"SELECT  TU6E  TYPE 


14A"  .DoldS 


CMOOTH" " " 
FINNEDf  19/ IN  ) 
HIGH  FLUX  """ 
TUR60-B  ""• 


IF  In-0  THEN 

PRINT  USING  "4X 

PRINT  USING  "6X 

PRINT  USING  "6X 

PRINT  USING  *SX 

PRINT  USING  "6X 

PRINT  USING  "SX 

PRINT  USING  "SX 

PRINT  USING  "SX 

ITT-TUBE  TYPE 

INPUT  Itt 

OUTPUT  9Filel  .  Itt 
END  IF 

PRINTER  IS  701 
Itt-2 
PRINT  USING  '16X,""Tube  Type 


Select  tube  t  ype  "  * " 

0  Smooth  '"" 

1  FINNED  19/IN  (DEFAULT)'"' 

2  HI6H  FLUX""" 


TURBO-B' 
GROWTH" ' 
GROWTH"' 
GROWTH"' 


,  DO" . Itt 


BEEP 

Bop-0 

INPUT  "ENTER  BULK  OIL  X   < OEFAULT-0X )  ",Bop 

OUTPUT  SFilel.Bop 

PRINT  USING  •lSX.-'Bulk  oilX-"" , DO". Boo 

BEEP 
Ipo-1 

INPUT  "ENTER  POOL  HEIGHT  ABOVE  TOP  TUBE  (0-LESS  THEN  5cm,  l-5cm  OR  GREATER 
AULT)",Ipo 
OUTPUT  9Filel , Ipo 

PRINT  USING  "16X,""Pool  height-" *, DO"  .  Ipo 
BEEP 

NHT-NUMBER  OF  HEATED  TUBES 
Nht-5 

INPUT  "Enter  number  of  heated  instrumented  tubes* defaui t«5  )", Nht 
OUTPUT  ©Fi lei .Nht 

PRINT  USING  "16X  ,"  "Number  of  heated  instrumented  tubes" ""  ,00 ". Nht 
BEEP 
Ipos»  1 
IF  Nht-1  THEN  INPUT  "WHICH  POSITION  IS  THIS  SINGLE  TUBE  (1  TO  5   DEF=>1)",I 

Ip05"Ip0S-l 

BEEP 
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166- '  Natp"Number  of  active  aumm>  pairs 

16EE  Natp-0 

1E7Z  ]NDUT  "Enter  number  of  active  dummy  pairs  i De f aul t «0  ' '  , Na tc 

1676  OUTPUT  *JPi  lei .NatD 

1680  PRINT  USING  " 16X, " "Number  of  active  dummv  pairs- "-, DC" . Natp 

1684  BEEP 

1668' 

169Z  i  NRT-NUMSEF.  OP  ADDED  HEATED  TUBES  TO  ENHANCE  BUNDLE  EXPECT 

1E9E  N-t-0 

1700  INPUT  "Enter  number  of  added  heated  tubes  from  simulation  heaters < Defaul  t" 

0 ) " , Nr  t 

1704  OUTPUT  BFl lei .Nrt 

1708  PRINT  US1N6  " 16X , ""Number  of  added  heated  tubes(from  simulation  heaters*'* 

"  ,DD" ;Nrt 

1712  BEEP 

17161 

1720'  CORR  IS  CORRECTION  FOR  INSTRUMENTED  TUBE  HEIGHT 

1724  Corr-0 

1728  INPUT  "WANT  TO  CORRECT  TSAT  FOR  TUBE  HEI6HT  ( 0-YES< DEFAULT  ), 1-NO  )", Com 

1732       IF  Corr-0  THEN  PRINT  USING  "16X,""TSAT  is  corrected  instrumented  heat 

ed  tube  height"" 

1736       IF  Corr-1  THEN  PRINT  USIN6  "1EX,""TSAT  is  NOT  corrected  for  instrumen 

ted  heated  tube  height""" 

1740  OUTPUT  flFilel.Corr 

1744  BEEP 

1748'  ILQU-INPUT  MODE;  LIQUID,  UAP0R,0R  LIQUID  UAPOR  AVERAGE 

1752  Ilqv-0 

1756  INPUT  "SELECT  (  0-LIQ(  defaul t  ), 1-UAP , 2-( LIQ+UAP  )/2  )", I  lqv 

17601 

1764'  DlA«Diameter  at  thermocouple  positions  (meters) 

1768  DATA  .0122,0.0098,0.0106,0.0116,0,0,0 

1772  READ  Dla<»  ) 

1776  Dl-Dla(Itt) 

17801 

17841  D2*Diameter  to  base  of  fins  (outside  dia  for  smooth  )( met ers  ) 

1788  DATA  .0158,0.0125,0.0156,0.01415,0,0,0 

1792  READ  D2a( •  > 

179E  D2-D2a(Itt) 

1800' 

1804'  Di"Inside  diameter  of  unenhanced  ends  (meters) 

180S  DATA  .0132,0. 0109, 0.0116,0. 0127, 0,e,0 

1812  READ  DiaC  ) 

iei6  Di-Dia(Itt) 

18201 

1824'  Do'Outside  diameter  of  unenhanced  ends  (meters) 

182e  DATA  .015675,0.0125,0.015875.0.01415,0,0,0 

1832  READ  Doa<«  ) 

1636  Do-Doa(Itt) 

1840' 

1844'  L"Length  of  enhanced  surface  (meters) 

1848  DATA  .2032. .2032,  .2032, .2032, .2032, .2032, .2032 

ieS2  READ  La(«  ) 

185E  L-La(Itt) 

I860' 

1864'  Lu«C0RRECTED  Length  of  unenhanced  surface  at  the  end-  (rETEPE 

18S8'LU=LPIN  +  THICKNESS/2 

1672  DATA  .0261, .0254, .0264. C. 0258, C,C,e 

1876  READ  Luai ♦  ) 

1880  Lu=Luai Itt  ) 

1864' 
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1 68 9 •  L V=correc ted  length  o'  ~    inch  finneC  lite  end 

16S2  DIM  ^vatE  ) 

199E  DATA  . C76S , . C7EZ ,  .  0773 , G .C76S  ,  0 .0  ,  0 

19P.C  READ  Lvei  •  > 

1904  Lv-Lva<  Itt  ) 

1908'  Kcue'Therns  1  Conductivitv  of  tube 

1912 '  DATA  401 ,G,C,C,e,C,G 

1916 i  READ  Kcuai  •  ' 

1920'  Kcu-Kcua' Itt  ) 

1924  A-PI»(Do*2-Di"2>/4 

1928  F-PI»Do 

1932  J-l 

1936  5x-0 

1940  Sy-0 

1944   Sxs-0 

1946   S*y«0 

1952  Repeat  :  i 

19561 

i960   IF  In-0  THEN 

19641      Dtld^desired  teriperature  of  liquid 

1968i     Dtld-47.S     IR-113 

1969      Dtld-2.2      'P-114 

1972      Ido-2 

1976      ON  KEY  0,15  RECOVER  1952 

1980      PRINTER  IS  1 

1984 

1988 

1992 

1996 

2000 

2004  1 

2008 

2012 

20161 

2020      BEEP 

2024'      Set  default  value  for  input 

2028      IP  Ido)2  THEN  ido-2 

2032'      Take  data  option 

2035      IP  Ido-0  THEN  2440 

2040  i 

2044i  LOOP  TO  SET  HEAT  FLUX   (FOR  TOP  INSTRUMENTED  TUBE) 

2048      IF  Ido-1  THEN 

2052      Dqdp-100000 

205E  PRINT  USING  "4X,""0dp  OOP sin       Nrt      Gdpau> 

Qtof"' 
2060  PRINT  USING  '4X,  _"(W/n-2  )        (W/n*2)  (U/n"2> 

(  U  ) "  • " 
2064  Err-1 

2066'  Reset, read  channel  25-30 .autonat ic  scaling 

2072 i  Channel  25-au>  anps,2E*sin  anps,27«inst  vol  t  s  ,  26-s  in,  volts, 29-au> 

vol t s , 30-34=inst  anps 
2076  OUTPUT  709. "AR  AF25  AL34  VR5" 

2080  FOR  1-10  TO  11 

2084  OUTPUT  70S. "AS  SA" 

2088  ENTEF  70S. Anp'  I  ) 

2092  NEXT  I 

209E  FOF  I-e  TO  2 

2100  OUTPUT  70S.  "AS  SA" 

2104  ENTER  709.Uolt( I  ) 

2106  NEXT  I 


PRINT  USING 

"4X, 

■■SELECT  OPTION  "'" 

PRINT  USING 

•6X, 

••0-TAKE  DATA--" 

PRINT  USING 

"EX, 

-"1-5ET  HEAT  FLUX"" 

PRINT  USING 

"BX, 

■-2-SET  Tsat  (OEFAULT  SET 

PRINT  USING 

"4X, 

--NOTE^  KEY  0  -  ESCAPE--- 

I do-des ired 

opt l on 

BEEP 

INPUT  Ido 
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Z1CS  FOR  j-C  TC  4 

:: ::  outcut  709.  -as  sa" 

211E  ENTEF  70S.Amp< I  ) 

2117  NEXT  I 

2120'  Calculate  actual  heat  flu> 

2124  Q( 0 '-EC'Volt < 0 >*Amp( Ipos  ) 

212E  Oap<C >-O(0  )/(PI»D2»L  ) 

2132  Osim-E0«20«Uolt( 1  >«Amp( 11  ) 

213E  Qdpsim-Qsim/(FI*C2«  .2032*3  > 

2140  Oeu»-B0«20»Volt (2 >»Anp( 10  ) 

2144  Qdpau>-Qaux/<  PI • .0160* . 1778*4  ) 

2146  Qtot-O<0  )«Nht+Osim+Gaux 

2152  Nrt-Qdpsim/Qdp<0  ) 

2156  IF  ABSt Aqdp-Dqdp  )>Err  THEN 

2160  IF  Aqdp>Dqdp  THEN 

2164  BEEP  4000, .2 

2168  ELSE 

2172  BEEP  250, .2 

2176  END  IF 

2180  IF  NrK.l  THEN  Nrt-0 

2184  IF  Qdpaux<100  THEN  Qdpaux-0 

2188  IF  Qdpsim<100  THEN  Qdpsim-0 

2192  PRINT  USING  "4X , 2( M2 . 3DE , 2X  ) , 2X , ( MDD .DD  ) , 2X , 2( MZ . 3DE, 2X  )  * . Qdp 

<0),Qdpsim,Nrt,Qdpau*,Otot 

2196  WAIT  2 

2200  60T0  2076 

2204  END  IF 

2208      END  IF 

22121 

2216'     LOOP  TO  SET  Tsat 

2220      IF  Ido«2  THEN 

2224  IF  Ikdt-1  THEN  2240 

2228  BEEP 

22321  INPUT  "ENTER  DESIRED  Tsat  < DEFAULT-47 .5  C  -  R-l 13  ) " . Dt Id 

2233  INPUT  "ENTER  DESIRED  Tsat  (DEFAULT-2.2  C  -  R-l  14  ) " , Dt  Id 

2236  Hdt-1 

2240  Oidl-0 

2244  Oid2«0 

2246  Nn«l 

2252  Nrs«Nr,  MOD  15 

2256  Nn-Nn+1 

2260  IF  Nrj-l  THEN 

2264  PRINT  USING  -4X,'-   DTsat      Tldl      Tld2      Tlbb 

Tvat  Tvab      Tlav 

2268  END  IF 

2272 '  Read  thermocouple  voltages  for  vapor,  liquid 

2275  OUTPUT  70S."AR  AF0  AL5  UR5" 

22301  Sample  each  thermocouple  20  times  and  report  temp  for  each  the 
rmocouple,  vapor»0,l,2.   liquid-3&4 

2264  FOR  1-0  TC  E 

2286  Sum»0 

2292  OUTPUT  70S. "AS  SA" 

2295  FOR  Ji=l  TO  20 

2300  ENTER  709.Eliq 

2304  Sum=Sum+Eiiq 

2308  NExT  Ji 

2312  Emf! I  >«Sum/20 

2315  T(  I  i^NTvsv!  Emf (  ! 

2320  NEVT  I 

2324  1  Compute  average  temperature  of  liquid 
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:::e 

T]ev-< TC3>+T<4  i  >«.5 

2332  i 

Cocputt  average,  t enpca t w~e    of  vapc 

^  jj  j 

Tvev]-(T<0  KT(  1  )  )/2 

:334 

lvev>Ti  2  ) 

233E 

Tva v-f  T(  0  )  +  T(  1  >  +  T(  2  )  i-  Z 

2340 

IF  ABS(Tlav-DUd) '.I  THEN 

2344 

IF  Tlav>Dtld  THEN 

2348 

BEEP  4000, .2 

2352 

ELSE 

235E 

BEEP  25e, .2 

23E0 

END  IF 

2364 

ELSE 

2366 

IF  ABS<  Tlav-DUd)  .1  THEN 

2372 

IF  Atld'Otld  THEN 

2376 

EEEF  3000, .2 

2380 

E-SE 

2384 

BEEP  800,  .2 

2388 

END  IF 

2392 

END  IF 

2396 

END  IF 

2400 

Errl-Tlav-Oldl 

2404 

Oldl-Tlav 

2408 

Err2-Tvav-01d2 

2412 

01d2-Tvav 

2416 

PRINT  USIN6  "4X,7(MDDD.DD,3X>,;Dtld,T(3),T(4),T(5),Tvavl,Tvav 

V 

2420 

WAIT  2 

2424 

60T0  2252 

2428 

END 

IF 

2432i 

2436' 

TAK 

E  DATA  IF  In-0  LOOP 

2440 

IF 

Ikol-1  THEN  2452 

2444 

BEEF 

2448 

Ikol-1 

2452 

OUTPUT  709. "AR  AF0  AL5  UR5" 

2456 

FOR  1-0  TO  5 

2460 

OUTPUT  709. -AS  SA" 

2464 

Sum-0 

2468 

FOR  J i-i  TO  20 

2472 

ENTER  709. E 

247E 

Sun«Sufi  +  E 

2480 

IF  1>2  THEN  Et(Ji-l  )-E 

2484 

NEXT  Ji 

2488 

Kdl-0 

2492 

IF  1/2  THEN 

2496 

Eave-Sun/20 

2500 

Sun-0. 

2504 

FOP  Jk-0  TO  19 

2508 

IF  ABS(Et(Jl  >-EaveK5.0E-E  THEN 

2512 

Surn«Sun+Et  (  Jk  ) 

2516 

ELSE 

2520 

Kdl-k'dl  +  1 

2524 

END  IF 

2528 

NEXT  Jl 

2532 

IP  12  THEN  FPINT  USING  "4X,""kdl  -  "",DD";Kdl 

2E35 

IF  l  dl  10  THEN 

2540 

BEEP 

2544 

EEEP 

254  8 

PRINT  USING  "4X,""Too  nutr,  scattering  in  date 

pest  date 

set 

.Tie 
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:s5:  goto  isee 

2556  END  IF 

2560  END  IF 

2564  Emf;i  i«6um/i20-KdJ  ! 

2566  NEXT  I 

2572  OUTPUT  709. "AF  AF4C  ALE9  VP5" 

2576  FOF  1-6  TO  35 

2see  output  709. -ae  sa" 

2584  Sum-0 

2586  FOF  J i-l  TO  5 

2592  ENTEF  709, E 

2596  Suf=Suri  +  E 

2600  NEXT  Ji 

2604  Enf( I  >-Sum/5 

2606  NEXT  I 

26121 

2616!  READ  UOLTAGES  ( 27-Inst ,28«Siii,29-Aux  ) 

2620  OUTPUT  709. "AR  AF2"  Al29  UR5" 

2624  FOR  1-0  TO  2 

2628  OUTPUT  709. "AS  SA" 

2632  ENTER  709.Uolt<  I  ) 

2636  NEXT  I 

26401 

26441  READ  CURRENTS  (30-34«Inst  tubes . 3S-39-ACT IUE  Dumny  ) 

2648  OUTPUT  709. -AR  AF30  Al_39  UR5" 

2652  FOR  1-0  TO  9 

2656  OUTPUT  709.  "AS  SA" 

2660  ENTER  709.Amp( I  ) 

2664  NEXT  I 

26681      Read  Current s( 25-Aux  afips,26-Sin  amps) 

2672  OUTPUT  709; "AR  AF25  AL26  URS" 

2676  FOR  1=10  TO  11 

2680  OUTPUT  709;  "AS  SA" 

2684  ENTER  70S.Anp*  I  ) 

2688  NEXT  I 

2692  ELSE 

2696  ENTER  eri le2 , Enf ( •  ) , Uol t ( ♦  ) , Anp( •  ) 

2700      END  IF 

2704  1 

2708'      CONVERT  EMF'S  TO  TEMP , VOLT , CURRENT 

2712      FOR  1-0  TC  35 

2716  T(  I  )-FNTvsv(Enf ( I > ) 

2720  IF  I>5  AND  ldtc>0  THEN 

2724  FOP  Ii-0  TO  Idtc-1 

272e  IF  Ldtct li  1-1-5  +  39  THEN  T(I>— 99.99 

2732  NEXT  Ij 

2736  END  IF 

2740      NEXT  1 

27441      Ntc-nr  of  thermocouples 

2748  Ntc-6 

2749  Ip  Ipos  >0  THEN 

2750  Q(Ipos>-=60«Volt<0  >»Anp<  Ipos  ) 

2751  Tua' Ipos  -0 

2752  Jj-0 

2754  Ndtc-0 

2755  FOF  1=1  TO  Ntc 
275E  Nn*Ipos*6+o+ J j 
275-  Jj-Jj+1 

2755  Ic  AB5< T(Nn)  )  95  THEN 

2760  T(Nn  '  =  -99.99 
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2761  Ndtc«Ndtc+l 

:7g:  else 

Z76Z  Twe( Ipos  >=Twa>  Ipos  >  +  T(Nn  i 

2764  END  IF 

2765  NEXT  I 

2766  Tue'  Ipos  )«Twa< Ipos  >/( 6-Ndtc  ) 

2767  GOTO  2620 
276S  END  IF 

2765  FOR  Ii-e  TC  4 

2770  0-  I i  '-60«Volt ( G  i»flmp( I  l  ) 

27711  Tya-Average  temperature  of  the  wall 

2772  Tuial  Ii  1-0 

2773  Ndtc-0 

2774  FOF  1-1  TO  Ntc 

277G1  Nn  is  counter  in  temp  a-ray,  start  at  E  (this  is  the   first  th 

ernocouple  in  the  tube  bant  ) 

2780  Nn=Ii«6+I+5 

2784  IF  ftBS(T(Nn)  )>99  THEN 

2788  T(Nn)--99.99 

2792  Ndtc-Ndtc+1 

2796  ELSE 

2800  Twadi  )-Twa(Ii  )+T(Nn) 

2804  END  IF 

2808  NEXT  I 

2812  Twa< Ii >-Tua( Ii  )/(6-Ndtc  ) 

2816  NEXT  Ii 

2820  Tlav-< T< 3 >+T(4 >  )/2 

2821  Tvav-T(2) 

28241  Tvav(  T(0  )  +  T(  1  )  +  T(2  )  )/3 

2828' 

2829'  Tlav-T<5) 

2832  Tcu«Tuia(0' 

2836  Kcu-FNKcu(Tcu)      'THERMAL  CONDUCTIVITY  OF  COPPER 

2840'  'IF  CURVE  FIT  NOT  AVAIL  USE  ARRAY  KCU< •  ) 

2844'  FOURIER  CONDUCTION  EQUATION  WITH  CONTACT  RESISTANCE  NE6LECTED 

2848  FOR  1-0  TO  4 

2852  Tu(I  )-Tua( I  )-Q(I  )»L0G(D2/Dl  >/< 2*PI»Kcu»L ) 

2856  IF  Ilqv-0  THEN  Texs-Tlev 

2860  IF  Ilqv-1  THEN  Te*s-Tvav 

2864  IF  Ilqv-2  THEN  Texs-(TIav+T<2  )  )• .5 

2868  IF  Corr-1  THEN  The t ab ( I  >-Tu< I  )-Texs 

2872  IF  Corr-0  THEN  Thet ab (1  )-Tw( I  >-( Texs* . 056  + 1 • . 129  )   'R-114 

2876'  IF  Corr-0  AND  Ipo-1  THEN  The t ab ( I >-Tu( I  )-( Texs  +  .  054+ I « . 144  )   'R-ll 

2877'  IF  Corr-0  AND  Ipo-0  THEN  Thetab ! I  >-Tu( I  )-' Texs- 1 . 078+ . 147« I  )    IR- 

113 

2880  NEXT  I 

2884  1 

2886'  COMPUTE  VARIOUS  PROPERTIES 

2e92  Tf iin-(Tu<e >  +  Texs  )«.5   'FILM  TEMPERATURE 

2896  Rho-FNRho' Tf ilm )        'DENSITY 

2900  Mu-FNMu(Tf i In )          'VISCOSITY 

2904  K-FNK(Tf ilm)            'THERMAL  CONDUCTIVITY 

2908  Cp-FNCp;Tf  ilm  >          'SPECIFIC  HEAT 

2S12  6eta=FN6eta(  Tf  ilfn  )      'THEPMAL  EXPANSION 

2516  Ni-Mu/Rhc               IKINEMATIC  VISCOSITY 

2920  nlpha-K/(Rho*Cp )       'THERMAL  DIFFUEK'ITY 

2924  Rr=Ni/Alphs             'PRANDTL 

2928' 

2932'  COMPUTE  NATURAL-CONVECTI VE  HEAT-TRANSFER  COEFFICIENT 
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:°3B' 

2940 

2944 

294B 

2952 

2956 

29E0 

2964 

2966 

'97'' 

2976 

2980 

2984 

29881 

2992  1 

2996 

3000 

3004 

3006 

3012 

301E 

3020 

3024 

3028 

3032 

3036 

3040 

3044 

3048 

3052 

3056 

3060 

3064  1 

3068' 

3072 

3076 

3080 

3084 

3085' 

3087 

3088 

3089 

3090 

30911 

3092' 

3095 

3100 

3104 

3108 

3112 

Tvav    T 

3116 

3120 

3124 

H 
3128 

2)   (W/m* 
3129 
3130 
3132 


FOF  UNENHANCEC  END<  S  I 

LU'LUB-  ltt  ) 

Hbar=190 
Fe*<Hbar»P/(Kcu»f\))' -5«Lu 

Tanh-FNTanh(Fe  ) 

Theta(  Ipos  >-Thetab( Ipos  >«Tanh'Fe 

Xx-(9.61«Beta»Thetab'Ipos  >«Dc">T  anh'  (Fe«Ni« Alpha  )  v  .  1E6SE"7 

Yy-i  l  +  <  .  559/Pr  r<  9  ■'  IE  )  )"(  e/27  ) 

Hbarc-K/Do»(  .6+ . 387»Xx/Yy  !*2 

IF    AES*  <HDe<--Hbarc  )/Hbarc  )'  .001    THEN 

Hbar-( Hbar+Hbarc  )•  .5 

60TC  2946 
END  IF 

COMPUTE  HEAT  LOSS  RATE  THROUSH  UNENHANCED  END(S) 
Ql(0)«<Thetab(Ipos>»Tanh)»((Hbar«P»Kcu«A;".5> 
Oq-OKO  )  +  0q 
2  =  2+ 1 

IF  Z-l  THEN 
Lu-Lv 
GOTO  2944 
END  IF 
Z-0 

Qipct-Qg/Q< Ipos ) 
Oq-0 

As-PI«D2*L 
FOR  I i-0  TO  4 

OK  Ii  >-Qlpct»Q(Ii  ) 

Odp(  1 1  >-<  0(  1 1  >-QK  Ii  )  >/As 

Htubef Ii  >-Qdp( I  i  >/Thetab( Ii  ) 
NEXT  Ii 
PRINTER  IS  701 

RECORD  TIME  OF  DATA  TAKING 
IF  Im-0  THEN 

OUTPUT  709. -TD" 

ENTER  70S.Tcld$ 
END  IF 

CHURCH1LL/CHU  CORRELATION  FOR  NATURAL  CONVECTION  REGION 
Raa-9.81«Beta«Thetab( Ipos  )•( D2  >* 3 •Rho/t Mu« Alpha ) 
Denon«=(  1  .  +  (  .  559/Pr  ) "  (  9/  16  )  )  *  (  16  /9  ) 
Nuch«(  .6+.387«(Raa/Denot>i  )"(  1/6  )  )"2 
Qch-K«Nuch«Thetab( Ipos  )/(D2  ) 

OUTPUT  DATA  TO  PRINTER 

PRINTER  IS  701 

PRINT 

PRINT  USING  "10X, "-Date  Set  Number 

PRINT 

PRINT  USING  "10X,""   Tvl      Tv2 
Idav  """ 

PRINT  USING  "10X ,8(MDD.DD,2X  )"  ;T( 0  ),T( 1 ) , T< 2 ) ,T( 3 ) ,T( 4 ) ,T(5 ) .Tvav.TJav 

PRINT 

PRINT  USING  "6>.,*"Tube     Wall  Temperatures  ( Deg  C>      Tnave      Odp 
The  tat" " " 

PRINT  USING  "6X,""*     12     3     4     5     6   ( De&  C>   CU/m* 
2  .  K  )    (  K  ! 

IF  lpos^e  THEN 
Jj-0 
FOR  Ji-0  TO  5 


"'■,DDD,2X,  14A"  .  J.ToidS 

Tv3  Tldl  Tld2  Tld: 
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3133  Tp< Ji  J»TClpos«6+Jj+6  ! 

3134  .'j-Jj  +  l 
313E  NEXT  Ji 
3137          Tnn=lpos+l 

3136         PRINT  USING  "6X, D, IX ,7( MDD.DD ) , IX . 2! MI .3DE  ) , IX . 1 ( MOD. DO  ) " ; Tnn. Tpt C 
,  Tp  !  1  )  ,Tp(  2  )  ,Tp(  3  ) , Tp<  4  ;,Tp(5;,Twa!  Ipos  ) ,  Odp  (  Ipos  i  , HtuDe(  Ipos  / ,  The  t  ab  ■.  ipos  ) 

3139  60TC  317? 

3140  ENC  IF 

3141  Jj-0 

3143  FOP  li«0  TC  Nht-.l 

3143  FOP  Ji-e  TO  5 

3144  Tpt  Ji  )-T( 1 i»5+Jj+E  ) 
3148  Jj-Jj+1 

3152  NEXT  Ji 

3156  Jj-Ii+J 

3160  FOF  Ji-0  TO  4 

3164  Tn( Ji  >-l+Ji 

3166  NEXT  Ji 

3172  PRINT  USING  "EX ,D, IX , 7(  MOD. DO ) , IX , 2( MZ . 3DE  ) ,  IX , 1 < MOD . DO  ) " . Tn( 1 1  ) , T 

p(0  ),Tp(  1  ),Tp(2  ).Tp(3  ),Tp(4  ),Tp(5  ),Tua< Ii  ),Odp( Ii  ;,Htube( I i  ) , The t ab<  1 1  > 

3176  NEXT  Ii 

3177  PRINT 

3179  PRINT  USIN6  "EX,"  Heat  Flux  and  Tdel  from  Churchi 1 1 /Chu  Correlation  j 
s  ""  ,  KMZ.3DE  ),2X,  I  (MDD.DD  )'.  Och,  The  tab  (  Ipos  ) 

3180  PRINT 
3162  Ok.-l 

3184  IF  Im-0  THEN 

3168  BEEP 

3192  INPUT  "OK  TO  STORE  THIS  DATA  SET  ( 1«Y( default > ,0-N  }?" ,0k 

3196  END  IF 

32001  J«the  counter  for  data  sets 

3204  IF  Ok-1  OR  Im»l  THEN  J-J+l 

3208  IF  Ok-1  AND  In-0  THEN  OUTPUT  SFilel ;Emf ( * >, Volt ( •  ) , Anp( •  ) 

3212  IF  In»l  OR  Ok-1  THEN  OUTPUT  gPl ot , Odp< *  ) , Htube' *  i , Thetab( •  > 

3216  6o_on-l 

3220  IF  In-0  THEN 

3224  BEEF 

3228  INPUT  "WILL  THERE  BE  ANOTHER  RUN  (  1-Y(  default ), 0-N )7" ,Go_on 

3232  Nrun-J 

3236  IF  Go_on-0  THEN  3272 

3240  IF  Bo_on<>0  THEN  Repeat 

3244  ELSE 

3248  IF  J<Nrun+l  THEN  Repeat 

3252  END  IF 

3256  St-1 

3260  BEEC 

3264  INCUT  "ARE  YOU  SURE  YOUF  READY  TO  TERMINATE  ( 1 -Y( DEFAULT i , 0-NO  )' ", St 

3268  Go_on-l 

3272  IF  St>0  THEN  3280 

3276  IF  St-0  THEN  GOTO  3240 

3280  IF  Im-0  THEN 

3284  BEEP 

3288  PRINT 

3292  PFINT  USING  " 10X , " "NOTE :  " "  ,  ZZ,""  data  runs  were  stored  ;r  file  "" 

, 10A" . J-i ,D2fileS 

3296  ASSIGN  gFi iel  TO  • 

3300  0JTPUT  gFile2;Nrun-l 

3304  ASSIGN  ©Filel  TO  DlfileS 

3305  ENTER  @c  i  1  e  I  .  Dst  eS  ,  i_dt  c'  •  ; ,  1 1 1  ,  Bop  ,  Nnt  ,  Nate  .  Nr *  ,  Co^- 
3312  OUTPUT  g-i  le2  ,  DateS  .Ldt  o  ♦  ),  Itt  ,Eop  ,Nht  .Netp.Nrt  ,Ccrr 
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3316 

FOP  1-1  TO  Nrun-1 

3320 

ENTEF  &~  i lei  ;Emf < ♦  >,VgI 1'  •  ;.ftmp<  ♦  ) 

3324 

OUTPUT  SFileZiEnf  (  •  ) , Ocl  t  f  •  i.hno'-  *  ' 

3328 

NEXT  I 

3332 

ASSIGN  PFilel  TO  • 

3335 

PURGE  ■OUMMv" 

35412 

END  IF 

3344 

BEEF 

334e 

PFINT 

3352 

PRINT  USING  -10X , " "NOTE   "",22,""  X-Y  pairs  utre  sto-ed  in  plot  aaia 

lie  - 

",10ft'; J-l ,Pf ileS 

3356 

ftSSIGN  eFile2  TO  • 

3350 

ASSIGN  PPlot  TO  ♦ 

3364 

BEEF 

336B 

SUBEND 

3372  1 

3376  1 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

33801 

3384 

DEF  FNKcu(Tcu  ) 

33881 

OFHC  COPPER 

3392 

Tk-Tcu+273.15     IC  TO  K 

3396 

Kcu-434-.112«Tk      '250-300K  USE  FOR  R-114  S2.2  C 

3400i 

Kcu-433.0-.l«Tk      '200-400K  USE  FOR  R-113  S47.5  C 

3404 

RETURN  Kcu 

3408 

FNEND 

3412' 

3416 

DEF  FNMu(T ) 

3420i 

CURUE  FIT  OF  VISCOSITY 

3424 

Tk-T+273.15     'C  TO  K 

3426 

Mu-EXP(-4.4E36+< 1011 .47/Tk ) )«1 .0E-3   'R-114  170-360  K 

3432' 

Mu-. 00001 34 •  < 10*<503/(  Tk-2.15)>)     'R113 

3436 

RETURN  Mu 

3440 

FNEND 

3444' 

3448 

DEF  FNCp(T  ) 

3452' 

CURUE  FIT  OF  Cp 

3456 

Tk-T+273.15     IC  TC  K 

3460 

Cp-. 40188+1 .ES007E-3«n  +  l .Sl494E-E«Tl  "2-E  .  678S3E- 10«U  "3  'R-114  180-400  K 

3464  1 

Cp-( 925+1 .e3«T )•. 001    'R-113 

3468 

Cp-Cp«1000 

3472 

RETURN  Cp 

3476 

FNENO 

3480' 

3484 

DEF  FNRho(T ) 

3488 

Tk-T+273.15     !C  TO  K 

3492 

X-l-( 1  .8*Tk/753.95 )   !K  TC  R 

3496 

Re- 36. 32+61 . 146414.x' ; 1/3  )+16.4l8015«x  +  l7.<7683e»X" .5+1 .119828«X'2 

3500 

Ro-Ro/ .062428      'R-114 

3504  ' 

Ro=l  .6207479E  +  3-T.{2.218E346  +  T.2.3578291E-3 )     !  R-113 

3506 

RETURN  Ro 

3512 

PNEND 

3516' 

3520 

DEF  FNPr(T)     'GOOD  FOP  R-114/R-11Z 

3524 

Pr-FNCp'T  >«FNMu<  T )/FNK(T ) 

3526 

RETURN  Pr 

3532 

FNEND 

3536' 

3540 

DEF  FNK( T  ) 

3544  ' 

T<360  K  WITH  T  IN  C 

354  = 

K«.071-.00026l»T 
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255Z  RETURN  t- 

3556  FNEND 

3S6P' 

3564  DEF  FNTanht  Fe  I 

3568  P-EXP(Fe) 

3572  O-EXP'-Fe) 

357E  Tanh»( P-0  )/( P  +  0  ) 

3580  RETURN  Tanh 

3584  FNEND 

3588' 

3593  DEF  FNTvsv(U) 

359E  COM  /Cc/  C(7  ) 

3600  T-C<0> 

3604  FOR  1-1  TO  7 

3608  T-T+C(I)«V-I 

3612  NEXT  I 

3616  RETURN  T 

3620  FNEND 

3624  1 

3628  DEF  FNBeta(T ) 

3632  Rop-FNRho(T4. 1 ) 

3636  Ror>-FNRho(  T-.l  ) 

3640  Beta--2/( Rop+Ron  )»(Rop-Ro«  )/  .2 

3644  RETURN  Beta 

3648  FNEND 

3652  DEF  FNPoly(X  ) 

3656  COr  /Cply/  A< 10 , 10  ) , C( 10  ) , B< 4  ) , Nop , Iprnt , Opo , 1 1 og 

3660  Xl-X 

3664  Foly-B<0) 

3668  FOR  1-1  TO  Nop 

3672      IF  Ilog-1  THEN  X1«L0G(X) 

3676      Poly-Poly+BCI  )*X1"I 

3680  NEXT  I 

3684  IF  IIoq-1  THEN  Pol y-EXP( Pol y > 

3688  RETURN  Poly 

3692  FNEND 

3696' 

3700 'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XXX XXXXXXXXX XXX XXXXXX 

3704' 

3708  SUE  Foly<Df ile$< *  >,Np, Itn  ) 

3712  DIM  R( 10 >,S< 10  ),5y( 12  >,Sx<  12  >,Xx< 100  ), Yy( 100  ),Xy( 17  ) 

3716  COM  /Cply/  A( 10 , 10  1 , C( 10  ) , B< 4  ) , N  ,  Iprnt , Opo, I  log 

3720  COM  /Xxyy/  Xp(5),Yp(5) 

3724  FOR  1-0  TO  4 

3728  Bd)«0 

3732  NEXT  I 

3735  In-1 

3740  BEEP 

3744  INPUT  "ENTER  DATA  FILE  NAME  '  ,  Df  i  le$<  0  > 

3748  BEEP 

3752  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS", Np 

3756  BEEP 

3760  INPUT  "LIKE  TO  EXCLUDE  DATA  PAIRS  <  1  -Y  ,  0-n:  DEFAULT  )  >"> ",  led 

3764  IF  Ied-1  THEN 

3768     BEEP 

3772      INPUT  "ENTER  NUMBER  OF  PAIRS  TC  EE  EXCLUDED" , Ipe> 

3776  END  IF 

3780  ASSIGN    6>Fjle    TC    Dfilel'0) 

37e4  N=2 

378S  BEEP 
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379Z 

INPUT  "ENTER  THE  ORDEF  OF 

PC^YNOMIA, 

3796 

FOE 

1-0  TO  N 

3S0C 

Sy( 1  >-0 

3304 

Sx(  I  »-0 

3S0£ 

NEXT  I 

3612 

IE 

led-1  AND  lm-i  THEN 

361E 

FOR  1-1  TO  Ipe> 

3B20 

ENTER  «Fi  le.Xy. •  ! 

3824 

1 

NEXT  I 

3e:e 

END 

IF 

3832 

EOF 

1-1  TO  Np-Ipe> 

3836 

ENTER  ©Fjle. Xy( • ) 

3840 

IF  Opo-0  THEN 

3844 

Y-Xy( Itn-1  » 

3846 

X-Xy(  11  +  Itn  ) 

3eS2 

END  IF 

3e5E 

IF  Opo-1  THEN 

3860 

Y-Xy(5+Itn  > 

3864 

X-Xy(  11  +  Hn) 

3868 

END  IF 

3872 

IF  Opo«2  THEN 

3876 

Y-Xy(5+Itn  ) 

3880 

X-Xy( Itn-1 ) 

3884 

END  IF 

3888 

IF  Ilog-1  THEN 

3892 

X-L06(X  ) 

3e96 

Y-LOG( Y  ) 

3900 

END  IF 

3904 

Xx<  I  )-X 

3908 

Yy(  I  )«Y 

3912 

R(0  )-Y 

3916 

Sy(0  )-Sy(G )  +  Y 

3920 

S(  1  )-X 

3924 

Sx< 1  )-Sx( 1 >+X 

3928 

FOR  J-l  TO  N 

3932 

R( J  )«R< J-l  >*X 

3936 

Sy< J  >-Sy( J  )+R(J  ) 

3940 

NEXT  J 

3944 

FOR  J-2  TO  N»2 

3946 

S« J  )-S< J-l >«x 

3952 

Sx( J  )-Sx( J  )+S< J  ) 

3956 

NEXT  J 

3950 

NEXT  I 

3964 

Sx<0  )-Np 

3966 

FOR 

1-0  TO  N 

3572 

C(  I  )-5y(  I  ) 

3576 

FOR  J-0  TO  N 

3980 

A(I,J  )»Sx(  I  +  J  ) 

3984 

NEXT  J 

3986 

NEXT  I 

3992 

FOR 

1-0  TO  N-l 

3995 

CALL  Divide'  I  > 

4000 

CALL  Subtract* 1  +  1  ! 

4004 

NEXT  I 

4006 

B(N 

>  =  C<  N  )/A(N,N  ) 

4012 

FOP 

1=0  TO  N-l 

4016 

BCN-l-I  i-CCN-1-I  1 

4020 

FOP  J-0  TO  I 

4024 

B( N-l -I  »-B(N-l-I )- 

■h'h-l-l ,N- 

4C2E 

NEXT  J 

' DEFAULT. 
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403Z       E  (  N-  1  -  I  i-B(N-l-I  ).•  •AiN-l-I.N-l-I  ) 

4036  NEXT  I 

404C  IF  IprnW  THEN 

4044      PRINT  USING  "12X,"  "EXPONENT     COEFFICIENT 

4046      FOF  1-0  TO  N 

4052         PRINT  U5INC-  -15X',DC,5>.,MD.7DE-  .  I  ,B<  I  ) 

405E      NEXT  I 

40E0      PRINT  ■'  " 

4064     PRINT  USING  "12X, ""DATA  POINT      >  >      y  (  CALCULATED  !  03 

SCREPANCY"  •'■ 

4066      FOR  1-1  TO  Np 

4072         Yc-Bi0) 

4076         FOR  J-l  TO  N 

4080  Yc«Yc+B< J  >«Xx(  I  )"J 

4084         NEXT  J 

40es         D=Yy; I  )-Yc 

4092         PRINT  USIN6  "15X,3D,4X ,4( MD.5DE, IX  )" ; 1 ,Xx( I  ) , Yy( I  5 , Yc ,0 

4096      NEXT  I 

4100  END  IF 

4104  ASSI6N  CFile  TO  • 

4108  SUBEND 

4112' 

4116  SUE  Divide(M) 

4120  COM  /Cply/  A: 10, 10 ),C( 10 ),B(4 ),N, Iprnt ,0po, Ilog 

4124  FOR  I-M  TO  N 

4128      Ao«A(I,M) 

4132      FOR  J-M  TO  N 

4136  A( I  ,  J  )-A( I , J  )/Ao 

4140      NEXT  J 

4144       C( I  )»C( I  )/Ao 

4148  NEXT  I 

4152  SUBEND 

41561 

4160  SUB  Subtracts:  ) 

41E4  COM  /Cply/  A( 10, 10  ),C( 10  ),B( 4  ),N, Iprnt .Opo, Ilog 

4168  FOR  I-K  TO  N 

4172       FOR  J-K-l  TO  N 

4176  A(  I, J  i-A<K-l,J  >-A! 1,  J  ) 

4180      NEXT  J 

4184      C< I  >-C(K-l  )-C( I  ) 

4188  NEXT  I 

4192  SUBEND 

4196' 

4200  SUB  PI  in 

4204  COM  /Cply/  A(  10,  10  ),C(  10  >,B<:  4  >,N,  Iprnt  ,0po,  Ilog 

4208  COM  /Xyyy/  Xx(5>,Yy!S> 

4212  PRINTER  IS  705 

4216  BEEP 

4220  INPUT  "SELECT  (0«h/h0;:  sane  t  ube  ,  l-h<  Hc  )/h(  sr,  ) "  ,  Irt 

4224  BEEF 

4226  INPUT  "WHICH  Tsat  ( 1=6 . 7 , 0  =  -2 . 2  ) " , 1 se t 

4232  Xr.in-0 

4236  Xr.ax'10 

4240  Xstep=2 

4244  IR  Irt-0  THEN 

4245  Yr,in  =  C 
4252     Ymax-1.4 
4256     >step=.2 
426C     ELSE 
42&4      Yr.in«0 
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4:ge 

Ymax-15 

iZ"2 

Ystep-5 

4276 

END  IF 

4280 

BEEF 

4284 

PRINT  "IN;5pl.IF  2300,2200,6500,6600." 

4286 

PRINT  "SC  0, 100,0, 100. TL  2,0." 

4292 

Sf >«100/( Xmax-Xmin  ) 

4296 

Sfy-100''  Ymax-Yr-ur,  1 

4  300 

PRINT  •'PL'  0,0  PD" 

4304 

FOR  Xe-Xrun  TO  Xma>  STEP  Xstep 

4308 

X-( Xe-Xnun  )»Sfx 

4312 

PRINT  "Pn";X,",0;  XT." 

4316 

NEXT  Xs 

4320 

PRINT  "PA  100,0;PU.  " 

4324 

PRINT  "PU  PA  0,0  PO" 

4328 

FOR  Ya-Yrun  TO  Yma.  STEP  Ystep 

4332 

Y-( Ya-Yrun  )»Sfy 

4336 

PRINT  'PA  0, " . Y, *YT' 

4340 

NEXT  Ya 

4344 

PRINT  "PA  0,100  TL  0  2" 

4348 

FOR  Xa-Xmin  TO  Xmax  STEP  Xstep 

4352 

X»(Xa-Xnjn  >»Sfx 

4356 

PRINT  "PA" ;X,",100.  XT" 

4360 

NEXT  Xa 

4364 

PRINT  "PA  100,100  PU  PA  100,0  PD" 

4368 

FOR  Ya-Ynin  TO  Ymax  STEP  Ystep 

4372 

Y-( Ya-Yrun  )»Sf y 

4376 

PRINT  "PD  PA  100, " ,Y, "YT" 

4380 

NEXT  Ya 

4384 

PRINT  "PA  100,100  PIT 

4388 

PRINT  "PA  0,-2  SR  1.5,2" 

4392 

FOR  Xa=Xrun  TO  Xmax  STEP  Xstep 

4396 

X=»(  Xa-Xnm  >»Sfx 

4400 

PRINT  "PA"  ;X,  "  ,0.  " 

4404 

PRINT  "CP  -2,-1 .LB" ;Xa; "" 

4408 

NEXT  Xa 

4412 

PRINT  "PU  PA  0,0" 

4416 

FOR  Ye«Yr,:n  TO  Ymax  STEP  Ystep 

4420 

IF  AES: Ya  X 1 .E-5  THEN  Ye«0 

4424 

Y»{Ya-Ycun )*Sfy 

4428 

PRINT  "PA  0,  " . Y, " " 

4432 

PRINT  "CP  -4, -.25. LB" . Ya . " " 

4436 

NEXT  Ya 

4440 

Xlabels-'Oil  Percent" 

4444 

IF  Irt-0  THEN 

4446 

Ylabel8-"h/h0Z" 

4452 

ELSE 

4456 

YlabelS="h/hsmooth" 

4460 

END  IF 

4464 

PRINT  "SF  1.5.2.PU  PA  50,-10  CP  "  . -LEN( X iabelS >/ 

4468 

PRINT  "PA  -11,50  CF  0, " ,-LEN( YlabelS  )/2»5/6. "Dl 

4472 

PPINT  "CP  0,0" 

4475 

Ipn-B 

448e 

EEEC 

4  484 

INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  ',  1=\,P  =  N)-' 

448E 

Icn-0 

44S2 

IF  01 p=l  THEN 

4  4  9  = 

BEEP 

4500 

INPUT  "ENTEF  THE  NAME  OF  THE  DATA  FILE " , D_f 1 leS 

4504 

EEEP 

"CLE"  .X  labels. 
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4  5Cc 

INPUT  "SELECT  <0=LINEAF,  1-L0G<  X  ,'>:",  I  log 

45:2 

ASSIGN  PFile  TC  C_f 1 leS 

451E 

BEEF 

452C 

INPUT  "ENTEF  THE  EEC-INNING  RUN  NUHEER'.Md 

45:-: 

EEEF 

4526 

INPUT  "ENTER  THE  N'JMBEF  OF  X-Y  PAIRS  STORED" 

1  Noa  it 

4532 

BEEP 

453£ 

INPU"  "ENTER  DESIRED  HEAT  FLUX",& 

4540 

BEEP 

4544 

PRINTER  IS  1 

4546 

PRINT  USIN6  "4X,  ""Select  a  symbol  

455: 

PRINT  USING  "4X,-"1   Star    2   Pluj  sign""" 

455S 

PRINT  USING  "4X,""3   Circle  4   Square-"" 

45Be 

PRINT  USING  "4X.--5   Rombus""" 

4554 

PRINT  USING  "4X.--6   Right -s i de-uc  triangle-1 

4566 

PRINT  USING  "4X,""7   Up-£ 1 de-down  triangle""' 

4572 

INPUT  Syr, 

4576 

PRINTER  IS  705 

4580 

PRINT  "PU  DI" 

4584 

IF  Syrt-1  THEN  PRINT  "SM»" 

4588 

IF  Syn-2  THEN  PRINT  -Sn+" 

4592 

IF  Syw-3  THEN  PRINT  "SMo- 

4596 

Nn-4 

4600 

IF  Ilog-1  THEN  Nn-1 

4604 

IF  Md>!  THEN 

4608 

FOR  1-1  TO  (Md-1  ) 

4612 

ENTER  3File.Xa.Ya 

4616 

NEXT  I 

4620 

END  IF 

4624 

01-0 

4528 

IF  Ilog-1  THEN  Q«L06(Q> 

4632 

FOP  1-1  TO  Npairs 

4636 

ENTER  @File.Xa,B(»  ) 

4640 

Ya-E(0  ) 

4644 

FOR  K-l  TO  Nn 

4648 

Ya-Ya+B(K  )»0*K 

4652 

NEXT  K 

4656 

IF  Ilog-1  THEN  Ya-EXP( Ya  ) 

4660 

IF  Ilog-0  THEN  Ya-Ql/Ye 

4664 

IF  Irt-0  THEN 

4668 

IF  Xa-0  THEN 

4672 

Yo-Ya 

4676 

Ya-1 

4580 

ELSE 

4684 

Ya-Ya/Yo 

468c 

END  IF 

4ES2 

ELSE 

4696 

Hsm-FNHsmootM C,Xa, Isat  ) 

4700 

Ya-Ya/Hsm 

4704 

END  Ic 

4708 

Xx(  1-1  >-Xa 

4712 

Yy(  1-1  )-Yfl 

4716 

X  =  (  Xa-Xr,in  )»Sf  >. 

4720 

Y-'  Ya-Ymin  >«Sfy 

4724 

IF  Syn>3  THEN  PRINT  "SM" 

4~2c 

IF  Synv4  THEN  PRINT  " SF  1.4,2.4" 

47Z2 

PRINT  "PA",X,Y,-" 

4~3E 

IF  Syr.  3  THEN  PRINT  "SR  1.2.1.6" 

4740 

IP  Syn-4  THEN  PRINT  " UC2 . 4 . 9= , 0 , -8 , -4  ,  0  .  G . 6  ,  4 

,0,  .  ' 

4744 

IF  Sym=5  THEN  PRINT  "UC3 , e , 99 , -3 ,-E , -3 , E ,3 , E, 

3 .  -£  ■ 
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IF  Sym-E 

THEN 

PRINT 

"UCC, 

5.; 

i.99,3, 

-6,- 

-6 

,e. 

3,1 

3;  " 

IF  5y«--" 

THEN 

PRINT 

"UCO, 

_c 

.2  99  - 

3,6, 

,6 

,e. 

-z 

.-£ 

next  I 

BEEP 

ASSI6N  PF 

lie  TC  • 

END  IF 

PRINT  "PL 

i  sr" 

BEEP 

INPUT  "UANT  TO 

PLOT 

A  POl 

.YNOMIAL  < 

l->  , 

0. 

•N  ) 

1  " 

,01 

IF  CHd-1 

THEN 

BEEP 

INPUT 

■SELE 

CT  (0- 

■LINEAR,: 

i-LOG(X 

,  Y  > 

I" , 

,n 

og 

ipr-nt" 

1 

4""4£ 

4753 

4756 

4760 

4764 

4766 

4772 

4776 

4780 

4764 

4788 

4792 

4796 

4800      CALL  Foly<  Itn  ) 

4804     FOP  Xe-Xmin  TO  Xma>  5TEC  Xstep/25 

4808  lcn-lcn+1 

4612  Ya  =  FNc'oly(  Xe  ) 

4ei6  V-t Ye-Ywin)«Sfy 

4620  X»(  Xa-Xrun  )*Sf  x 

4824  IF    Y<0    THEN    Y-0 

4828  IP    Y)100    THEN    GOTO    4B68 

4632  Pu-0 

4836         IF  Ipn-1  THEN  Idf-Icn  MOD  2 

4640  IF  Ipn-2  THEN  Idf-Icn  MOD  4 

4844         IF  Ipn-3  THEN  Idf-Icn  MOD  B 

4846         IF  Ipn-4  THEN  Idf-Icn  MOD  16 

4852         IF  Ipn-5  THEN  Idf-Icn  MOD  32 

4856         IF  Idf-1  THEN  Fu-1 

4860  IF  Pu-0  THEN  PRINT  "PA" , X , Y , *PD ' 

4864  IF  Pu-1  THEN  PRINT  "PA" , X , Y , "PU- 

4666     NEXT  Xa 

4872     PRINT  -PU' 

4876      Ipn-Ipn+1 

4880      60T0  4480 

4884   END  IF 

4888   BEEP 

4892   INPUT  "WANT  TO  QUIT  ' 1  =  Y , 0-N >' "  ,  Iqui  t 

4895   IF  Ipuit-i  THEN  4904 

4900   GOTO  4480 

4904   PRINT  "PU  SP0" 

4906   SUBEND 

4912   SUB  Stats 

4916   PRINTER  IS  701 

4920   J-0 

4924   K-G 

4928   BEEP 

4932   INPUT  "PLOT  FILE  TO  ANALYZE? ', F l leS 

4936   ASSIGN  §File  TO  FjleS 

4940   BEEP 

4944       INPUT     "LAST    RUN    No?' 0=OUI T ) " , Nn 

4948      IF    Nn=e    THEN   5092 

4952      Nn-Nn-J 

4956      Sx-0 

4960        Sy-0 

4954  5:-0 

496e  Sxj-e 

4972  5ys=0 

4976  S:s-0 

4980  FOP    1=1    TO   Nn 

4964  J-J+l 
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4985  ENTEF  @"ile.C,T 

4992  H'Q  "I 

499E  5>-Sx+Q 

5000  Sxs-Sxs+O*: 

5004  Sy-Sy+T 

5008  Syi-Sys+T~2 

5012  S:-S:+H 

5016  Szs-Srs+H": 

5020  NEXT  I 

5024  Qave-S//Nn 

5028  Tave-Sy/Nn 

5032  Have-S:/Nn 

5035  Sdevq-SQP( ABS< ( Nn*SxS-S» *2  l/(Nn*< Nn-1  )  )  >  i 

5040  Sdevt-SQR(ABS((Nn«Sys-S>~2  )/<Nn»(Nn-l  )  > ) ) 

504  4  Sdevh-SQR< A8S< (Nn«S;s-Sr"2  >/(Nn«(Nn-l  )  )  )  ) 

5048  Sh»l00«5devn/Have 

5052  Sq-100«Sdevq/Qave 

5055  St-100«Sdevt/Tave 

5060  IF  K-l  THEN  5084 

50B4  PRINT 

5068  PRINT  USING  "11X,""DATA  FILE ■ " " , 14A" . F i leS 

5072  PRINT 

5076  PRINT  USING  "11X,""RUN  Htube     SdevH   Qdp        SdevQ   Thetab  SdevT" 

5080  K-l 

5084  PRINT  USING  •UX,DD,2<2X,0.3DE,1X,3D.2D),2X,DD.3D,1X,3D.2D-; J .Have , Sh , Qave 

, Sq , Tave , St 

5088  GOTO  4940 

5092  ASSIGN  gFilel  TO  • 

5096  PRINTER  IS  1 

5100  SUBEND 

5104  SUB  Coef 

5108  COM  /Cply/  A<  10,10 >,C<  10  !,B<4  ),N, Iprnt.ODO, Iloc 

5112  BEEP 

5116  INPUT  "GIUE  A  NAME  FOR  CROSS-PLOT  FILE",Cpfl 

5120  BEEP 

5124  INPUT  "OUTPUT  TYPE  (0-q  vs  Dt ,  1-h  vs  Dt ,  2-h  vj  q)",0pc 

5126  CREATE  BOAT  CpfS.6 

5132  ASSIGN  ©File  TO  CpfS 

5136  BEEP 

5140  INPUT  "SELECT  ( 0-L INEAP , 1 «L06( X , Y  )  )" , I  log 

5144  BEEP 

5148  INPUT  "ENTER  OIL  PERCENT  (- 1 -STOP >", Bop 

5152  BEEP 

5156  INPUT  "ENTER  TUBE  NUMBER  (1,  2,  3,  4,  OR  5)",Itn 

5160  IF  Bop'.0  THEN  5176 

5164  CALL  Polyv Itn) 

5168  OUTPUT    PFjle.6op.Bc  ) 

5172  GOTO  5144 

5176  ASSIGN  BFile  TO  * 

5180  SUBEND 

51841 

5188'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

5192' 

5196  SUE  Plot 

5200  COM  /Cply/  A'.ie,lE>,C(  10),E(4),Nop,Iprnt,Opc,Ilog 

5204  DIM  Xy! 17) 

5208  INTEGEF  Ii 

5212  PRINTER  IE  1 

5216  BEEP 
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s::e 
sz:< 

5226 
5:3: 

5236 

5r*e 

5244 

5248 
5252 
S2SE 
5250 
5264 
52Be 
5272 
5276 
5280 
5284 
5288 
5292 
5295 
5300 
5304 
5306 
5312 
5316 
5320 
5324 
5328 
5332 
5336 
5340 

5344 

5348 
5352 
5356 
5360 
5354 
5368 
5372 
5376 
5230 
5384 
538e 
5392 
5396 

5400 
5404 

5406 
5412 
5416 
5420 
5424 
5428 
5432 
5436 

5440 

5444 
5446 
5452 
5456 


"4X  , 

"BX, 

■ex, 

"EX, 


" Se  lee  t  Opt  1  or, 


Idv-1 

INPUT  "LIKE  DEFAULT  VALUES  FOF  PLOT  I  l»=Yl  DEFAULT  1 ,0-N 
Opo-0 
BEEP 

PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
INPUT  Opo 
BEEP 

INPUT  "SELECT  UNITS  ( 0=51 ( DEFAULT  ), 1-ENGLlSH  ) ' 
PRINTER  IS  705 
IF  IdvOl  THEN 
BEEP 

ENTER  NUMBER  OF  CYCLES  FOR 


■0 

■  1 

-  -> 


versus 
versus 
versus 


delte-T( DEFAUL'  >' 
delte-T"  "  " 
Q""  " 


,  Iun 


X-AXIS" ,Ca 
■ENTER  NUMBER  OF  CYCLES  FOR  Y-AXIS", Cy 
'ENTER  MIN  X-VALUE  (MULTIPLE  OF  10)-,Xrun 


INPUT 
BEEP 
INPUT 
BEEP 
INPUT 
BEEP 

INPUT  'ENTER  MIN  Y-VALUE  (MULTIPLE  OF  10)",Ynin 
ELSE 

IF  Opo-0  THEN 
Cy-2 
Cx-2 
Xrun- 1 
Ymm-1000 
END  IF 

IF  Opo-1  THEN 
Cy-2 
Cx-2 
Xnin"l 
Ymin-100 
END  IF 

IF  0po=2  THEN 
Cy-2 
Cx«Z 

Xmm=1000 
Ymin=100 
END  IF 
ENC  IF 
BEEF 

PRINT  "IN.SP1;IP  2300,2200,6300,6800,- 
PRINT  "SC  0,100,0, 100, TL  2,0: " 
Sf>-100/Cx 
Sfy-100/Cy 
EEEP 

INPUT  "UANT  TO  BY-PASS  CAGE  (1-Y,  0-NO( DEFAULT  )", Ibp 
IF  Ibp=l  THEN  5906 
PRINT  "PU  0,0  PD" 
Nn=9 
FOF  1=1  TO  Cx+1 

Xat-Xmin»10'  ! 1-1  ) 
IF  I-Cx+I  THEN  Nn-1 
COR  J-I  TC  Nn 

IF  J-l  THEN  PRINT  "TL 

IF  J-2  THEN  PRINT  "TL 

Xa-Xat«J 

X»LGT( Xa/Xwin  )*Sf> 

PRINT  "PA",>, -,  e,  XT.' 


2  C" 

1  0" 
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S46C 

ne>:t  J 

5454 

NEx~  3 

S46E 

PRINT  "PA  10C.  CPU.  " 

54?: 

PRIN*  "PU  PA  C,0  PD" 

547G 

Nn-9 

5480 

C0F  1-1  TO  Cy+1 

5484 

Yat-Ymm*  l0"t  1-1  ) 

54S& 

IF  1-Cv+l  THEN  Nn-1 

5492 

FOF  J-l  TO  Nn 

549B 

IF  J-l  THEN  PRINT  ' 

'TL 

:  e" 

5500 

IF  J-2  THEN  PRINT  ' 

TL 

1  0" 

5504 

Y*-Yat  «J 

5508 

Y-LGT< Ya/Ynin  >»Sfy 

5512 

PR3NT  'PA  0, '  . V  ,  "Yl 

5516 

NEXT  J 

5520 

NEXT  3 

5524 

PR3NT  "PA  0,100  TL  0  2" 

5526 

Nn«9 

5532 

FOR  1-1  TO  Cx  +  1 

5536 

Xat-X«un«10"< 1-1  ) 

5540 

IF  I-Cx+1  THEN  Nn-1 

5544 

FOR  J-l  TO  Nn 

5548 

IF  J-l  THEN  PRINT  ' 

TL 

0  2" 

5552 

IF  J>1  THEN  PRINT  " 

'TL 

0  1" 

5556 

Xa-Xat»J 

5550 

X-L6T(Xa/Xmin  >*Sfx 

5564 

PRINT  -PA" ; X, " , 100, 

XT" 

5568 

NEXT  J 

5572 

NEXT  I 

5576 

PRINT  "PA  100, 100  PU  PA  10C 

),0 

PD" 

5580 

Nn-9 

5584 

FOR  1-1  TO  Cy+1 

5588 

Yat-Ymin»10-<  1-1  ) 

5592 

IF  I-Cy+1  THEN  Nn-1 

5596 

FOP  J-l  TO  Nn 

5500 

IF  J-l  THEN  PRINT  ' 

TL 

0  2' 

5504 

IF  J, 1  THEN  PRINT  ' 

'TL 

0  1" 

560E 

Ya-Yat«J 

5612 

Y«LGT<  Ya/Ymin  >*Sfy 

5616 

PRINT  "PD  PA  100," , 

V,  ' 

•YT" 

5620 

NEXT  J 

5624 

NEXT  I 

5528 

PRINT  'PA  100,100  PU" 

5632 

PPINT  "PA  0,-2  SR  1.5,2" 

5636 

I i«LGT( Xmin ) 

5640 

FOF  3-1  TC  Cx+1 

5644 

Xa->m:n»10-<  1-1  ) 

5648 

X-LGT( Xa/Xmin  >»Sfx 

5652 

PRINT  "PA"  .  X,  " ,0.  " 

5656 

IF  Ii>-0  THEN  PRINT  "CF 

>  _■ 

:,-2.leic.pf  -: 

,2; LB"; 

5B50 

IF  1 1 <e  THEN  PRINT  "CP 

-2 

,-2;LB10iPR  e,2 

.LB"  .  li 

5664 

Ii-Ii+1 

5SGE 

NEXT  I 

5E72 

PPINT  "PU  PA  0,0" 

5E7E 

1  i=_C-T<  Ymin  ) 

5E8e 

Y10-10 

5564 

FO11  1-1  TC  Cy+1 

SEBe 

Ye-Ymin*10"< 3-1  ) 

56S2       Y«LGTt Ya /Ymin >*Sf y 
5E9E       PPINT  "PA  0,"  ,  Y,  "•' 


171 


S70C  print  -cp  -i,-.:s.LEie.PF  -:,:.le" . i i . " ■ 

5"?04  Il-Il  +  l 

5~0S  NEXT    I 

b?i:  beep 

5-1E  Idl-l 

572e  INPUT  "WANT  USE  DEFAULT  LABELS  <  1-Y  (  DEFAULT  )  .  0-N  )"> "  ,  I  dl 

5724  IF  Idl :s 1  THEN 

5726  BEEP 

5^32  IN^UT  "ENTEF  X-LABE^",XiebelS 

S~3E  BEEP 

5740  INPUT  "ENTER  v-LABEL' ,>label$ 

5744  END  IF 

574E  IF  Opo-:2  THEN 

5752  PRINT  "SP  1,2; PU  PA  40,-14." 

5756  PRINT  -LB(T,PF  -1.6,3  PD  PR  1.2,0  PU.PF  . 5 , -4 , LBuo  .  PR  .5,1;" 

5760  PRINT  "LE-T;PR  .5 , -1 ; LBsat ; PR  .5,1." 

5764  IF  Iun-B  THEN 

5768  PRINT  "LB  )  /  ( K > " 

5772  ELSE 

5776  PRINT  "LB  )  /  (F  >" 

5780  END  IF 

5794  END  IF 

5788  IF  Opo-2  THEN 

5732  IF  Iun-0  THEN 

5796  PRINT  "SR  1.5,2;PU  PA  40,-14;LBq  /  <U/n;SR  1,1. 5. PR  0 . 5 , 1 ; LB2 . SR  1 

.5,2;PR  0.5,-1 .LB)" 

5800  ELSE 

5804  PRINT  "SR  l.E,2.PU  PA  34,-14.LBq  /  <Btu/hr;PR  .5..5.LE..PP  .5, -.5. 

5808  PRINT  "LBftiPR  .5,1;SR  1.1.5.LB2.SR  1.5, 2. PR  .5,-l.LB)." 

5ei2  END  IF 

5816  END  IF 

5820  IF  Opo-0  THEN 

5824  IF  Iun-0  THEN 

5828  PRINT  "SR  1.5.2.PU  PA  -12.40.DI  0,1. LBq  /  (U/m.PR  -1.0.S.SR  1.1.5.L 

62. SR  1.5, 2. PR  1 , .S:LB>" 

5S32  ElSE 

5S3E  PRINT  "SF  1.5, 2. PL  PA  -12,32.01  0,1. LBq  /  (Btu/hr.PR  - . 5 , . 5 . LB . . PR 

.5. .5." 

5640  PRINT  "LBftiSR  1,1. 5. PR  -1..5.LB2.PF  1..5.SR  1.5, 2. LB)" 

5844  END  IF 

5848  END  IF 

5eS2  IF  Opo>0  THEN 

585E  IF  lun>0  THEN 

5860  PRINT  "SR  1.5.2.PU  PA  -12,38;DI  0,1. LBh  /  (U'n.PP  -1..5.SR  1.1.5.LE 
2.SP  1 .5,2. PF  .5. .5.  " 

5864  PRINT  "LB. .PR  .5,0, LBK  )" 

5858  ELSE 

5872  PRINT  "SR  1.5,2:PU  PA  -12.28.DI  0,1. LBh  /  (Btu/hr.PR  - . 5  ,  . 5 . LE .  . PR 

.  5  ,  .  5  .  " 

5676  PRINT  "LBft.PP  -  1  .  . 5 . SF  1,1.5;LB2;SR  1.5,2;PP  .5,.5.LE..PR  .5, .5. 
LBF)" 

seee  end  if 

53E4  END  IF 

SSSt  IF  IdI-0  THEN 

5e92  PRINT  "SF  1.5,2;PU  FA  5C.-1E  CP" ; -LEN( XlabelS  )/2 ; "0; LE" ; XlebelS . " " 

599E  PRINT  "PA  -14,50  CF  0, " ; -LEN(  YlabelS  ).'2*5/E  ;  "Dl  0, 1  ;LB" ;  YlebelS; " " 

5S00  PRINT  "CF  0.0  Dl" 

5904  END  IF 

590E  Ipn-0 
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ssi: 

S91E 

ss:e 
sszc 
59:e 
59;: 

5936 
5940 
5944 
5948 
5952 
5956 
5960 
E954 
5958 
5972 
5975 
5980 
5984 
5988 
5992 
5995 
6000 
6004 
6008 
6012 
6016 
6020 
6024 
6028 
6032 
6036 
6040 
6044 
6048 
6052 
6056 
6060 
6064 
60SS 
6072 
6077 
6080 
6084 
6086 
6052 
6096 
6100 
6104 
6108 
6112 
6116 
6120 
6124 
6126 
6132 
E136 
6140 

E14C 

6146 


synbol 


Plus  sign' 
Square" " " 


Repeat  ' 
X  1 1-1  .E  +  E 
Xul  —  1  .E  +  E 
lcn-0 
BEEP 
CH-1 

INPUT  "WANT  TO  PLOT  OATA  FROF  A  FILE  (  1  -Y  1  DEFAULT  >,  0«N  >?",  01 
IP  Ot-  »1  THEN 
BEEF 

INPUT  "ENTEF  THE  NAME  OF  THE  DATA  F ILE " , Df 1 le$i 0 > 
ASSIGN  BFile  TO  Df ileS(0  ) 
BEEP 

Npairs*20 

INPUT  "ENTEF  THE  NUMBEF  OF  X-Y  PAIRS  STORED! DEFAULT  =  20  )", Npa irs 
BEEP 

Itn-Itn+1 

INPUT  "ENTEF  TUBE  NUMBER  (1,  2,  3,  4,  OF  5>",Itn 
BEEP 

PRINTER  IS  1 

INPUT  "WANT  DEFAULT  SYMBOLS'?  (YES-0  (  DEFAULT  )  ,NO«l  )",  Symb 
Syn-Itn+2 
IF  Synb-0  THEN 

60T0  6036 
END  IF 

PRINT  USING  "4X 
PRINT  USING  "6X 
PRINT  USING  "GX 
PRINT  USING  "BX,"'5 
PRINT  USIN6  "6X,"6 
PRINT  USING  "BX,"*7 
INPUT  Syn 
PRINTER  IS  705 
PRINT  -PU  DI" 
IF  Sym«=l  THEN  PRINT 
IF  Sym>2  THEN  PRINT 
IF  Sym-3  THEN  PRINT 
FOF  1-1  TC  Npairs 
ENTEF  @~iie,Xy< «  ) 
IF  Opo-0  THEN 

Ya  =  Xy( Itn-1  ) 

Xe«Xy< 1 1  +  Itn  ) 
END  IF 
IF  Opo-1  THEN 

Ya-Xy(5+Itn ) 

Xe-Xy< 1 1+1 tn ) 
END  IF 
IF  Opo-2  THEN 

Ya=Xy(5+Itn ) 

Xa«Xy< Itn-1 ) 
END  IF 

IF  Xe'Xll  THEN  X11=X* 
IF  Xa:-Xul  THEN  Xul=Xe 
IF  Iun=l  THEN 

IF  Opo.2  THEN  Xa«Xe»l  .E 

IF    Opc0    THEN    \e=te>  .1761 

IF    Op o=C-   THEN    Ye=Ye*.317 

IF    Opo=2    THEN    Xe-Xa«  .31" 
END    IF 

>=L6T; Xa  'Xmin  l»Sf> 
V-^6T' Ya/Ymin >»£fy 


Select    e 
1      Star 
3      Circle 

5  Rombus " " ' 

6  Right-side-up    triangle" 

7  Up-s  1  ae-dowr,    triangle""' 


'SM«  " 
'SM+" 
'SMo" 
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Eis:  Kj«e 

6  155  CALL  5ynb(X,>  ,5yr,,Icl,Kj  ) 

Ei6c  goto  e:i: 

E164  IF  Syr.>3  THEN  FR1MT  "SM" 

5158  IF  Syrr-4  THEN  PRINT  "SR  1.4,2.4' 

B170  IF  Icl-0  THEN 

E17E  PRINT  "PA",X,Y, "" 

Else  else 

E184  PRINT  "PA" ,X.Y, "PD" 

6188  END  IF 

61S0  IF  Sym/3  THEN  PRINT  "SR  1.2,1.6" 

519E  IF  Syn-4  THEN  PRINT  "UC2 , 4 , 99  ,0 ,-6 , -4 , 0 , e . E , 4 , 0 .  " 

E200  IF  Syn-5  THEN  PRINT  "UC3 , 0 , 99 , -3 , -6 , -3 , E , 3 , E , 3  ,  -6 . " 

E204  IF  Sym-E  THEN  PRINT  "UC0 ,5 . 3 , 99,3 , -8 , -E ,0, 3 , 8 ; " 

E208  IF  Syr,*?  THEN  PRINT  "UC0  ,  -5  .  3  ,  99 .  -3  ,  8  ,  E  ,  0  ,  -3  ,  -8  .  " 

6210  NEXT  I 

6216  PRINT  "PIT 

6220  BEEP 

6224  Ilab-1 

6228  INPUT  "WANT  TO  LABEL"?  (  1-Y  (  DEFAULT  )  ,0-N  )",  I  lab 

6232  IF  Ilab-1  THEN 

6236  PRINT  "SP0.SP2" 

6240  BEEP 

6244  IF  Klab-0  THEN 

6248  Xlab-65 

6252  Ylab=eS 

6256  INPUT  "ENTER  INITIAL  X,Y  LOCATIONS- , X  lab , Ylab 

6260  Xtt-Xlab-5 

6264  Ytt-Ylab+8 

6268  PRINT  "SR  1,1.5" 

6272  PRINT  "SM;PA*,Xtt,Ytt,"LB     Tube   X        File" 

6276  Ytt-Ytt-3 

6280  PRINT  -PA-.Xtt ,Ytt , -LE      No    Oil   Name" 

5284  IF  Syn-=1  THEN  PRINT  "Sn«" 

6288  IF  Syn-2  THEN  PRINT  "  SM+ " 

6292  IF  Syn-3  THEN  PRINT  "Sflo" 

6296  Klab-1 

6300  END  IF 

6304  Kj-1 

6308  CALL  5yi"ib<Xlab,Ylab,Sym,Icl,Kj  ) 

6312  PRINT  "SR  l.l.SiSM" 

6316  IF  Syn<4  THEN  PRINT  "PR  2,0" 

6320  PRINT  "PR  2,0. LB''  .  Itn.  "" 

6324  BEEP 

E32E  INPUT  "ENTER  B0P< 0-DEPAULT >" , Bop 

E322  IF  Bop  10  THEN  PRINT  "PR  3  ,  0 ; LB " . Bop . " " 

6336  IF  Bop>9  THEN  PRINT  "PP  1 .5,0;LB" :Bop ; " " 

6340  PRINT  "PR  2.0. LB" ;Df ile$(0  ); "" 

6344  PRINT  "SD0.SD1.SR  1.5,2" 

634E  Ylab=Ylab-E 

6252  END  IF 

6355  EEEP 

£360  ASSIGN  @File  TO  • 

E354  Xll-Xll/1.2 

£365  Xul-Xul»1.2 

6Z~:  i  GOTO  6040 

E3^E  END  IF 

£380  PRINT  "PU  SM" 

£364  BEEP 

E36E  Gc  on=l 
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E3S:  INPUT  "UANT  TO  FLOT  h    POLYNOMIAL  •  1  ■  '  '  DEFAULT  i  ,  P  =  N  :■  °  '  .  &c  _  on 

EZ9E  IF  Gc_on-i  THEN 

6400  EEE~F 

6404  PRINTEF  IS  I 

640E  INPUT  "WANT  DEFAULT  LINE  TYPE9  ^  YES-C  (  DEFAULT  ),  NO-  1  ;  "  ,  Lr, 

6412  Ipn-Itn 

641E  IF  Ln-0  THEN 

6420  60TC-  6448 

6424  END  IF 

642E  PRINT  USING  "4X, ""Select  line  type:""* 

6432  PRINT  USING  *6X,""0     Solid  line""" 

6436  PRINT  USING  "BX,""1      Dashed"" 

6440  PRINT  USING  "6X,""2,,,5  Longer  line  -  dash 

6444  INPUT  Ipn 

6446  PRINTER  IS  705 

6452  BEEP 

645E  I log-1 

6460  INPUT  "SELECT  ( 0-LIN , 1-LOG! DEFAULT  ))", I  log 

6464  Iprnt-1 

6466  CALL  PolylDf ile$< •  ),Npairs, Itn) 

6472  FOR  X*-0  TO  Cx  STEP  Cx/200 

6476  Xa-Xmin«10*Xx 

6480  IF  Xa<Xll  OR  Xe>Xul  THEN  6572 

6484  Icn-Icn+1 

6488  Pu-0 

6492  IF  Ipn-1  THEN  Idf-Icn  MOD  2 

6496  IF  Ipn-2  THEN  Idf-Icn  MOD  4 

6500  IF  Ipn-3  THEN  Idf-Icn  MOD  8 

6504  IF  Ipn-4  THEN  Idf-Icn  MOD  16 

6508  IF  Ipn-5  THEN  Idf-Icn  MOD  28 

6512  IF  Idf-1  THEN  Pu-1 

6516  Ya-FNPoly(Xa > 

6520  IF  Ya<Yrun  THEN  6572 

6524  IF  Iun-1  THEN 

6528  IF  0po<2  THEN  Xe-XaM.8 

6532  IF  Opo>0  THEN  Ya-Ya«.1761 

6536  IF  Opo-0  THEN  Ya-Ya«.3l7 

6540  IF  Gpo=2  THEN  Xa-Xa«.317 

6544  END  IF 

6548  Y-L6T( Ya/Ymin  )«Sfy 

6552  X«L6T(Xe/Xpiin  )«Sf  x 

6556  IF  Y<0  THEN  Y-0 

6560  IF  Y>100  THEN  GOTO  6572 

6564  IF  Pu-0  THEN  PRINT  "PA* , X , Y, "PD" 

6566  IF  Pu«l  THEN  PRINT  "PA"  ,X , Y , "PU" 

6572  NEXT  X> 

6576  PRINT  "PU" 

6580  END  IF 

6584  eEEP 

E588  INPUT  "WANT  TO  QUIT  ( 1 -Y , P  =  N( DEFAULT  )  ) " , Iqt 

£592  JF  Iqt-1  THEN  6600 

E59E  GOTO  5916 

6600  PRINT  "PU  PA  0,0  SP0" 

6504  SUBEND 

6E0EI 

6612IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

66161 

6620  SUE    Syeib<X,Y,Sym,Icl  ,Kj > 

6624  IF  Syr.,'3  THEN  PRINT  "  SM  " 

6626  IF  Syn<4  THEN  PRINT  "SF  1.4,2.4" 
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6E32 

Vad-0 

6E36 

IF  Kj»l  THEN  ^ad= . B 

6E40 

IF  Ic 1-0  THEN 

6644 

PRINT  "PA" ,X, Y+Yad, "" 

6646 

ELSE 

665Z 

PRINT  "PA" ,X ,Y+Yad, "PD" 

6656 

END  IF 

6660 

IF  Sym>3  THEN  PRINT  "SR  1.2,1.6" 

6664 

IF  5yn-4  THEr.  PRINT  "UC2  ,  4  .  SS  ,  0  ,  -E  .  -4  ,  0  ,  G  ,  6  ,  4  ,  0  ;  " 

6666 

IF  Syr>«5  THEN  PRINT  "UCZ  ,  0  ,  95  ,  -3  ,  -6  ,  -3  ,  E  ,  3  ,  6  ,  3  ,  -6  .  " 

6672 

IF  Syn-6  THEN  PRINT  "UC0 , 5 . 3 , 95 , 3 , -E , -6 , 0 , 3  ,  8  .  " 

6676 

IF  Sym-7  THEN  PRINT  "UC0 , -5 . 3 . 99 , -3 , 8 , 6 . 0 , -3 , -8 .  " 

6660 

IF  Kj-1  THEN  PRINT  "SM.PR  0,-.8" 

6684 

SUBEND 

66B8' 

6692' 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

66961 

6700 

SUE  Fixup 

6704  1 

FILE   FIXUP 

67081 

6712 

DIM  Enf (34  ),Amp( ll),Volt(4),Ldtc<4) 

E716 

BEEP 

6720 

INPUT  "OLD  FILE  TO  F  IXUP" , D2f i leS 

6724 

ASSIGN  grile2  TO  D2fileS 

6728 

DlfileS-'TEST" 

6732 

CREATE  BDAT  Dlfile$,60 

6736 

ASSIGN  gFilel  TO  DlfileS 

6740 

ENTER  »File2;Nrun,Date$,Ldtc< •  ) . 1 1 1 , Bop ,Nht ,Na tp , Nrt ,Corr 

6744 

OUTPUT  OF i lei ; Nrun .DateS ,Ldtc( •  ), Itt.Bop.Nht ,Natp,Nrt ,Corr 

6748 

FOP  1-1  TO  Nrun 

6752 

ENTER  8File2.ToldS,Emf < •  ),Vol t< •  >,Afip< •  ) 

6756 

IF  1-1  THEN  E764 

6760 

OUTPUT  PFilel  . Bop , Tol dS ,Emf  (  •  >,Uol  t(  •  ),Ai«ip<  •  ) 

6764 

NEXT  I 

6768 

ASSIGN  gFale2  TO  • 

6772 

ASSIGN  gFilel  TO  • 

E776I 

RENAME  "TEST"  TO  D2_fileS 

6780 

BEEP  2000,  .2 

6784 

BEEF  4000,  .2 

6788 

BEEP  4000,  .2 

E792 

SUBEND 

6796  1 

6800' 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

6804  1 

6808 

SUE  Move 

68121 

FILE  NAME-  MOUE 

68161 

6620 

1 

DIM  A (66  ),B(66  \C<66 . ,D<66  ;,E(6E  ),F(G6  ),G<66  ),H<66  ) ,  J ( 56  :.K(66  ),L<6 

1 

6824 

DIM  N(66  ),Emf ( 34  J.Uolt ( 2 >,Amp< 1 1  >,Ldtc<  4  ) 

6828 

BEEF 

6832 

INPUT  "OLD  FILE  TO  MOiJE",C:  f;lel 

683E 

ASSIGN  gFile:  TC  D2_f:iei 

6840 

ENTER  BFile2;Nrur,,ToldS,Ldlc<  •  -  , 1 1 1 . Bop , Nht , Netc . Nrt , Cc-r 

6844 

FOP  1-1  TO  Nrur 

6B48 

ENTER  PFile2.ToldS 

B852 

I  ;,N( 

6856 

ENTER  SFileliA!  I  >,B(  I  ),C  I  >,D<  I  >,E<  I  >,F<  I  ),G<  I  >,H(  I  !,J'  I  .< ,  h  ■  I  !,l 
I  ) 

ENTER  §Fi le2 .Enf ( ♦  :  , Ucl  t ' •  >,Amp< •  ) 

6860 

NEXT  1 

,M(66 


l:  :  ',n< 
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ES6C 

ASSIGN  eFiie2  TC  • 

6  656 

BEEC 

BE": 

INPUT  "SHI^T  DI Sh  AND  HIT  CONT 1 NUE " , 0» 

E87S 

EEEC 

688C 

INPUT  "INPUT  BOAT  SIZE'. Size 

6864 

CREATE  60AT  D2_f i let , Si ze 

6886 

ASSIGN  e-ilel  TO  D2_fi leS 

E89Z 

OUTPUT  @Fi lei .Nrun. DateS. Ldtc< • ). Itt ,6op,Nht , Nate, Nrt , Corr 

6896 

FOP  1=1  TC  Nrun 

£900 

OUTPUT  g^ilel ,TcldS 

6904 

OUTPUT  SFileliAd  ),B(I  ),C(I  >,D<I),E(I  ),F<I  ),6(1  ),H(I  ),J(I),K< 

1),L<  I  ),M 

(  I  ),N 

(I  ) 

6908 

OUTPUT  eFilel;Emf(«),Volt<«)  ,  Amp<  •  ) 

6912 

NEXT  I 

6916 

ASSI6N  gFilel  TO  • 

6920' 

RENAME  "TEST"  TO  D2_file$ 

6924 

BEEF  2000, .2 

6926 

BEEP  4000,  .2 

6932 

BEEP  4000,. 2 

6936 

PRINT  'DATA  FILE  MOVED" 

6940 

SUBEND 

6944' 

6948' 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

6952' 

6956 

SUB  PurQ 

6960 

BEEP 

6964 

INPUT  "ENTER  FILE  NAME  TO  BE  DELETED" , Fi leS 

6968 

PURGE  FileS 

6972 

GOTO  6960 

6976 

SUBEND 

6980' 

6  984 IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

6988' 


6992 

6996' 

7000' 

7004 

7008 

7012 

7016 

7020 

7024 

7028 

7032 

7036 

7040 

7044 

7046 

7052 

705E 

7060' 

7064 

70S8 

7072 

7076 

70S0 

7064 

7068 

7092 

70  96 


SUB  Conb 
FILE  NAME 


COMB 


54  ),Volt<2  \Anp(  11  ),Ldtc<  4  ) 


FIXUP-,D2_fileS 
f  lies 


DIM   EmfC 
BEEF 

INPUT    "OLD   FILE    TC 

ASSIGN   @File2    TO   D2_ 

Dl_f  ileS«="TEST" 

CREATE  BDAT  Dl_fileS,30 

ASSIGN  tFilel  TO  Dl_file$ 

ENTER  PFile2 .Nrun, DateS ,Ldtc(«  > , I i t  , Bop , Nht , Na tp , Nr t ,Corr 

IF  K-0  THEN  OUTPUT  SFi 1  el . Nrun , DateS , Ldtc ( •  ) , 1 1 1 ,Bop .Nht , Nat p ,Nr1 

FOR  1-1  TO  N'un 

Vol t< •  > , Amp( ♦  ) 

, Vol t ( •  ) , Amp ( ♦  ) 


,  Corr 


eF;le2.Bop,ToldS ,Enf ( •  !, 
r  gFi lei ;Bop,Tc!dI,Emf <  • 


ENTER 

OUTPU 

NEXT  I 

ASSIGN  BFile2  TO  « 

RENAME  "TEST"  TO  D2_fileS 

BEEP  4000. .2 

BEEP 

01-  e  =  i 

INPUT  "UANT  TO  ADD  ANOTHER  FILE  ( 1= Y , 0=Ni de f au 

IF  01 e=l  THEN 

K"l 

BEEP 

INPUT  "C-IVE  NEU  FILE  NAME" ,Nf ileS 

ASSIGN  &ci it2  TC  NfileS 
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71CC  GOTO  7032 

7104  END  IF 

7105  ASSIGN  PFi leZ  TO  • 
7112  SUBEND 

7116' 

7120IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXVXXXXXXXXX 

71241 

712e  SUE  Readplot 

7132  DIM  Qdp<5 ),Htube(S ),Thetabi5 > 

713E  PRINTER  IS  701 

7140  INPUT  'ENTER  FILE  NAME'.FileS 

7144  INPUT  "ENTER  THE  NUMBER  OF  DATA  PAIRS", Nmn 

7146  ASSI6N  SFilel  TO  File$ 

7152  FOR  1-1  TO  Nrun 

715G      ENTER  PFi le 1 . Odp< •  ) , Htubef •  ) , Thet ab< •  ) 

71G0      PRINT  Qdp(«  ) 

71G4      PRINT 

7166      PRINT  Htube( • > 

7172      PRINT 

7176      PRINT  Thetab<« ) 

7180      PRINT 

7184      PRINT 

7188  NEXT  I 

7192  SUBENO 
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